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ABSTRACT 



Aims. To study the structure of nearby (< 500 pc) dense starless and star-forming cores with the particular goal to identify and 
understand evolutionary trends in core properties, and to explore the nature of Very Low Luminosity Objects (< 0.1 L Q ; VeLLOs). 
Methods. Using the MAMBO bolometer array, we create maps unusually sensitive to faint (few mjy per beam) extended (« 5') 
thermal dust continuum emission at 1.2 mm wavelength. Complementary information on embedded stars is obtained from Spitzer, 
IRAS, and 2MASS. [The journal article will come with FITS files of the maps.] 

Results. Our maps are very rich in structure, and we characterize extended emission features ("subcores") and compact intensity 
peaks in our data separately to pay attention to this complexity. We derive, e.g., sizes, masses, and aspect ratios for the subcores, as 
well as column densities and related properties for the peaks. Combination with archival infrared data then enables the derivation of 
bolometric luminosities and temperatures, as well as envelope masses, for the young embedded stars. 

Conclusions. Starless and star-forming cores occupy the same parameter space in many core properties; a picture of dense core 
evolution in which any dense core begins to actively form stars once it exceeds some fixed limit in, e.g., mass, density, or both, 
is inconsistent with our data. A concept of necessary conditions for star formation appears to provide a better description: dense 
cores fulfilling certain conditions can form stars, but they do not need to, respectively have not done so yet. Comparison of various 
evolutionary indicators for young stellar objects in our sample (e.g., bolometric temperatures) reveals inconsistencies between some 
of them, possibly suggesting a revision of some of these indicators. Finally, we challenge the notion that VeLLOs form in cores not 
expected to actively form stars, and we present a first systematic study revealing evidence for structural differences between starless 
and candidate VeLLO cores. 

Key words. Stars: formation - ISM: evolution - ISM: structure - ISM: dust - ISM: clouds 



1. Introduction 

Stars form from dense gas (» 10 5 cm -3 ). Such gas is, e.g., 
found in discrete nearby (< 500 pc) cold (« 10 K ) small-scale 
(< 0. 1 pc) condensations referred to as dense cores (My ers et alJ 
1983). These are thought to be the sites where low-mass stars 
form dMvers & Bensonlll983l) . Their properties provide the ini- 
tial conditions for star formation. It is thus necessary to under- 
stand the physical state of dense cores in order to be able to un- 
derstand the star formation process in detail. 

Some of the many physical parameters of dense cores are 
their sizes, masses, and densities. Given that thermal dust con- 
tinuum emission at a; 1 mm wavel ength is a better m ass tracer 
than molecular emission lines (e.g. jTafalla et al.ll2.Q02l) , and that 
modern bolometer cameras enable complete maps of the dust 
emission from dense cores to be made within a short time, dust 
emission maps are a prime tool to study the above core param- 
eters. Dust emission surveys thus enable a systematic overview 
of dense core properties. 

In recent years several surveys studied dense cores in 
large (up to several 10 pc) and massi ve (up to s everal 
1O 3 M ) complexes of molecular clouds (M otte et al.l [1998; 



Johnstone & Ballvl 120061: [Johnstone et all 12004 121 
2000; Hatched et al.ll2005b Enoch et al.ll2007ll2006tfi 



2001, 200' 
Young et al 
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* Based on observations with the IRAM 30m- telescope. 



2006b; Stanke et al.l l2006l). Most cores in these clouds are in 
groups and clusters (e.g.. Enoch et al.ll2007l) . Although it is be- 
lieved that most stars in the galaxy form in such environments 
(Mag nani et al.l[T995h . two issues make it difficult to study the 
fundamental physics of star formation in such environments. 
First, cores can be confused in crowded regions forming clus- 
ters; it bec omes difficu l t to id entify and study individual cores 
(see, e.g., iMotte et all [l998; see also IWard-Thompson et al.l 
|2007| for illustrative examples). Second, cores might interact 
in crowded regions; more parameters than relevant in isolated 
cores influence the evolution of dense cores. These problems 
can be avoided by studying dense cores in r elative isolation 
(IClemens & Barvainis 1988; Bourke et al.ll9 95) instead of cores 
hosted in larger (up to > 100 pc) complexes of giant molecular 
clouds (lBlitzlll993h . 



The present study therefore focuses on isolated cores. We 
have observed both starless cores and cores with embedded stars, 
covering a large range in size, mass, central density and (sup- 
posedly) evolutionary status. This creates a unique database for 
detailed studies of dependencies between various dense core pa- 
rameters and their relation to dense core evolution. 
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1.1. The Need for a new Dust Emission Survey 

Several surveys of isolated cores with and without active 
star formation have been conducted in the past » 15 years 
(e.g.. iMotte & Andrell200U IVisser et al.l l200l iKirk et alJ f2005. 



lYoung et al.ll2006al IWu et al.ll2007l) . Thev have covered a large 
number of dense cores in a broad variety of physical states. 
While these surveys in principle should yield a comprehensive 
database on the properties of dense cores, this is in fact true only 
to a limited extent. Many existing studies were relatively insen- 
sitive to faint emission (because of noise) and structure on scales 
beyond 1' to 2' (due to artifacts, e.g., due to 'jiggle maps'). 

State-of-the-art bolometer cameras nowadays allow for the 
reliable imaging of faint large-scale (> 5') extensions of sources 
initially — if at all — only detected in their brightest intensity 
peaks. This urged us to conduct a new comprehensive dust emis- 
sion survey towards isolated dense cores at a; 1 mm wavelength 
using up-to-date instrumentation. For this we u se the Max- 
Planc k Millimetre Bolometer (MAMBO) arrays dKreysa et alJ 
1999) at the IRAM 30m-telescope. This combination of tele- 
scope and receiver was probably the most sensitive facility for 
mapping of extended dust emission available before the recent 
emergence of LABOCA (when considering the levels of ex- 
tended artifacts and statistical noise). The research presented 
here is the first census with MAMBO that covers > 10 starless 
dense cores. 

A further motivation of our survey is the demand for com- 
plementary data on dense cores for the recent Spitzer Space 
Telescope imaging surveys of dense cores. Such data is needed 
to, e.g., put the properties of the stellar content probed by Spitzer 
in context with dense core properties like density and mass; the 
latter can usually not be derived from the Spitzer data. In particu- 
lar the Spitzer Legacy Program "From Molecular Cor es to Planet 
Form ing Disks" (AKA. "Cores to Disks" or "c2d"; Evan s et all 
[2003b . which imaged 82 nearby (< 500 pc) isolated dense cores 
as part of its agenda, stimulated a number of coordinated surveys 
of the c2d cores covering many spectral windows. Our survey — 
the c2d MAMBO survey — is one of four new dust continuum 
emission s urveys of isolated cores carried out in the c2d frame- 
work (with lYoung et al.ll2006al [\Vu et alj|2007l and Brede et al., 
in prep., being the other three). 



1.2. Specific Aim of the c2d MAMBO Survey 

Dust emission data can be used to address a number of open 
issues in star formation research. In general, stars are thought 
to form from very dense molecular cores (H2 densities > 
10 5 cm -3 ). In particular, the spatial mass distribution is thought 
to govern t he stability of dense cores against gravitational col- 
lapse (e.g. jMcKee & HollimaD| [T999). One might thus hope to 
derive criteria to be fulfilled to allow star formation in dense 
cores by comparing the mass distributions of starless cores with 
those of cores actively forming stars. This could lead to a better 
understanding of the underlying star formation physics. The c2d 
MAMBO survey is ideally suited for such comparative studies 
since it probes the mass distribution in cores with and without 
ongoing star formation at high sensitivity (see Sec. |5.3.4l for de- 
tails). 

Also, the evolution of young stellar objects (or "protostars"; 
hereafter YSOs), i.e., stars that are surrounded by significant 
amounts of the matter from which they did or do form, can be 
studied with dust continuum emission data. These allow the in- 
vestigation of the relation between the broadband spectral prop- 



erties, the structure of the circumstellar envelope, and present 
models of evolutionary schemes (see Sec. |5.4| for details). 

Our sample includes four candidate Very Low 
Lu minosity Objects (VeLL Os; see iKauffmann et al.l I2005I 
and Idi Francesco et al.l I2007I) . enigmatic objects of apparent 
protostellar nature but unu sually low luminos i ty that were 
recen t ly discovered by c2d dYoung et alJ l2004al iBourke et al.l 
I2005L IDunham et all l2006h cf. [Andre et aL| [1999J). They are 
defined as objects with internal luminosities (i.e., not including 
the dense core luminosity from interstellar heating due to, e.g., 
the interstellar radiation field) < 0. 1 L Q that are embedded in 
dense cores. Some of the se objects have be en interpreted as 
very young (< 10 4 years: lAndre et al] 1 1999b YSOs of stellar 
final mass (i.e., > 0.08 M Q ). Based on their low accretion rates, 
others are thought to be young objects in the process of growing 
to su bstellar final mass dYoung et al.ll2004allBourke et alj |2005; 
see IDunham et al.l 120061 for a summarizing discussion). The 
c2d MAMBO survey is the first study with a sensitive and 
homogeneous data set on several starless, YSO, and VeLLO 
natal cores. This allows comparative studies between these core 
types to better understand how VeLLO natal cores differ from 
others (see Sec. |5.5| for details). 



1.3. Structure of the present Study 

In Sec. [2] we begin with an introduction to our sample, provide 
a summary of the observations, and detail our data reduction. 
This includes a description of a new iterative data reduction 
scheme designed to mitigate problems affecting the reconstruc- 
tion of maps of weak extended emission. Section [3] provides a 
first discussion of the source properties and an overview of the 
source identification and quantification strategies. YSO prop- 
erties from supplemental data are derived in Sec. [4] The main 
discussion of the data follows in Sec. [5] where we in particular 
study the structure and evolution of dense cores in general (Sec. 
I53T) . YSO cores (Sec. E3), and VeLLOs (Sec. |53). We sum- 
marise our study in Sec. [6] Appendix lAl discusses the standard 
dust emission properties adopted by the c2d collaboration to de- 
rive the mass distribution from observed dust emission maps. 



2. Sample, Observations, and Data Reduction 

2. 1 . Sample Selection 

Our sample of 38 cores (TableQ~|i is drawn from the initial target 
list of the c2d Spitzer survey of isolated dense cores dEvans et al.l 
2003). This list contains dense cores within about 400 pc from 
the sun that are smaller than w 5' and were mapped in dense gas 
tracers before the start of the c2d surveys. These cores are nearby 
enough to allow for the detection of intrinsically faint embedded 
sources. Their small apparent size allows them to be mapped 
in a reasonable tim e. This sample was drawn from the compila- 
tions a nd surveys of Jiii na et ail d 19991 forNH3 maps ), Lee et al.l 



dljuTl for CS and N 2 H + maps), and ICaselh et all (|2002|, for 
N2I-L maps), and completed by adding individual cores. Due to 
time constraints not all sources in our sample were finally ob- 
served by Spitzer as part of c2d. Some further cores have or will 
be covered as part of the Spitzer Gould Belt survejQ, or have 
been covered in other programs (see TableQ]for details). 

The initial c2d target list contains about 150 cores. At 
the start of the c2d survey, MAMBO was the most efficient 
bolometer camera available for sensitive mapping of large fields. 



http ://w w w. cf a. harvard, edu/gouldbelt/index .html 
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Therefore, those cores in this list that were expected or known 
to have faint or extended dust continuum emission, or both, and 
are visible from Spain, were mapped by MAMBO. The other 
cores were observed in other surveys with SCU BA on the JCMT 
(lYoungetal.ll2006ah . SHARC n on the CSO dWu et al.ll2007l) . 
and with SIMBA on the SEST (Brede et al., in prep.). We ob- 
tained usable MAMBO data for 38 cores, which are listed in 
TableQ] Of these cores 21 are starless to our present knowledge, 
4 contain VeLLO candidates (out of which 2, those in L1014 
and L1521F, can be considered proven VeLLOs), and 13 contain 
YSOs (see Sec. 14. II for a discussion of the association between 
YSOs and dense cores). Spitzer data is available for 32 cores 
(i.e., 3/4) of the MAMBO sample (see Table Q] for Spitzer ob- 
servation references). The MAMBO sample covers 25 out of the 
82 cores observed by c2d, and thus roughly 1/3 of the total c2d 
sample. In total 13 out of the 17 YSO and VeLLO candidate 
cores observed by us using MAMBO were targeted by Spitzer. 
However, only for 8 of these (less than half) there is c2d data and 
therefore homogeneous Spitzer photometry. 

The observed cores reside in very different environments, 
ranging from regions of isolated low-mass star formation, like 
Taurus, over more turbulent regions, like the Cyg OB7 molecu- 
lar cloud complex and the Cepheus Flare, to the high-mass star- 
forming site of Orion. In this respect the c2d MAMBO sample 
does not constitute a homogeneous sample, but a cross section 
of the different star-forming clouds in the solar neighborhood. 
Correspondingly our sample includes cores from near (x: 100 pc) 
to intermediate (« 400 pc) distances. 

2.2. Observations 

Continuum observations of the 1.2 mm thermal dust emission 
were done with the IRAM 30m-telescope on Pico Veleta (Spain) 
using the 37-receiver MAMBO- 1 (projected array diameter of 
120") and 117-receiver MAMBO-2 bolometer (240" diame- 
ter) cameras of the Max-Planck-Institut fur Radioastronomie 
dKreysa et al] fl999). A total of 122 usable maps was taken be- 
tween the summer of 2002 and the winter of 2003/2004 in the 
framework of a flexible observing pool (Table [2]). Only 6 maps 
were acquired using the MAMBO- 1 array. The weather condi- 
tions were good, with zenith optical depths between 0.1 and 0.3 
for most of the time, and above this for * 20% of the time. Some 
of the maps were affected by clouds. All but 13 maps (i.e., 11%) 
were taken at an elevation above 40°. This yields typical line-of- 
sight optical depths much bel ow 0.5. Two maps ar e affected by 
strong anomalous refraction dAltenhoffet al.l l 1987b one map on 
the southern part of L1041-2, one map on LI 148). 

The beam size on the sky was 11", and the effective fre- 
quency 250 GHz with half sensitivity limits at 210 and 290 GHz. 
Pointing and focus position were usually checked before and af- 
ter each map. The pointing and focus corrections were usually 
below 3" and 0.3 mm, respectively. The zenith optical depth 
was typically measured with a skydip at least once within an 
hour from start or end of a map. The sources were observed 
with the standard on-the-fly technique, where the telescope sec- 
ondary was chopping by * 40" to * 70" at a rate of 2 Hz par- 
allel to the scanning direction of the telescope. The telescope 
was constantly scanning in azimuth for up to 90 s at a speed 
of 6 to 8" s before turning around, except for one map of 
L1521B-2 for which the telescope moved in azimuth and eleva- 
tion (then using the "rotated" chopping secondary). Where pos- 
sible, each source was mapped with varying scanning directions 
(in equatorial coordinates) and chop throws. Several times per 
week Mars and Uranus were observed for absolute flux calibra- 



Table 2. MAMBO map details. For each core the products 'num- 
ber of maps' x 'array type' list how often a core was observed 
with which array ('37' for MAMBO- 1 and ' 1 17' for MAMBO- 
2). Several products are given if a dense core was observed with 
both arrays. The last column lists the minimum noise level in a 
map after smoothing to 20" resolution. 



c2d Name 


MAMBO Maps Minimum Noise Level 
mjy per 11" beam 


LI JJJ 


2 x 117 


1.6 


T 1 S91R-9 


5 x 117 


1.1 


L-iL jz ir 


2 x 117 


2.0 


T 1 S9 1 9 


4x 117 


1.3 


T 1 S9A A 


2 x 117 


1.4 


R 1 & 1 
IJ 1 O 1 


4x 117 


1.1 


TMC-2 


4x 117 


0.9 


B18-2 


2 x 117 


1.3 


RIQ 4 
lJ 1 O "H* 


3 x 117" 


1.1 


R 1 R S 


3 x 117" 


1.1 


TVIP 1 P 


6 x 117 


1.2 


TA4P 1 


6 x 117 


1.7 


L1JU / t\ 


5 x 117 


1.6 


T 1 SR9 A 


3 x 117 


1.3 


Hvrt.OuJt 1 J 


3 x 117 


1.2 


T 169? A 


7 x 117 


1.1 


T 


4 x 117 


1.4 


T zHR 

1_#tO O 


2 x 117 


1.2 


T zLQ9 


3 x 117 


0.8 


PR 1 

l_-I3 loo 


1 x 37,2 x 117 0.8 


LD / J- 1 


2 x 117 


1.3 


T fV7S 

LU / J 


1 x 37,2 x 117 0.9 


T 1 1 DO 


1x37,2x117 1.3 


T 10/11 9 


4 x 117 


1.4 


LI 148 


5 x 117 


1.1 


L1155E 


2 x 117 


1.5 


L1082C-2 


2 x 117 


1.8 


L1082C 


2 x 117 


1.5 


L1082A 


4 x 117 


1.1 


L1228 


3 x 117 


1.3 


Bern48 


1x37,1x117 1.0 


L1172A 


2 x 117 


1.5 


L1177 


1 x 117 


1.3 


L1021 


3 x 117 


1.0 


L1014 


2 x 37 


1.3 


LI 103-2 


3 x 117 


1.4 


L1251A 


3 x 117 


1.2 


LI 197 


6x 117 


0.9 



Notes: a) B18-4 and B18-5 were observed together in one field. 



tion. "Secondary" calibrators of constant but not a priori known 
flux density were observed every few hours. The scatter in their 
retrieved flux densities suggests a relative map-to-map calibra- 
tion uncertainty < 15%, suggesting an absolute uncertainty of 
order 20%. 

Pa rt of the dat a were already pre sented by Crapsi ^taT] 
(120041. forL 1521FllCrapsietaT1 J200ll for TMC-2 and L4 92), 
Youn g et all d2004al for L1014). and lSchnee et all d2007bl for 
TMC-1C). The raw data for these sources are included in this 
work and are partially complemented by additional mapfl 



Some of these data was kindly provided by A. Crapsi. 
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Table 1. Dense cores covered by the c2d MAMBO survey. The table lists, from left to right, the core name adopted within the c2d 
collaboration, names for the cores adopted by other authors, the core distance, the type of the core (i.e., starless, YSO, or VeLLO), 
the approximate galactic environment of the dense core, and whether Spitzer data is available for the core. The last column gives 
a rough description of the morph ology of the sub cores (i.e., extended substructures; see Sec. 13. 2\ within each core using the c2d 
standard morphology keywords (lEnoch et al.ll2.006l: see Sec. l3.3l for a description). 



c2d Name 


Other Names 


Distance" 
pc 


lype 


Region' 


Spitzer Data 


Morphology Keywords'" 


L1355 




200 + 


50 1 


s.l. 


Cepheus Flare 


c2d 


elongated, extended 


L1521B-2 


B217B 


140 ± 


10 2 


s.l. 


Taurus 


c2d 


multiple, elongated, extended 


L1521F 




140 ± 


10 2 


Vc. 


Taurus 


c2d 


elongated, extended 


L1521-2 




140 ± 


10 2 


s.l. 


Taurus 


c2d 


round, extended 


L1524-4 




140 ± 


10 2 


s.l. 


Taurus 


c2d 


multiple, round to elongated, extended 


B18-1 


TMC-2A 


140 ± 


10 2 


s.l. 


Taurus 


c2d 


multiple, round to elongated, extended 


TMC-2 


LI 529 


140 ± 


10 2 


s.l. 


Taurus 


c2d 


multiple, elongated, extended 


B18-2 




140 ± 


10 2 


s.l. 


Taurus 


c2d 


multiple, elongated, extended 


B18-4 


TMC-3 


140 ± 


10 2 


YSO 


Taurus 


IRAC, MIPS 


multiple, elongated, extended 


B18-5 




140 ± 


10 2 


s.l. 


Taurus 


c2d 


multiple, elongated, extended 


TMC-1C 


B220 


140 ± 


10 2 


s.l. 


Taurus 


c2d 


multiple, elongated, extended 


TMC-1 


L1534 


140 ± 


10 2 


s.l. 


Taurus 


c2d 


multiple, elongated, extended 


LI 507 A 




140 ± 


10 2 


s.l. 


Taurus 


c2d 


multiple, round to elongated, extended 


L1582A 


B32 


400 ± 


40 3 


YSO 


A Orionis 


IRAC, MIPS 


multiple, elongated, extended 


IRAS05413 


HH212-mm 


450 ± 


50 4 


YSO 


Orion B 


IRAC, MIPS 


multiple, round to elongated, extended 


LI 622 A 




450 ± 


50 4 


s.l. 


Orion East 


IRAC, MIPS 


multiple, elongated, extended 


L183 


L134N 


110 ± 


10 5 


s.l. 


north of Ophiuchus 


IRAC, MIPS 


multiple, round to elongated, extended 


L438 


CB119 


270 ± 


50 6 


s.l. 


Aquila Rift 


c2d 


round, extended 


L492 


CB128 


270 ± 


50 6 


s.l. 


Aquila Rift 


c2d 


elongated, extended 


CB188 


L673-1 


300 ± 


100 7 


YSO 


Aquila Rift, Cloud B 


c2d 


elongated, extended 


L673-7 




300 ± 


100 7 


Ve. 


Aquila Rift, Cloud B 


c2d 


multiple, round to elongated, extended 


L675 


CB193 


300 ± 


100 7 


s.l. 


Aquila Rift, Cloud B 


c2d 


elongated, weak 


L1100 


CB224 


400 ± 


50 s 


YSO 


northern Cyg OB7 




elongated, extended 


L1041-2 




400 ± 


50 s 


YSO 


northern Cyg OB7 




multiple, round to elongated, extended 


L1148 


LI 147 


325 ± 


25 9 


Ve. 


Cepheus Flare 


c2d 


multiple, elongated, extended 


L1155E 


L1158 


325 ± 


25" 


s.l. 


Cepheus Flare 


c2d 


multiple, elongated, extended 


L1082C-2 




400 ± 


50 10 


s.l. 


Cyg OB7 




elongated, extended 


L1082C 




400 ± 


50 10 


YSO 


Cyg OB7 




elongated, extended 


LI 082 A 


B150 


400 ± 


50 10 


YSO 


Cyg OB7 




multiple, round to elongated, extended 


LI 228 




200 + 


50 11 


YSO 


Cepheus Flare 


c2d 


round, extended 


Bern48 


RN0129 


200 ± 


50" 


YSO 


Cepheus Flare 


c2d 


round, extended 


L1172A 




288 ± 


25 9 


YSO 


Cepheus Flare 


IRAC, MIPS 


multiple, round, extended 


L1177 


CB230 


288 ± 


25 9 


YSO 


Cepheus Flare 


IRAC 


round, extended 


L1021 




250 ± 


50 12 


s.l. 


Cyg OB7 


c2d 


elongated, extended 


L1014 




250 ± 


50 12 


Vc. 


Cyg OB7 


c2d 


elongated, extended 


LI 103-2 


L1106 


250 ± 


50 12 


s.l. 


Cyg OB7 




elongated, extended 


L1251A 




300 ± 


50 13 


YSO 


Cepheus Flare 


c2d 


multiple, round to elongated, extended 


LI 197 




300 ± 


50 14 


s.l. 


Cepheus Flare 


c2d 


elongated, extended 



Notes: 

a) Standard c2d distances (Bourke, priv. co mm.), derived by red dening methods or ass ociation (in position- veloc ity space) to objects of k nown 
distance. The ind i vidual references are: [1] Obavash i et alj d 19981). vi a ass ociation; [ 2] Kenyon et al.l dl994h: r3llMurdin& Penston dl977l) : [4] 
iGenzel & Stutzkil dl989h , via a ssociation to L1630: 151 iFrancol dl989l): T61 IStraizvs et alJ d2003l) : 171 iHerbig & Jonesl dl983l) . via association to 
L673 ; [8] Dobashi et al. ( 1994), and associa tion to L1041; r9l|Straizvs et al] d!992l) . and association with LI 172 ; [10] via reddening , following 
Bour ke et al] dl995h: Ull iKunl dl998h : H21 iPagani & Breart de Boisanged dl996h . via association to B164; [13] iKun & Prustil dl993l) : and [14] 
lYonekura et al.ldl997h . via association to LI 192/L1200. 

b) The abbreviations mark starless cores ('s.l.'; these have no known star within the core boundaries defined in Sec. 13. 2\ , YSO cores (YSO'; YSOs 
within the core boundaries), and VeLLO cores ('Ve.'; VeLLOs within the core boundaries). The cores L1521-2, B18-1, TMC-1, and L1622A are 
excluded from the YSO category since the YSOs are off the main core body. See Sec. 14. ll for a discussion of the association between YSOs and 
dense cores. 

c) The designation of regions roug hly f ollows iDame et alj dl987l) . lMaddalena et al] d!986h . lYonekura et al] d 19971) . and lDobashi et alj d 19941) . 

d) Cores with data from the c2d project dEvans et al.l2003l) are labeled 'c2d'. For other cores existing Spitzer imaging data is indicated in the table, 
and the Spitzer program numbers and Pis are: GO 3584 (Padgett) and GO 20302 (Andre) for B18-4; GO 20339 (Stauffer) for L1582A; GO 3315 
(Noriega-Crespo) and GTO 47 (Fazio) for IRAS05413; GTO 47 (Fazio) for L1622A; GTO 94 (Lawrence) and GTO 53 (Rieke) for L183; L 30574 
(Allen) and IOC 717 (Rieke) for LI 172A; and GTO 124 (Gehrz) for LI 177 

e) The c2d standard keywords are adopted from the c2d Bolocam survey by Enoch et al. ( 2006). See Sec. l3.3l for a description. 
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2.3. Data Reduction 

2.3.1 . Classical Limitations & A new Strategy 

The reduction of bolometer data for extended sources is still an 
evolving art, in particular when dealing with faint emission. The 
basic problem is that the brightness of the sky usually exceeds 
the intensity of the astronomical sources by 2 to 4 orders of mag- 
nitude (e.g., at 1.2 mm wavelength the intensity from a 300 K 
atmosphere of optical depth > 0.1 exceeds the intensity from in- 
terstellar dust at 10 K with an IT column density < 10 23 cm -2 by 
a factor > 640). The sky brightness (and if possible also its rapid 
fluctuations known as "skynoise", that contribute to the noise) 
thus has to be subtracted during the data reduction. This subtrac- 
tion unfortunately can introduce artifacts to the derived maps. 

Specifically, classical data reduction algorithms filter a sig- 
nificant fraction of the source flux when sources have a size com- 
parable to the size of the bolometer array (relevant when sub- 
tracting skynoise), a size similar to an individual bolometer map 
(i.e., before mosaicing; relevant when subtracting the mean sky 
brightness as, e.g, the minimum intensity in the map), or both. 
(Note that these fundamental problems also affect maps from 
modern total-power bolometers like LABOCA0, SCUBA-20, 
and BolocarrQ A d etailed analy s is of t his problem for Bolocam 
has been given by Eno ch et all (120061) .) This is the case for a 
large fraction of the sources studied in our sample. We there- 
fore implement an innovative data reduction scheme that par- 
tially mitigates classical reconstruction problems. 

The central element of our data reduction algorithm is the 
subtraction of the expected signal due to the astronomical source 
(based on a source model) from the raw data before determin- 
ing and then subtracting from it the mean sky brightness and 
skynoise. This processed raw data is then turned into a map. The 
sum of this "residual" map (since source emission is subtracted) 
and the source model yields an improved source model that can 
be used in further iterations. This approach has the advantage 
that the resulting source model suffers from lesser artifacts due 
to skynoise and sky subtraction than if no existing knowledge 
of the source structure was used during the subtraction. The ini- 
tial source model for the iterative scheme comes from a map 
created using classical (i.e., biased) algorithms. It can be shown 
analytically that, in the ideal case without noise, sources much 
larger than in the classical approach can be recovered. In our iter- 
ative scheme, the residual intensities decrease with a factor 2~ n , 
where n is the number of iterations. To be specific, full recovery 
of the source structure requires the source diameter, Source, to be 
smaller than the effective map size, 

^source ^ ^scan ^chop > ( 1 ) 

where € scm and 4hop are (following IRAM terminology) the scan 
length and the chop throw (of the secondary mirror), respec- 
tively. Slightly larger sources can still be recovered, but conver- 
gence is slower than 2~ n . This size can be increased by mosaic- 
ing. Classical methods, on the contrary, require for full recovery 
(with ff array being the array diameter) 

^source ^ ^scan ^chop — ^array (2) 

when not subtracting skynoise, and have the additional con- 
straint 

^source ^ ^array (3) 

3 http://www.apex-telescope.org/bolometer/laboca/ 

4 http://www.jach.hawaii.edu/JCMT/continuum/scuba2.html 

5 http://casa.colorado.edu/ aikinr/research/bolocam.html 



if doing so. Mosaicing will furthermore bias intensities. 
Experiments with simulated raw data including noise show that 
the simulated maps are consistent with the above ideal analysis. 

In a formal sense, the improvements due to our approach thus 
mean that observations are only limited by Eq. (Q}, while con- 
ventional approaches are limited by Eqs. (0 and (0. The gain in 
mapping parameter space is illustrated in Fig. IB. II in Appendix 
IB1 There, we also present a direct comparison between maps de- 
rived using conventional methods to those from our approach 
(Fig.EH). 

For typical MAMBO map sizes i sc . dn = 360" to 540" and 
chop throws £ c hop = 60" sources with diameters Source < 6' to 9' 
can thus be recovered without significant bias when using our 
new approach. For comparison, when using the larger MAMBO- 
2 array, classical algorithms would only recover sources with 
a diameter <i SO urce < 3' to 6', where skynoise filtering sets the 
additional constraint that d sous:ce < 4' even for much longer scan 
lengths. 

2.3.2. Implementation 

Our iterative data reduction scheme is implemented using the 
MOPSI software package, which is developed and maintained 
by R. Zylka (IRAM Grenoble); for L1521B-2 we used MOPSIC, 
an upgrade of MOPSI also able to handle data obtained us- 
ing the "rotated" wobbler. These packages provide subroutines 
to estimate and remove the correlated skynoise, to reconstruct 
maps from chopped bolometer data using the EKH algorithm 
(Eme rson et alj|l979l) . to simulate chopped intensity difference 
maps for given source models, and many other tools. These 
routines are combined to carry out our iterative data reduction 
scheme. 

During the first iteration no skynoise is removed. The 
chopped intensity difference maps derived for each epoch are 
co-added and then restored using EKH reconstruction. In the re- 
sulting map the area with emission is then marked by hand. The 
intensity for the area outside the area with emission is set to zero, 
as is done for all pixels with negative intensity inside the area 
with emission. The modified intensity map is then smoothed to 
an effective resolution of 20", except for visually selected areas 
with strong emission, where the data is smoothed to 14" res- 
olution. This smoothing increases the signal-to-noise ratio for 
emission more extended than the beam. This map is then taken 
as the source model for the second iteration. 

In the further iterations the simulated chopped intensity dif- 
ference map for the source model is first subtracted from the raw 
data for every single observation. Then the correlated skynoise 
is estimated and subtracted in these residual maps. We use a 
correlation exclusion radius of 50" for all sources with data 
from the 37 receiver MAMBO- 1 camera (i.e., for each array re- 
ceiver a neighborhood of 50" radius is excluded when estimating 
skynoise, since this reduces intensity biases), and 100" other- 
wise. To remove offsets due to sky intensity and instabilities, a 
constant intensity (i.e., a baseline of order 0) is then subtracted 
from each receiver signal for every scan lag. A co-added inten- 
sity map including all observations is then derived using EKH 
restoration. This map is then added to the source model, and a 
source model for the next iteration is derived as done above. In 
total 10 iterations are performed. Note that the smoothing of the 
source model does not reduce the resolution in the reconstructed 
maps: the source map derived at the end of every iteration is 
corrected for differences between the actual and model source 
structure due to smoothing since these show up in the residual 
maps finally added to the source model. 
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The final maps presented here are usually smoothed to an 
effective resolution of 20", though they are available for resolu- 
tions down to 11". This smoothing improves the signal-to-noise 
level for emission more extended than the beam. In order to fa- 
cilitate comparison with maps of different resolutions the data 
are calibrated in mJy per 11" beam. Further weight maps give 
the variation of the effective integration time per pixel across the 
map. They are calculated during the data reduction and used to 
calculate noise maps. The noise increases with increasing dis- 
tance from the map center due to a decrease in exposure time 
towards the map boundaries. 

3. Results 

3.1. Maps 

Figure Q]presents the intensity maps of the dense cores covered 
by our MAMBO observations. The minimum noise level in the 
maps ranges from 0.8 to 2.0 mJy per 11" beam and is listed for 
each source in Table [2] To facilitate comparison, all maps are 
drawn with the same contouring levels and to the same angular 
scale. 

In most sources the maps reveal extended emission on scales 
of several arc-minutes, and in many fields the presence of com- 
pact emission with a full width at half maximum of only a few 
beams. Compared with many other bolometer surveys that often 
detect only the most prominent intensity peaks, our MAMBO 
maps are very rich in structure. 

In spite of our careful data reduction, some maps are still af- 
fected by artifacts such as stripes in the scanning direction and 
extended areas of negative intensities. Sources significantly af- 
fected by this are given lesser weight in our analysis. These are in 
particular objects with obvious scanning stripes in their bound- 
aries (e.g., L183 C4; see below for object names), i.e., series of 
intensity peaks in the scanning direction separated by the chop 
throw, those located at map boundaries (e.g., L183 C6), or with a 
low intensity contrast to the surrounding at the object boundary 
(e.g., L1082C-2 CI; often due to artifacts). 

[The journal article will be accompanied by the maps in 
FITS format.] 

3.2. Source Identification and Quantification 

The MAMBO maps reveal extended emission and compact 
emission peaks with no direct correspondence between these. 
We therefore separately identify and quantify extended features, 
which we term 'subcores', emission peaks in these subcores, and 
additional 'significant peaks' that fulfill some significance crite- 
rion and are not necessarily located within the subcore bound- 
aries. 

For some properties it is not possible to derive their un- 
certainties by Gaussian error propagation. In these cases we 
run Monte-Carlo simulations with artificial noise with proper- 
ties identical to the actual one to estimate the uncertainties. 
Uncertainties are then derived as the standard deviation of the 
derived values. To reduce the noise level for the extended emis- 
sion, all quantities derived for the cores and peaks are derived 
from maps smoothed to 20" resolution unless noted otherwise. 

3.2.1. Extended Structures 

Identification Our source identification method aims at a formal 
definition of what a human investigator would intuitively iden- 



tify as a source. Our criteria are therefore not mathematically 
exact, but adopted from a series of trial-and-error experiments. 

To identify sources we generate signal-to-noise ratio (SNR) 
maps for a beamsize of 20", which are then smoothed to a 
spatial resolution of 60". Regions above an SNR of 2 in these 
maps are taken as source candidates. Obvious artifacts at the rim 
of the maps are excluded. The remaining area above an SNR 
of 2 is then hand-divided into 'subcores'. For this we search 
the 60" SNR maps for saddle points and draw core boundaries 
by following the gradient. This scheme is conceptually sim- 
ila r to, but not identica l with, the "CLUMPFIND" algorithm 
bv lWilliams et al.l dl994). Our scheme thus assigns all emission 
above an SNR of 2 to some 'subcore'. 

Quantification Table [3] lists general information on the identi- 
fied subcores. For each subcore we list the centroid position, de- 
fined as the mean position of the pixels in each subcore, where 
the pixels are not weighted by their intensity. The quoted uncer- 
tainties reflect the influence of statistical noise, while systematic 
contributions are not folded in. For each field mapped the sub- 
cores are labeled by a 'C, followed by a number in order of 
increasing right ascension. In the tables a star marks those cores 
associated with a YSO (see below). Cores with uncertain prop- 
erties due to an unreliable reconstruction are marked by square 
brackets. 

We then list the area assigned to the subcore, A, and the cor- 
responding effective radius, 

r eff = (A/tt) 1/2 . (4) 

Furthermore, we derive the flux densities, F, and the correspond- 
ing masses, M, for the total subcore's area (subscripts 'tot') and 
for the area above 50% of the peak intensity for the respective 
subcore (subscripts '50%'). The masses are derived from the 
flux densities as described in Sec. I3.2.3I Finally, Table [3] lists 
the maximum SNR for each subcore, Sjjf™. The uncertainties of 
F to t and M tot are derived by Gaussian error propagation of the 
measured noise level. For Fsq% and M$q% the uncertainties are 
estimated by the Monte-Carlo experiments described above. 

Table|4]lists geometrical information on the identified cores. 
This is done by fitting ellipsoids to the area above 50% of the 
peak intensity for the respective core (no weighting by intensity 
is appliedfj- We give the major and minor axis, a and b, their 
ratio, and the position angle (east of north). The filling factor, /, 
gives the fraction of the fitted ellipse filled with emission above 
50% peak intensity. Low filling factors indicate cases in which 
ellipses are a bad fit to the observed intensity distributions. We 
also list the effective diameter for the area above 70% of the 
peak intensity, r-j^% (ICrapsi et al. 2005j). No attempt has been 
made to separately handle subcores in which the emission from 
two or more peaks contributes to rjo%. In these cases 7-70% is 
overestimated by an uncertain amount. 

3.2.2. Peaks 

Identification In order to quantify all intensity peaks in our 
maps we study the peak position of every identified subcore. In 
addition we search each map manually for 'significant peaks'. 
These have an SNR of at least 4 and can be spatially separated 
from already identified peaks. To identify them we search the 
map for local intensity peaks with an SNR > 4 and draw the 



6 This is implemented by using the IDL routine FitEllipse by David 
Fanning. 



Jens Kauffmann et al.: MAMBO Mapping of Spitzer c2d Small Clouds and Cores 



7 



surrounding contour at the peak intensity minus three times the 
local noise level. Such a peak is taken to be significant if no 
brighter peak is situated in this contour. In the area outside iden- 
tified cores we also require that there is also no fainter local peak 
within the contour. This additional restriction excludes most ar- 
tifacts. 



Quantification Peak properties are listed in Table[5] The quoted 
uncertainties reflect the influence of statistical noise, while sys- 
tematic contributions are not folded in. For each mapped field 
peaks are labeled by a 'P', followed by a number in order of 
increasing right ascension. Uncertain peaks, which are not well 
reconstructed, are marked by square brackets. These were iden- 
tified by searching the maps for peaks close to map boundaries, 
and those affected by scanning stripes. 

Besides the peak position we give the subcore to which the 
peak belongs, if applicable. The peak intensity, Fjjjjjj* 1 , is listed 
besides the corresponding peak H2 column density and visual 
extinction, NQI2) and Ay- These are derived using the c2d stan- 
dard assumptions on dust properties, as outlined in Appendix 
lAl but they are only listed for starless cores for which there 
is no internal heating. In YSO cores the dust temperature ex- 
ceeds 10 K and needs to be calculated separately for each source 
(Sec. l3.2.3l l. Also given are the flux density and the correspond- 
ing mass within an aperture of 4200 AU radius, F4200AU and 
M4200AU, as derived in Sec. 13.2.31 This aperture of constant phys- 
ical diameter is a better measure of the mass distribution than the 
central column density since, for close enough sources, beam 
smearing does not p lay a role. The specific radius of 4 200 AU is 
somewhat arbitrarv; lMotte & Andrei d200 if) used it before, given 
that it matches 30" at the distance of Taurus (« 140 pc) and 
is about the radius of an envelope's volume available for YSO 
formation (Sec. 15.3.3b . We finally list the flux density for the 
c2d standard apertu res of 20", 40", 80 ", 120" diameter (e.g., 
lYoung et al] l2006al and IWu et alj [2007b . In deviation to most 
other parts of this work the flux densities for the apertures are 
derived from unsmoothed maps since the smoothed maps some- 
times have a resolution comparable to the aperture size. 



3.2.3. Mass and Column Density Estimates 

Mass and column density estimates for starless cores are derived 
from the observed flux densities assuming the c2d standard dust 
emission properties, i.e., a dust temperature of 10 K and an opac- 
ity of 0.0102 cm 2 g- { (per gram of ISM) at MAMBO's effective 
observing wavelength. Appendix lAl presents a detailed discus- 
sion of the conversion between mass and dust emission. The re- 
lated discussion of uncertainties in temperature and dust opac- 
ities suggest a systematic uncertainty for mass estimates of an 
order of a factor two. Dust near YSOs is, however, heated by the 
YSO radiation, and for later YSO stages likely to be more co- 
agulated than interstellar dust. Therefore, we adopt higher dust 
temperatures and opacities when deriving masses for YSO cores. 
The temperature gradients also prevent us from deriving mean- 
ingful estimates of the peak column density of YSO cores; there- 
fore we do not quote values. 

The dust temperature near YSOs can be estimated from ana- 
lytic models of temperature profiles for YSO envelopes cooling 
via dust emission. For opacities with power-law dependence on 
wavelength, k cc A~@, where f3 is the emissivity spectral index, 



the dust temperature profile in a YSO envelope only heated ra- 
diatively by the star is 



(5) 



where q = 2/(4 + IS), L is the luminosi ty, and r is the 
distance from the star dTer ebev et al] 1 1993b . The numerical 
values hol d for /3 - 1, whic h is consistent with observa- 
tions (e.g.. lBelloche et a"Dl2006h and assumed in the following. 
Rearrangement allows for an estimate of the radius at a particu- 
lar dust temperature, 



^(r d ) = 2890Au(lVL) te) 



(6) 



However, in absence of nearby stars, the interstellar dust is 
heated by the interstellar radiation field and in collisions with 
gas (which is heated by cosmic rays and the UV-part of the in- 
terstellar radiation field), whic h produces a typ ical equilibrium 
temperature of order 10 K (e. g..lGoldsmithl200ll: see lEvans et al] 
1200 ll iGalli et al] [20021 and (Young et al]|2004bl for actual core 
temperature models). The dust temperature cannot drop below 
this value, so that 



r d (r) = max(10 K, Tf[r]) 



(7) 



For this temperature profile iBelloche et al] d2006l) derive the 
mass-weighted dust temperature. This assumes optically thin 
dust emission, which applies to the majority of the mass. Only 
the densest and hottest parts of an envelope might be optically 
thick, which leads to an insignificant bias towards slightly over- 
estimated average dust temperatures. For this calculation we as- 
sume the density profile in a subcore to be roughly described by 
a power-law, i.e. g cc r~ p , where g is the density and r is the 
distance from peak center. Then the mass-weighted dust temper- 
ature within an aperture of radius R is 



(T d ) = 



3-P ^ 
3—q—p ' d 



if T r , dd (R) > 10 K, 



j_2_ [z^SOS J 3 "" + A . io K , otherwise . 



(8) 



For YSO peaks and subcores we use this average temperature to 
derive masses from the flux densities of the whole subcore, of the 
area above 50% peak intensity, and of an aperture of 4 200 AU 
radius. In these cases the aperture radius R is set to be the ef- 
fective core radius of the whole core, or the geometric mean of 
the major and minor axis at 50% peak intensity, or 4200 AU, 
respectively. We assume a power-law exponent of p = 2 fo r the 
densi ty profile (i.e., profile for a singular isothermal sphere: IShul 
1977). The mean dust temperatures are thus overestimated for 
density profiles shallower than g cc r~ 2 . The YSO luminosities 
are taken to be identical to the bolometric luminosities derived in 
Sec.gH For the VeLLOs associated with L1521F PI, LI 148 PI 
(a VeLLO candidate), and L1014 PI, however, w e use the bet- 
ter co nstrained "internal luminosities" de rived by iBourke et alj 
(120061) , Kauffmann et al. (in prep.) and lYoung et alj d2004al) . 
Given the observed luminosities, the temperatures assumed for 
class I sources are (15 + 3) K, while for class sources they reach 
from 10 to 24 K and 21 K is derived for the only class II source 
(in L1021 PI; see below for YSOs and their classes). 

The uncertain density profile leads to an uncertain mass es- 
timate. We consider a range i n the expo nent p from 2 to 3/2, as 
expected for YSO envelopes (IShulll977h . in order to gauge the 
uncertainties in the mass estimate. This shows that the aperture 
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averaged temperatures in power-law envelopes, and correspond- 
ingly the derived masses, are uncertain by 30% and less. 

For the dust near class II YSOs we adopt an opacity of 

0. 02 cm 2 g" 1 because of an expected enhanced co agulation 
dMotte & Andrei 1200 ll lOssenkopf & Henningl |1994|) . For the 
other YSO and starless cores we adopt the standard c2d opacity, 
which is 0.0102 cm 2 g -1 at MAMBO's observing wavelength. 

Circumstellar disks might also contribute to the millimetre 
continuum emission. Their contribution to apertures with radii 
> 4200 AU is estimated to be < 10% for YSOs of the classes 

and I dMotte & A ndre 2001). The emission from the disks 
does therefore not significantly bias our envelope mass estimates 
for the youngest YSOs. Class II YSOs, in contra st, are not ex- 
pected to have envelopes (e.g.. lAndre et al.l l2000 and references 
therein); their emission is dominated by compact (< 1 000 AU) 
disks with a temperature above the one estimated from Eq. ([8]). 
Our procedure will thus overestimate the mass for YSOs with 
significant disk emission. 

3.3. Overview of Source Properties 

We u se the morphology keywords adopted by lEnoch etal] 
(2006) to describe our maps. These are summarised in Table Q] 
and can be employed to compare the dense core morphologies 
revealed by our maps with those found in other c2d bolometer 
surveys. 

All dense cores that were observed are detected in our maps, 

1. e., the maximum SNR in the map is at least 4. For one source 
only 'weak' emission is detected, meaning a peak SNR equal to 
or below 5; all other sources are brighter. In total 21 cores (55% 
of all cores) have at least two peaks separated by less than 3', 
and we consider them to be 'multiple'. All cores are 'extended', 
as they contain at least one subcore with an equivalent radius ex- 
ceeding 30". We take all subcores with a major-to-minor axis ra- 
tio exceeding 1.2 to be 'elongated', and the others to be 'round'. 
Then 6 cores (16%) only contain round subcores, 21 cores (55%) 
only contain elongated subcores, and 9 cores (24%) contain both 
round and elongated subcores. 

Figure [3] gives an overview over the sizes, masses, column 
densities, and ellipticities of well reconstructed subcores. These 
have typical values in the range (1 to 6) TO 4 AU = 0.05 to 0.30 pc 
for the effective radius, 0.5 to 20 M for the subcore masses, 
0.1 to 1.0 M Q for the mass within 4200 AU from the peak, and 
major-to-minor axes ratios at 50% peak intensity ranging from 

1 to about 4. Implications from these properties for the physical 
state and evolution of dense cores and YSOs are discussed in 
Sec.g] 



4. Supplemental YSO Data 

For the c2d MAMBO survey to be fully exploited, it needs to 
be complemented with information on the YSOs covered by the 
maps. Most of our analysis below is based on IRAS data since 
c2d finally observed only half of the YSO (& VeLLO) cores in 
our sample (i.e, 8 out of 17) and a homogeneous characterization 
is not possible based on c2d Spitzer data only. Also, we only use 
c2d photometry for Spitzer source characterization to guarantee 
homogeneity, implying that some images in the Spitzer archive 
cannot be used for the present study. All YSOs visible in Spitzer 
images are, however, either characterized by IRAS or Spitzer, 
with the singular exception of a Spitzer source in L1582A not 
covered by c2d data and not detected by IRAS. 



Table 7. IRAS sources in the MAMBO maps that are not de- 
tected as MAMBO emission peaks. For every core in the c2d 
MAMBO survey the table lists sources from the IRAS Point 
Source Catalogue and the IRAS Faint Source Catalogue (pre- 
ceded by an 'F'), and the subcores in which the IRAS sources 
are located, if applicable. Objects from the IRAS Point Source 
and Faint Source catalogues present in both catalogues are con- 
nected. The last column gives the spectral index between 12 and 
25 jt/m wavelength. 



Core 


IRAS Source 


Asso. 


25 fjan 


L1521-2 


F04262+2654 


CI 


0.27 ± 0.43 


L1524-4 


04274+2420 




-0.85 + 0.27 




L F04274+2420 




-0.83 ± 0.19 


B18-1 


04292+2427 




< -1.5 




L F0429 1+2427 




< -1.9 


TMC-2 


04294+2413 




-0.9 ±0.15 




L F04294+2413 




-0.89 ±0.11 


B18-4 


04326+2405 


CI 


< -1.68 




L F04326+2405 




< -0.75 


TMC-1C 


04380+2553 




-0.51 ±0.15 




L F04380+2553 




-0.76 ±0.12 


1ML-1 


r04383+2549 




< —1.12 




04392+2529 






L1582A 


F05290+1229 




-0.51 ±0.17 


L 1622 A 


05517+0151 




-2.14 ±0.16 




LF05517+0151 




-2.25 ±0.15 




05519+0148 


[C2] 






05519+0157 






L438 


18116-0707 




< -1.33 


L492 


18130-0341 








18132-0350 




-1.82 ± 0.15 


CB188 


19180+1127 




< -1.31 


L675 


19217+1103 




-1.81 ±0.25 


L1148 


20395+6714 








20410+6710" 


C2 




L1228 


F20598+7728 




-0.26 ± 0.34 


Ll 103-2 


21399+5632 




> -1.47 


L1251A 


F22282+7454 




< -2.54 



Notes: a) Probably an artifact since no corresponding Spitzer source 
exists. 



4.1. Associated IRAS, 2MASS, and Spitzer Sources 

Table [7] l ists sources from the IRAS Point Source Catalogue 
(PSC peichman et al][l9 88) and IRAS Faint Source Catalogue 
(FSC: lMoshir & et al.lll990» not associated with intensity peaks 
in the MAMBO maps, and Table[8]lists those that are associated. 
The sources are considered to be associated with subcores and 
peaks detected by MAMBO if the separation between an IRAS 
source and a MAMBO subcore or peak is less than the uncer- 
tainty of the separation at the 2cr level (i.e., the IRAS uncer- 
tainty ellipse plus the MAMBO pointing error). These sources 
are most likely physically associated with the dense cores and 
are assumed to be young stars. 

We include IRAS sources of any quality. This in principle 
could lead to considerable bias of our sample, since also extra- 
galactic sources and late type stars would be picked up by these 
criteria. However, only a single IRAS source projected onto well 
detected emission is not detected as a MAMBO intensity peak 
(i.e., in L1521-2 CI; see Table O; all other IRAS sources that 
lie within subcore boundaries but are not detected as MAMBO 
peaks are either detected as extensions in the MAMBO intensity 



Table 8. Properties of MAMBO-detected IRAS sources. For every core in the c2d MAMBO survey the table lists dust emission peaks and, if applicable, the related subcore 
associated with an IRAS source (Sec. |4.1| discusses the identification of associated sources; FSC sources are preceded by an 'F'). Associated 2MASS sources are listed too. The 
table further lists the spectral index between 12 and 25 yum wavelength, the bolometric temperature and luminosity, the submillimetre-to-bolometric luminosity ratio, the mass 
within a peak-centered aperture of 4 200 AU radius, and the infrared SED class. 



Core 


Asso. 


IRAS Source 


2MASS Source 


25 urn 
12 fim 


T IRAS 
' bol 


r IRAS 
Hoi 


r IRAS it IRAS 
submm ' bol 


m yso 

' H 4200AU 


Class 


Jens 












K 


Le 








Kauffma 


B18-1 


[C2], [P5] 


04292+2422 
& F04292+2422 


04321540+2428597 


0.84 + 0.13 


< 602 


> 1.09 


> 0.014 


> 0.179 


I 


B18-4 


CI, PI 


04325+2402 
& F04325+2402 


04353539+2408194 


> 3.76 


< 83 


0.42 to 0.79 


0.053 to 0.099 


0.522 to 0.589 


I 


nn et 


TMC-1C 


P3 


04385+2550 
& F04385+2550 


04413882+2556267 


0.41 ±0.13 


> 596 


0.23 to 0.47 


0.007 to 0.015 


0.042 to 0.047 


I 


B- 


TMC-1 


C2, P2 


04381+2540 


04411267+2546354 


1.52 ±0.14 


131 


0.67 ± 0.02 


0.043 ± 0.017 


0.379 ± 0.010 


I 


AME 






&F0438 1+2540 
















IRAS05413 


C2, P2 


05412-0105 


05434630-0104439 


> 0.25 


> 67 


0.52 to 7.26 


0.008 to 0.115 


0.308 to 0.597 





O 






&F0541 1-0106 
















Mapping of S 




C3, P3 


05413-0104 


none 


> -0.70 


< 55 


~ 10.5 


« 0.017 


~ 0.95 





CB188 


CI, PI 


19179+1129 


19201494+1135400 


> 0.88 


< 248 


1.74 to 2.14 


0.015 to 0.018 


0.272 to 0.288 


I 


LI 100 


CI, PI 


20355+6343 
& F20355+6343 


none 


0.70 ±0.18 


< 103 


1.22 to 2.13 


0.023 to 0.041 


0.445 to 0.515 


I 


L1041-2 


C4, P4 


20361+5733 
&F20361+5733 


none 


> 3.12 


< 60 


3.91 to 5.07 


0.025 to 0.033 


0.900 to 0.966 





pitzer 


L1148 


CI, PI 


F20404+6712 


20405664+6723047 


> -0.07 


> 110 


0.06 to 0.40 


0.023 to 0.165 


0.115 to 0.146 


I 


o 

K> 


L1082C 


CI, PI 


20503+6006 


none 


> 1.05 


> 62 


0.32 to 1.73 


0.026 to 0.141 


0.402 to 0.576 





CL 

in 
3 






& F20503+6007 
















L1082A 


CI, PI 


20520+6003 


20531346+6014425 


0.39 ± 0.20 


< 198 


1.00 to 2.49 


0.018 to 0.044 


0.388 to 0.491 


I 


B, 




C3, P5 


20526+5958 


none 




< 119 


0.84 to 1.36 


0.038 to 0.061 


0.538 to 0.602 


I 


Q 


LI 228 


CI, PI 


20582+7724 
& F20582+7724 


20571294+7735437 


0.35 ±0.11 


307 


2.27 ± 0.04 


0.035 ± 0.008 


0.782 ± 0.007 


I 


ouds : 


Bern48 


CI, PI 


21004+7811 
&F21004+7811 


20591408+7823040 


-0.17 ±0.10 


628 


11.05 ±0.20 


0.008 ± 0.001 


0.483 ± 0.003 


I 


andC 


L1172A 


C2, P2 


21017+6742 


21022122+6754202 


> -0.13 


< 100 


0.37 to 0.59 


0.049 to 0.077 


0.370 to 0.403 


I 


ores 






&F21017+6742 
















L1177 


CI, P2 


21169+6804 
&F21 168+6804 


21173862+6817340 


1.26 ±0.27 


< 127 


2.50 to 3.32 


0.038 to 0.051 


1.030 to 1.114 


I 




L1021 


PI 


21197+5046 


21212751+5059475 


-0.82 ±0.13 


> 990 


4.09 to 6.91 


0.001 to 0.002 


0.021 to 0.024 


II 




L1251A 


C2, P3 


22290+7458 
& F22290+7458 


22300004+7513578 


> 0.68 


> 88 


0.12 to 0.75 


0.033 to 0.216 


0.284 to 0.380 


I 
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contours (in B 18-4 CI), projected onto noisy parts of a MAMBO 
map (in LI 622 A C2), or of dubious quality and discarded any- 
way (in LI 148 C2). 

The criterion of association is relaxed for IRAS05413 P2, 
LI 100 PI, and L1082C PI. For these sources inconsistencies 
between the IRAS PSF and FSC indicate positional errors of 
order 1', and we associate them with neighboring MAMBO 
peaks since these would be unusually bright and compact for 
starless peaks. For IRAS 20410+6710, which is projected onto 
LI 148 C2, our more sensitive Spitzer maps show no correspond- 
ing source in the MIPS bands. We therefore consider this IRAS 
source to be an artifact. 

For those IRAS sources associated with dust emission peaks 
it is possible to significantly improve the accuracy of their po- 
sition by adopting the position of the MAMBO peak. This then 
allows for a search of 2MASS counterparts of these sources. We 
do so by searching for 2MASS sources less than 10" away from 
the MAMBO peak. If several 2MASS sources are found we as- 
sume the one closest to the MAMBO peak to be the counterpart. 
The identified counterparts are listed in Table[8] If no counterpart 
is found, 2MASS upper limits are taken to be similar to those of 
nearby 2MASS sources with upper limits. 

This combined data yields a spectral coverage with data 
near 1 /mi wavelength (from 2MASS), a well sampled range 
from 12 to 100 fim (from IRAS), and information at 1 200 fim 
(from MAMBO) . Fo r consistency with p revious work (e.g., 
lAndre et al.ll 19991 and iMotte & Andrell200ll) . and to avoid prob- 
lems of distance bias, we use dust emission flux densities for the 
4200 AU aperture to study the spectral energy distributions. 

In Table [9] we also present a list of Spitzer point sources from 
the 3rd c2d data deliverjQ associated with MAMBO dust emis- 
sion peaks. To identify associated objects we search our maps 
for Spitzer sources that are detected as point sources at 24 fim 
wavelength, have a 24 fim flux > 2 mJy, and are offset by less 
than 20" from a dust emission peak. We use the 24 fim data for 
our search since in this band Spitzer reliably picks up sources 
with infrared emission apparently in excess of the photospheric 
emission. 

The choice of the flux and offset cuts is guided by previous 
knowledge about sources in our sample. The 24 fim flux cut for 
our search must be < 2 mJy, since manual inspection of our 
maps reveals YSOs of such flux in our sample (e.g., SSTc2d 
J223105.6+751337 in L1251A). The offset cut must be > 10", 
larger than offset s recently observe d for YSOs in our sample 
(L1014-IRS; see iHuard etall 120061) . Using the c2d processed 
source catalogues (documented as part of c2d's 4th data deliv- 
ery) for the area surveyed by the Spitzer Wide-Are a Infrared 
Extragalactic survey (SWIRE; Lonsdale et al. 2003), with our 
criteria we expect to find 1.07 chance alignments between ex- 
tragalactic background objects and one of the about 110 dust 
emission peaks in our MAMBO maps. 

Table [10] furthermore gives a list of prominent YSOs de- 
tected in Spitzer data not taken by c2d, defined as previously 
known IRAS YSOs. This excludes sources in L1622A, since 
these are not well covered by our MAMBO maps. One source 
in L1582A that has apparent outflow nebulosity and was not de- 
tected by IRAS is also included. We cannot, however, associate 
a MAMBO peak to this star, and we therefore have to exclude it 
from the further analysis. Manual inspection of the images did 
not reveal further probable YSO candidates. No c2d photom- 
etry exists for these images, and the Spitzer source properties 

7 See http://ssc.spitzer.caltech.edu/legacy/c2dhistory.html for a de- 
tailed documentation. 



Table 10. Prominent YSOs detected in Spitzer data not taken by 
c2d, and related IRAS sources, if existing. 



Core 


IRAS Counterpart 


B18-4 


04325+2402 


L1582A" 


none 


IRAS05413 


05412-0105 




05413-0104 


L1172A 


21017+6742 


L1177 


21169+6804 



Notes: a) MIPS source at 05 32 02.9, +12 31 05 (J2000.0) 

Table 11. Cores in the c2d MAMBO survey with YSOs uncov- 
ered by c2d that were believed to be starless at the beginning of 
the c2d survey. Comments on outflows refer to jet-like extended 
(» 1') features in IRAC images. 



MAMBO Core 


Comments 


L1521F 


CI 


VeLLO. Bourke et al. (2006) 


L1582A 


C2° 


no c2d Spitzer data, possible outflow 


L673-7 


CI" 


VeLLO, Dunham et al., in prep. 


LI 148 


CI 


VeLLO, Kauffmann et al. 2005 and in prep. 


L1014 


CI 


VeLLO. Young et al. (2004a) 


L1251A 


C3" 


new class 0, prominent outflow 




C4" 


new class 0, possible outflow 



Notes: a) reported here for the first time 



are not further discussed in the following, except when refining 
YSO positions (Sec. l5.4.61 >. Including the object in L1582A, the 
combined Spitzer data thus shows 7 YSOs in cores previously 
believed to be starless (Table ITTb. Four out of these (including 
the VeLLO candidate in L673-7) are reported here for the first 
time. 

The data on the MAMBO-associated Spitzer sources is com- 
plemented by 2MASS data where possible. Given the nominal 
positional uncertainties, 2MASS sources within 2" from Spitzer 
sources are assumed to be counterparts of these. For Spitzer 
sources without 2MASS counterparts upper limits to their flux 
densities in the 2MASS bands are derived from upper limits for 
nearby 2MASS sources not detected in all filters. The source 
emission is thus probed in the 1 to 70 fim wavelength range (by 
2MASS and Spitzer; no c2d Spitzer data is available at 160 fim) 
and at 1 200 fim (from MAMBO for the 4 200 AU aperture). 

Note that probably many more Spitzer sources are associated 
with the MAMBO cores observed in our study. However, here 
we are only interested in sources directly associated with fea- 
tures (i.e., peaks) seen in our MAMBO maps. The discussion of 
the other sources — which are expected to belong to the YSO 
classes II and III, since they lack significant millimetric dust 
emission — is deferred to a later paper. Also note that the emis- 
sion of some YSOs might be spatially confused with dense core 
emission. For example, the PSC source 04326+2405 in B18-4 
(marked by the only circle) might manifest as an extension in 
the MAMBO dust emission contours, which however cannot be 
uniquely separated from the dense core emission. 



Table 9. Properties of MAMBO detected c2d Spitzer sources. Like Table[HJ but now listing sources from Spitzer point source catalogues and giving the spectral index between 
2.2 and 24 /im wavelength. VeLLO candidates are marked in the class designation. 



Field 


Asso. 


c2d Spitzer Source 


2MASS Source 


24 pi 
2.2 




7"SST 
1 bol 

K 


rSST 
Hoi 

Le 




^SST IJ_SST 
submm ' bol 


m yso 

J ™4200AU 

M 


Class 


L1521F 


CI, PI 


SSTc2d J042839.0+265135 


none 


1.56 ± 


0.08 


< 30 


0.06 


to 0.08 


0.400 to 0.545 


0.729 to 0.749 


0, VeLLO 


B18-1 


[C2], [P5] 


SSTc2d J043215.4+242859 


04321540+2428597 


-0.25 


±0.06 


no IRAC data 


available 








TMC1-1C 


P3 


SSTc2d J044138. 8+255627 


04413882+2556267 


-0.04 


±0.05 


681 


0.37 


±0.03 


0.008 ± 0.023 


0.022 ± 0.005 


II 


TMC-1 


C2, P2 


SSTc2d J0441 12.7+254635 


04411267+2546354 


0.96 ± 


0.06 


no IRAC data 


available 








CB188 


CI, PI 


SSTc2d J192014.9+1 13540 


19201494+1135400 


0.25 ± 


0.05 


433 


1.12 


±0.09 


0.024 ± 0.022 


0.322 ± 0.007 


I 


L673-7 


CI, PI 


SSTc2d J192134.8+112123 


none 


0.97 ± 


0.08 


< 55 


0.03 


to 0.06 


0.284 to 0.573 


0.375 to 0.394 


0, VeLLO 


LI 148 


CI, PI 


SSTc2d J204056.7+672305 


20405664+6723047 


0.73 ± 


0.05 


< 140 


0.05 


to 0.10 


0.074 to 0.162 


0.139 to 0.148 


I, VeLLO 


LI 228 


CI, PI 


SSTc2d J205712.9+773544 


20571294+7735437 


0.14 ± 


0.05 


388 


1.82 


±0.13 


0.042 ± 0.027 


0.832 ± 0.007 


I 


Bern48 


CI, PI 


SSTc2d J205914.0+782304 


20591408+7823040 


-0.42 


± 0.10 


742 


8.35 


±0.51 


0.009 ± 0.013 


0.265 ± 0.002 


II 


L1014 


CI, PI 


SSTc2d J212407.5+495909 


none 


0.33 ± 


0.07 


< 163 


0.04 


to 0.15 


0.118 to 0.411 


0.318 to 0.357 


I, VeLLO 


L1251A 


C2, P3 


SSTc2d J222959.5+75 1404 


22300004+7513578 


0.51 ± 


0.05 


< 197 


0.11 


to 0.35 


0.060 to 0.187 


0.327 to 0.382 


I 




C3,P4 


SSTc2d J223031. 8+751409 


none 


0.57 ± 


0.08 


< 44 


0.17 


to 0.30 


0.170 to 0.301 


0.831 to 0.900 







C4, P5 


SSTc2d J223105. 6+751337 


none 


-0.23 


±0.07 


< 63 


0.08 


to 0.17 


0.201 to 0.421 


0.630 to 0.683 
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4.2. YSO Properties 

Below we characterise the sources associated with dust emission 
peaks. All our target cores were covered by the extensive but 
insensitive IRAS and 2MASS surveys. These data thus allows 
for a homogeneous quantification of all sources in our survey 
that are brighter than a few times 0. 1 L . Fainter sources could 
only be detected by Spitzer, from which data exists only for a 
fraction of our dense cores. 

Given the differences in the spectral bands probed by these 
instruments, it is difficult to compare a source only detected by 
2MASS and IRAS to one only detected by 2MASS and Spitzer. 
In order to explore and suppress related biases in estimates of 
source properties, we therefore analyse the Spitzer and IRAS 
data separately. To better distinguish results from the different 
approaches, in the following properties are labeled by the data 
source used to derive them (superscripts TRAS' and 'SST', the 
latter for the Spitzer Space Telescope). 

4.2.1. Estimates from IRAS and 2MASS Data 



Following iLadal d 1987b . the spectral properties of young stars 
can be characterised by the spectral index between 12 and 25 /mi 
wavelength, 



a 



25 /jm 
12 /jm 



log(12 pm ■ F 2 5 nm/[25 (im ■ F n ^ m ]) 
log(25 /mi/ 12 pm) 



(9) 



Spectral indices are listed in Tables|7]and[8]for all IRAS sources 
in the MAMBO maps, if detected in these bands. They roughly 
probe whether the spectral energy distribution (SED) is dom- 
inated by photospheric or envelope emission and are some- 
times use d to classif y observed YSO SEDs within evolutionary 
schemes dLadalll987l) . 

For the IRAS sources associated with MAMBO peaks for 
which 2MASS data is available, we calcul ate the bolometric 
temperature defined by My ers & Laddl ( 1 19931) . 



Tbol - 



m h(v) 

4<r(5) * B 



(10) 



where (v) = vF v dv / J^°° F v dv is the flux- weighted mean fre- 
quency, ( is the Riemann zeta function, and h and k B are Planck's 
and Boltzmann's constant. We integrate across the SED by inter- 
polation between observed 2MASS and IRAS bands. For this 
we use piecewise power laws matching the flux densities, as 
indicated in Fig. [2] At 500 /im wavelength these connect to a 
power-law representing a modified blackbody of opacity oc A -2 
and 15 K temperature that is tuned to fit the MAMBO data; 
piecewise powerlaws are adopted to hold between A = 500 /im 
and 1200 pm (which is MAMBO's central wavelength), and 
A = 1200 //m and oo (adopting a flux density distribution oc A~ A ). 
For opacity-modified blackbodies of 10 to 20 K temperature this 
approximates the actual flux density at A — 500 /im within a 
factor of 1.5. The related uncertainty on integrated properties is 
much lower though, since the submillimetre part of the SED usu- 
ally only contributes a few percent to the total luminosity (see, 
e.g., Table[8]l. The integrati on extends from the 2MASS bands at 
ss 1 /im to infinity. Unlike iMvers & Laddl we treat all flux den- 
sity upper limits like actual detections. Therefore, the derived 
bolometric temperatures must be interpreted with some caution, 
if a source is not detected in some of the bands. In these cases 
lower and upper limits are assigned to the bolometric temper- 
ature, depending on whether adopting flux densities below the 



o 
o 



10" 



10" 



10" 



10" 



10" J r. 



o 
o 



10" 



10" 



10 



10" 



L1251A 

SSTc2d J222959. 5+751 404 



I ' ' — I 

l M = 0.11 to 0.35 L» 
T bol = <197 K 



Spitzer IRAC & MIPS , 




interpolated from 
MAMBO 



r 



1 1 matched blackbodies 
' ( (truncated for clarity) 



mean wavelength 



MAMBO 



B18-4 

IRAS 04325+2402 



2 MASS 



IRAS 



L bol = 0.42 to 0.79 L sun 
T„„, = 83 K 

. interpolated from 
MAMBO 




MAMBO 



mean wavelenght 



10 10 
wavelength [/im] 



10 J 



Fig. 2. Example SEDs for YSOs in the c2d MAMBO 
sample. The top panel shows the infrared star SSTc2d 
J222959.5+751404 associated with the P3 peak in the C2 sub- 
core of LI 251 A, while the bottom panel shows the data for IRAS 
04325+2402 associated with the PI peak in B18-4. The SEDs 
are composed from 2MASS, Spitzer or IRAS, and MAMBO 
data, as indicated by labels. Also indicated are the interpolated 
datapoints at 500 /im that is calculated from the MAMBO flux 
density assuming a dust opacity oc A~ 2 and a dust temperature of 
15 K. Straight lines indicate the flux density interpolation used to 
calculate the bolometric luminosity and temperature indicated in 
the upper right corners (as well as the mean wavelengths calcu- 
lated as c/(v) from Eq. IfTUl and indicated at the frame bottom). 
Only limits can be derived on Tboi and Lbol because of some 
upper limits in the 2MASS data; a conservative lower limit to 
the luminosity is set by the sum of opacity-modified blackbod- 
ies matched to the SED at short and long wavelength (dashed 
lines, drawn for opacities scaling oc y° to 2 ; see text for details). 
The "sign" of the bolometric temperature limits is set by the 
mean wavelength lying at longer wavelength than all photomet- 
ric bands with upper limits. 



upper limits would increase or decrease the calculated mean fre- 
quency, respectively. 

Based on the derived bolometric temperature the YSOs are sep- 
arated into four infrared classes. Following IChen et aT I (fl997l 
1995) all stars with 7/ bol < 70 K are considered to be class 
sources. Those with 70 < 7boi/K < 650 belong to class I, for 
650 < 7/boi/K < 2880 they belong to class II, and class III 
sources have 7boi > 2880 K. In case of upper or lower limits 
to the bolometric temperature, the class corresponding to the de- 
rived limit is given. This treatment of uncertainties in class as- 
signments is somewhat unsatisfying, however it does not influ- 
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ence our later analysis that is based on the numerical value of 

Sources in a given class (i.e., with similar bolometric tem- 
perature) are beli eved to be in a simil ar ev olutionary stage 
dChen et alJ 1 19951 and IChen et al.l Il997t see lYoung & Evansl 
2005 for illustrative examples). In their evolution from a deeply 
embedded object to a star surrounded by a remnant dust disk, 
YSOs are though to pass through phases to III in increasing 
order. Thus class objects would represent the youngest YSOs 
and class III objects would be in a stage just prior to the main 
sequence. Almost all of the MAMBO-detected stars covered by 
the present survey are in classes and I; only one class II and no 
class III object are among them. 

From the interpolated SED we derive an estimate for the bolo- 
metric luminosity, 



F v dv, 

o 



(11) 



for IRAS sources associated with MAMBO dust emission peaks. 
If the source is detected in a few bands only, an upper limit to the 
bolometric luminosity is calculated by integrating across the flux 
density upper limits, while a lower limit is given by the sum of 
the maximum 2MASS "monochromatic luminosity", vF Y , plus 
the maximum detected vF Y in the bands of longer wavelength 
(see Fig. |2] for an illustration). The latter builds on the fact that 
for blackbodies modified by an opacity scaling with frequency as 
v° to 2 , the bolometric luminosity is the peak value of vF v times a 
factor 1.05 to 1.36. In a similar fashion we derive the luminosity 
in the submillimetre wavelength range, 



•c/350 fim 



F v dv, 



(12) 



where c is the speed of light. For both luminosities we also quote 
uncertainties due to statistical noise, while we do not consider 
systematic contributions. The aforementioned uncertainties in 
the submillimetre SED near 500 yum wavelength do, however, 
render the submillimetre luminosity uncertain by a factor of or- 
der 2; no precision measurements of it from our data is possible. 

4.2.2. Estimates from Spitzer and 2MASS Data 

The YSO properties derived from Spitzer and 2MASS data 
(Table [9]) are calculated using methods similar to those adopted 
for the combined IRAS and 2MASS data. The spectral sam- 
pling at wavelengths > 25 jum is, however, much worse than 
for the combined 2MASS and IRAS dataset, resulting in more 
uncertain estimates for cold sources (with SED peaks at wave- 
lengths > 25 yum, i.e. with 7b i S 160 K). The spectral index 
for the 2MASS and Spitzer bands, a^i^ is derived by fitting 
the monochromatic luminosities in the 2.2 to 24 fim wavelength 
range by a power law. It is, thus, not directly comparable to the 
spectral index derived from IRAS data. 

Four stars in our sample might qualify as VeLLOs, given 
that they appear to be embedded in dense cores and have in- 
ternal luminosities likely below 0.1 L . This group includes the 
sources associated with L1521F PI, LI 148 PI, and L1014 PI 
which have been subject to detailed studies (|Bourke et"aT1 l2006; 
iKauffmann et alj|2005l and in prep.; I Young et al. 2004al) . These 
studies confirm the spectral properties found here. In particular, 
based on detailed radiative transfer SED modeling, they support 
the prevalence of internal luminosities < 0.1 L Q . This is an im- 
portant factor, since the SED analysis presented here has bad 



spectral sampling at long wavelength, leading to uncertain prop- 
erties for cool sources. It is also not suited to single out the in- 
ternal luminosity due to an embedded VeLLO. Note that SED 
modeling of the source in L673-7 has not been presented so far, 
leading to a more uncertain status of this source. 

Given the evidence for YSO outflows from scattered light 
nebulosity and outflows seen in CO, the stars in L1521F PI, 
L1014 PI, LI 148 PI, and L673-7 PI are proven to be associ- 
ated with the respective cores and are not unrelated background 
objects (see the above refere nces and additiona lly Dunham et 
al., in prep., for L673-7, and lBourke etal]|2005l for L1014). A 
recent search for YSOs with internal luminosities < 0.1 L Q by 
Dunham et al. (in prep.) lists three of these sources in their 
"group 1" of confirmed low-luminosity YSOs; L1448-IRS is 
listed in their "group 3" of good, but unconfirmed, candidates 
for low-luminosity YSOs. The star in L1251A P5 however, for 
which our above analysis suggests a low luminosity, is not likely 
to be a VeLLO since its luminosity limits appear to be too high. 

4.2.3. Properties chosen for the Analysis 

The above analysis yields two different sets of properties, one 
based on Spitzer data, one based on IRAS. For the bright sources 
of a few 0. 1 Lq, the bolometric temperatures and luminosities 
estimated from IRAS data are used for the further YSO analy- 
sis. This is the preferred option because IRAS data is available 
for all our cores, and because Spitzer fails to detect a signifi- 
cant fraction of the YSO emission at long wavelengths (see Sec. 
15.4.31 leading to biased estimates. We abstain from constructing 
SEDs combining flux densities from IRAS and Spitzer in order 
to better control mission-related systematic trends in the data. 
The IRAS data might, on the other hand, suffer from contami- 
nation of unresolved nearby sources. Such contamination is not 
apparent at a significant level in any of the cases where we have 
data from both IRAS and Spitzer though. 

Fainter sources are usually only detected by Spitzer. Thus, 
Spitzer data is used for the study of the sources in L673-7 
PI, and L1251A P4 and P5. For L1521F PI, LI 148 PI, and 
L1014 PI results fr o m the more involved Spitzer data analy- 
sis by | Bourke et al] d2006l) . Kauffmann et al. (in prep.), and 
lYoung et al.l d200 4a) are used. 



5. Analysis 

In this section we exploit the survey data in order to better un- 
derstand the state and evolution of starless and YSO dense cores. 
We discuss the general core properties before we turn to the dis- 
cussion of specific issues. Throughout the following discussion 
we exclude starless peaks and subcores with uncertain proper- 
ties, unless noted otherwise. 



5. 1 . General Core Properties 

Figure [3] shows the frequency distributions of some properties 
of well reconstructed dense cores. Some conclusions about the 
physical state of individual cores, and of dense cores in general, 
can be derived from these distributions. 

5.1.1. Radius 

The sizes of YSO cores are on average larger than those of star- 
less and candidate VeLLO cores. This could to some extent be 
an observational bias. Only one YSO subcore in our sample is at 
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Fig. 3. Overview of the dense core properties, excluding subcores and peaks with uncertain properties due to artifacts. The dashed 
line represents the properties of starless peaks and subcores, while the grey and black areas hold for YSO and candidate VeLLO 
cores, respectively. Bullets of the same colours, and the related numbers, give the respective median values of the distributions. 
Vertical lines in panel a) give the beam radius at distances of 100 pc and 450 pc, respectively. The vertical lines in panel c) 
show the range of masses within a 4200 AU aperture for singular isothermal spheres of 10 K temperature and truncation radii of 
4200 AU — > oo . The vertical line in panel d) marks the boundary between round and elongated sources. Typical uncertainties of 
the shown properties are > 10% for M tot , a few 0.01 M Q for M4200AU, and 10% to 20% for the axis ratio. 



a distance below 200 pc, while two dozen starless subcores are 
within this distance. Because confusion increases with distance, 
the subcores identified in the more distant YSO cores are thus 
likely biased towards larger sizes. Large cores might, however, 
indeed have better chance to form stars. 

At small radii the size distribution of starless cores is lim- 
ited by the beam size; their number steeply drops to zero at the 
resolution limit. We might miss a population of small subcores. 

5.1.2. Total Mass 

The mass of the subcore L1622A C3 exceeds the mass of all 
other subcores by a factor 1.8 and more. Given that L1622A is 
one out of only 3 dense cores located in the Orion star forming 
complex, in which special environmental conditions prevail, it 
might be that L1622A C3 has, e.g., a mass much larger than the 
other cores in our sample. However, the extreme mass contrast 
to all other cores does cast some doubt on this. 

It might be that L1622A is warmer than most other cores. 
By adopting a dust temperature of 10 K one would then over- 
estimate the true mass. Furthermore, substructure easily disen- 



tangled in more nearby cores may be confused in L1622A, the 
core with the largest distance in our sample (450 pc). Then the 
subcore area, and therefore the mass, would be biased towards 
higher values. This is supported by the fact that L1622A C3 has 
the largest effective radius of all subcores. Also, L1582A C2 and 
L1622A C4, which are second to L1622A C3 in mass and are 
also located in Orion, are second to L1622A C3 in the effective 
radius. This may hint that mass and size estimates of cores are 
indeed biased for larger distances. 

Thus, data on distant cores must be interpreted with some 
caution. They might well have extreme properties, but given the 
lack of sufficient spatial resolution, we cannot be sure about this. 
Properties of nearby and distant cores are not necessarily directly 
comparable. 

5.1.3. Mass within 4200 AU 

For L1082A P3 the mass M4200AU exceeds those of all other 
peaks by a factor 1.7 and more. This large aperture mass con- 
trast may hint on a biased mass estimate. The dust in L1082A P3 
might be significantly warmer than 10 K, if this core is heated by 
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an undetected embedded star. Following My ers et ail d!987h . the 
failure of IRAS to detect a point source in this peak implies (for 
a distance of 400 pc) an upper limit to the bolometric luminos- 
ity of ss 0.8 L to any embedded source. In this case the average 
dust temperature could be 15 K instead of the assumed 10 K, 
and the mass could be overestimated by us by up to a factor of 
2. Then M4200AU would not be unusual for a YSO core. Though 
this evidence is not conclusive, we, thus, suspect the presence of 
a heating YSO source in L1082A P3 that is too faint to be de- 
tected by IRAS. Unfortunately, no Spitzer data is available for 
this region to confirm this. 

The peak B18-4 P2 is separated from the other starless sub- 
cores. As Spitzer images show, no point sources are associated 
with this peak, and so internal heating cannot explain the un- 
usual value of M 42 ooau- A YSO of * 0.5 L Q resides in B18-4 PI, 
which is separated from B18-4 P2 by 10000 AU = 0.05 pc. At 
this distance a YSO of this low luminosity is unlikely to provide 
significant external heating to the peak P2. Thus the high value 
of M4200AU in B18-4 P2 is likely to be real. 

The M4200AU distributions for all types of subcores peak at a 
mass of * 0.35 M . The mean (0.42 M , 0.54 M , and 0.41 M 
for starless, YSO, and VeLLO cores, respectively) and median 
(0.35 M , 0.48 M , and 0.39 M , respectively) aperture mass is 
of the same order. Since this is a typical value for all kinds of 
subcores, the underlying physical process shaping the distribu- 
tions must be fundamental. Interestingly, this mass is similar to 
the mass within 4 200 AU for a (truncated) singular hydrostatic 
equilibrium isothermal sphere (SIS) of 10 K gas temperature, 



M4200AU =O.34M 0£ (7yiOK). 



(13) 



Here T g is the gas temperature and the factor e depends on the 
density distribution outside a radius of 4 200 AU; it is 1 if the 
density drops to zero outside 4 200 AU and becomes n/2 « 1 .57 
if it continues out to infinity. This mass constitutes a critical 
value for an SIS: no hydrostatic equilibrium solutions exist for 
higher values of A^ooau- It thus might be that this critical con- 
dition manifests in characteristic values of dense core properties. 

If this is true it has two interesting implications. First, dense 
cores would preferentially exist in a close-to-critical physical 
state. We note that interestingly such near-cr i tical c ore states 
have previo usly been inferred b y lAlves et alj ([2001 ) and par- 
ticularly by iKandori et al] (120051) . Second, the total pressure in 
the core, P, appears to be comparable to the thermal pressure. 
Contributions to the total pressure from the effect of turbulent 
gas motions or magnetic fields can not much exceed the ther- 
mal pressure. Otherwise the critical mass in the 4 200 AU aper- 
ture would be much larger since one would need to replace the 
gas temperature with some higher effective one, T g — > = 
2.33 Pmn/(gks) > T g , where ;«h is the hydrogen mass and the 
factor 2.33 holds for a gas mixture at cosmic abundance with 
most hydrogen in molecular form. 

5.1.4. Elongation 

Only 14% and 44% out of all starless and YSO subcores, respec- 
tively, have major- to-minor axis ratios < 1.2 (where the uncer- 
tainty is 10% to 20%) and can be considered round; most sub- 
cores are not round. Unless these subcores are shaped by mag- 
netic fields, which provide a non-isotropic supporting pressure, 
they can hardly on the whole be in a state of hydrostatic equilib- 
rium. 

The axis ratio distribution is continuous up to ratios of 2.0 
and 2.5 for YSO and starless subcores, respectively. Several 



subcores have much larger axis ratios. Interestingly, two of the 
most distant subcores have the largest elongations (L1082A C4 
and L1622A C3 at distances > 400 pc). This suggests that 
some cores appear elongated due to unresolved cores within the 
beam, an effect that on average increases with distance (also see 
lYoung et al] d2006b) for this effect in Bolocam maps). 

For starless cores, the median axis ratio that we derive, i.e. 



ersj 



1999 
1999 



1.9, is very similar to those derived from NH3 ( Jiiina et 
median of 1.5 to 2.2) and extinction maps dLee & Mver 
mean of 2.4 + 0. 1) of larger core samples. For YSO cores, how- 
ever, the median aspect ratio of 1.3 for our sample is at the low 
end of the aspect ratios derived by the above studies (medians of 
1 .4 to 2.0, respectively a mean of 2.2 + 0.2, for their YSO cores). 
This difference is probably an effect of the internal heating on 
the dust emission intensity distribution (see below). 

The axis-ratio distribution has in the past been used to 
study the intrinsic three-dimen s ional geometry of dense core 
shapes. Following My ers et ail dl99ll) . the median axis ratio 
of 1.9 for starless subcores (we exclude YSO cores, since in- 
ternal heating influences the observed shapes; see next para- 
graph) is the expected projected axis ratio of prolate ("cigar- 
like") spheroids with an intrinsic major-to-minor axis ratio of 
ss 2.5, and of oblate ("disk-like") ones with an axis ratio of 
> 4 (for an observed ratio > 1.7), respectively We cannot 
resolve this ambig uity. Some existing theoretical analysis fa- 
vors prolate cores dMvers et al.ll 1 99 lb [Rvden 1996), while other 
work prefers oblate o nes dJqnesetalhoOlM J ones & Basull2002l: 
iGoodwin et alj|2002tlTassisll2007l) . 

The median axis ratio of YSO subcores is smaller than that 
of starless and VeLLO cores. This could be because the YSO 
heating leads to an intensity peak of small physical size and low 
aspect ratio (since the heating, to first order, decreases radially 
away from the star; Eq. [5]), which is blurred into a structure of 
even lower aspect ratio because of beam smearing. On the other 
hand stars are expected to form in close-to-spherical densit y en- 
hance ments, which in projection show a low elongation (Shu 
[19771) . Then, YSO cores would be less elongated than the star- 
less cores with their complex shapes at lower intensity contours. 
Both effects would also explain the aforementioned decrease of 
the aspect ratio from the large-scale core structure probed by, 
e.g., NH3 and extinction maps to the core centers probed by dust 
emission. A more detailed analysis is needed to separate these 
effects. 

The YSO subcores CB188 CI, IRAS05413 C2, LI 100 CI, 
LI 251 A C2, and LI 251 A C4 are well separated from the other 
YSO cores in the elongation distribution. For IRAS05413 C2, 
however, the axis ratio is very uncertain (1.9+0.5). The observed 
elongation is likely to be an observational artifact. In L1251A, 
where the major axis exceeds the beam size by a factor of sev- 
eral, this elongation likely reflects the morphology of the dense 
core from which the embedded star formed. For CB188 and 
LI 100, however, where the major axis exceeds the beam size 
by a factor 3.2 and less, this elongation should reflect the mor- 
phology of the immediate density peak from which the young 
star accretes. This region could be shaped by the interaction 
with outflows from the c entral star. An outflow has indeed been 
detected toward CB188 dYun & Clemens! 1 1994 and for LI 100 
there is some eviden ce for an outflow from broad line wings 
dDe Vries et al]|2002l) . For CB188 the position of the outflow 



8 This analysis assumes a random distribution of the core axes in 
space. This might not be justified for our sample, since most of our 
cores are situated in a small number of star-forming regions, in which 
preferred orientations might prevail. 
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axis derived bv lYun & Clemens] ([1994) and the major axis of the 
dust emission intensity distribution are neither aligned nor per- 
pendicular (position angles of m 15° and 16° ± 5°, respectively). 
If the elongated structure towards CB188 seen in the MAMBO 
maps is indeed related to the outflows, then the dust emission 
feature is likely to trace the wall of an outflow cavity. If it would 
trace the jet or a YSO disk, the axes are expected to be parallel or 
perpendicular. While these observations do not prove the inter- 
action between jets and the dense core, they motivate dedicated 
studies on this issue. A more detailed analysis of the dust emis- 
sion maps might yield more candidates for jet-core interactions. 

5.2. Correlations between Core Properties 

Correlations betwe en core par ameters, like e.g. the famous size- 
linewidth relation dLarsonl 19811) . can provide important hints on 
the nature of dense cores. Such correlations are indeed seen in 
our data (Fig.|4|i. Unfortunately, many of these cannot be cleanly 
separated from observational biases, and great care is needed if 
one wishes to exploit our data in this fashion. 

To give examples, our source identification scheme directly 
leads to a radius-dependent lower limit to the mass, since the 
intensity will be at least twice the noise level within any sub- 
core, and thus M > 2fiji 2 mjj7rN-RMsQiz)'ig (where fiu^ and mu 
are the mean molecular weight per hydrogen molecule and the 
mass of the hydrogen atom, respectively, and A^rms(H2) is the 
column density corresponding to the intensity RMS; also see 
Appendix lAl. Similarly, we cannot detect massive cores of low 
peak column density, A^ max (H2), in maps of limited radius, r map , 
since M < ^H 2 m H^N m!ix (H2)rl 1!ip . Also, we might miss low-mass 
(and therefore small) cores of high column density because of 
beam smearing. Similar considerations apply regarding M4200AU) 
which can be interpreted as an aperture-averaged column density 
(that, however, does not suffer from beam smearing). 

Given the above uncertainties, we are careful in drawing con- 
clusions from apparent correlations. In other words, parts of the 
parameter space e.g. explored in Fig.|4]might not be populated 
because cores in these regions would not be detected given our 
sensitivities. 



5.3. Dense Cores and their Evolution 

Our data not only reveals cores in various states of evolution but 
in combination with the complementary YSO data we can also 
make statements on the progress of dense cores through evo- 
lutionary stages. The latter suggests some necessary conditions 
for active star formation to be possible, while it also indicates 
diverse evolutionary tracks for individual cores. 



5.3.1. Densities from Masses and Radii 

A rough estimate of the core densities can be derived from the 
peak aperture masses when assuming specific density profiles. 
The estimates will certainly often be off from the actual value by 
a large factor, but this analysis is still illustrative to understand 
trends in the data. 

We here compare the observed aperture masses and radii at 
70% peak intensity to those expected for dense cores with den- 
sity profiles n(H2) = n c (H2)/(l + [r/ri/2] 2 ), which provides a 
good approximatio n to the density structure of starless cores 
dTafalla et al.ll2002 |). In this n c (H2) is the central particle den- 
sity, while r is the radius from the center, and ri/2 is the radius 
at which n(H2) = n c (H2)/2. Particle densities can be related to 
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Fig. 5. Relation between the intrinsic radius at 70% peak in- 
tensity and the mass within 4 200 AU radius from the peak for 
VeLLO candidate and starless cores. The plotted effective radii 
deviate from the ones listed in Table |4]in that the beam smear- 
ing is removed (7'.'2 at 70% peak). See Fig. [4] for an explana- 
tion of the symbols. The radius bias due to central heating of 
VeLLO cores (Sec. |5.5.2t is indicated by indicated by arrows. 
The dotted lines indicate curves of constant central density for 
some simple density profile of varying radius (see the text for de- 
tails). The labels give the related central H2 column density. The 
dashed line indicates the upper radius limit for ev olved dense 
cores, < 4 800 AU, suggested bv lCrapsi etai1(l2005l) . The candi- 
date VeLLOs observed all fall below this limit and are, for given 
aperture mass, smaller than most starless cores, though this in 
part might come from internal heating. Also, most VeLLO cores 
appear to have central densities exceeding those of most starless 
cores. VeLLO cores are thus possibly physically more evolved 
than starless cores. This appears to distinguish them from star- 
less cores. 



mass densities using the c2d standard H2-to-gas mass conver- 
sion factors documented in Appendix|A] Figure [5] illustrates this 
analysis. 

This analysis indicates typical central densities of order 

3 • 10 4 cm 3 , while for some extreme cores we de rive densities 
of up t o 10 6 cm~ 3 . S i milar figure s were derived bylTafalla et ail 
(12002b . ICrapsi et al l (120041) . and iKirk et ail d2005l) . Note, how- 
ever, that the obtained densities are crude estimates that suffer 
from significant systematic uncertainties. Anyway, to give ex- 
amples, for L1521F CI the density of w 4 ■ 10 5 cnT 3 derived 
from our simplified analysis is in good agreement with the den- 
sit y of » 5 • 10^ cm" 3 (for our choice of dust properties) derived 
bv iCrapsietaHd2004l) . 

Note that VeLLO candidate cores appear to have unusually 
large densities, when compared to starless cores. Section [5.5.21 
discusses this in detail. 

Our sample is highly biased and inhomogeneous, but it is still 
illustrative to derive relative dense core lifetimes from the rela- 
tive numbers of dense cores of different physical state. Such es- 
timates assume a constant core formation rate across the sample 
and common evolutionary paths for all cores, neither of which 
are likely correct for our cores. Anyway, we find about 35 star- 
less subcores in the 10 4 < n c (H2)/crrT 3 < 10 5 density range, and 

4 subcores in the 10 5 < n c (H2)/cmr 3 < 10 6 range. Under the 
above assumptions this implies that dense cores reside at densi- 
ties of 10 4 to 10 5 cm 3 for a time exceeding the one for densities 
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Fig. 4. Necessary conditions for active star formation from relations between dense core properties. Filled and empty circles give the 
properties of subcores with well determined and uncertain properties due to artifacts, respectively. Stars and asterisks are used for 
subcores hosting YSOs and candidate VeLLOs, respectively. The diagonal dashed line in panel a) indicates a selection limit; for a 
given total mass all subcores larger than this limit would be too faint to be identified in our maps. The diagonal dashed lines in panel 
c) mark where M = 0.1 M4200AU and M = 0.01 M4200AU: on ly a fraction of the dense core mass is available to form stars from a 
column density peak. The curved dashed lines show the critical stability limits for hydrostatic equilibria with pure thermal pressure. 
For given total mass, or mass within the half intensity contour, all subcores with pure thermal pressure and with radii smaller than the 
critical one — respectively column density or aperture mass exceeding the critical one — are unstable. These stability limits might 
manifest in the data as necessary conditions for active star formation; most dense cores presently forming stars exceed these limits. 
However, these limits do not give sufficient conditions for ongoing star formation. Many starless cores have properties exceeding 
the critical values. In the core evolution framework discussed in the text this could be understood if these subcores are in hydrostatic 
equilibrium but supported by significant non-thermal pressure, or are collapsing. 



of 10 5 to 10 6 cm 3 by a factor m 9. For reference, the free fall 
timescalfl 

T ff = 9.8 • 10 4 yr [n(H 2 )/(10 5 cnT 3 )r 1/2 , (14) 

only varies by a factor lO 1 ^ 2 « 3 between the density bins. Dense 
core evolution in free fall would thus not explain the relative 
number of cores in different density bins. 

5.3.2. Evolutionary States 

Crap si etaTJ ([2005) suggested some criteria to assess a starless 
core's evolutionary state. These also include criteria based on 

9 We give the collapse time for a homogeneous sphere of density 
n(H2), i.e. = (3n/[32Gg]) l/2 , where G is the constant of gravity 
and q is the volume-averaged density. Our calculation assumes a mean 
particle weight of 2.8 per hydrogen molecule. 



observations of molecular lines, which we cannot evaluate using 
our data. Table [T2l theref ore only evalu ates criteria based on the 
dust emission. Following ICrapsi et all evolved cores must have 
central densities (Sec. 15.3. T] recall that this analysis only yields 
rough estimates) exceeding 2.5 • 10 5 crrT 3 (when converting their 
results to our choice of dust properties) and radii at 70% peak in- 
tensity < 4 800 AU. Naively, one would also associate high col- 
umn densities with an advanced evolutionary state. Somewhat 
arbitrarily we thus choose a limit of 0.7 M to highlight those 
starless cores with unusually high column densities in our sam- 
ple. Table [12] also includes candidate VeLLO cores which are 
supposedly young and thus likely to not yet be significantly af- 
fected by YSO formation and evolution. 

This analysis r eveals 5 cores that fulfill both criteria by 
ICrapsi et alj ([2005), i.e., they are small in size and have a high 
central density. For two of these (L1521F CI and L183 C3) this 
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Table 12. Starless cores that fulfill at least one criteria for an 
advanced evolutionary state, as well as VeLLO candidate cores. 
Brackets for VeLLOs indicate that does not fulfill the crite- 
rion for evolved cores after correcting for internal luminosity. 



Core 




«c(H 2 ) 

> 2.5 • 10 5 cm- 3 


770% 

< 4 800 AU 


M4200AU 
> 0.7 M s 


evolved starless cores: 






L1521B-2 


C2 


- 


+ 




L1521-2 


C2 


- 


+ 


_ 


B18-1 


CI 


+ 


+ 


+ 


B18-2 


CI 


- 


+ 


- 


TMC-2 


CI 




— 


+ 


B18-4 


CI" 


+ 


+ 


+ 


TMC-1C 


C3 




— 


+ 


TMC-1 


CI 




+ 


— 


LI 507 A 


C3 




+ 




LI 622 A 


C3 






+ 


L183 


CI 




+ 






C2 


+ 


+ 


+ 




C5 




+ 




LI 082 A 


C2 


+ 


+ 


+ 


VeLLO cores 










L1521F 


CI 


+ 


+ 


+ 


L673-7 


CI 




+ 




LI 148 


CI 




(+) 




L1014 


CI 




(+) 





Notes: a) Numbers hold for the starless P2 peak in B18-4. 



was already shown by Crapsi et al., while our study gives 3 more 
candidate evolved cores matching his criteria (B18-1 CI, B18- 
4 CI, and L1082A C2). Note that these cores are actually the 
only ones with n c (H2) > 2.5 • 10 5 cm 3 and that they also have 
M4200AU > 0.7 M Q . Three further cores (TMC-2 CI, TMC-1C 
C3, and L1622A C3) stand out in their aperture masses, thereby 
suggesting a somewhat evolved state. They do not fulfill the 
Crap si et all criteria, but recent molecular screening of TMC-1 C 
C3 shows that at least one of thes e cores indeed seem s to be 
close to the onset of star formation (Schnee et al. 2007a). 

The c2d MAMBO survey thus reveals 3 previously unknown 
candidates for evo lved starless cores fulfilling the criteria by 
Crap si et alj (|2005), and 3 further ones likely to be evolved on 
basis of their aperture masses. Note, however, that — following 
the above criteria — some of these cores seem to be more ad- 
vanced in their evolution than some candidate VeLL O cores that 
actually already do form stars. This suggests that the Crapsfetal] 
criteria are somewhat biased: they are well suited to select some 
evolved cores, but not all of them. In this respect it is impor- 
tant to keep in mind that their criteria were tuned to select cores 
like the prototype evolved starless core L1544. The analysis pre- 
sented here suggests that evol ved cores with phys ical properties 
different from those derived bv lCrapsi et alj(l2005l) do exist. This 
is not surpri sing though, e.g. given the observed range in stellar 
masses. The lCrapsi et al.l d2005l) study might, e.g., preferentially 
select evolved rather massive cores on the verge to form stars 
of rather large mass. It would then be plausible to find evolved 
cores that are on the verge of forming lower mass stars (i.e., a 
few 0. 1 Mq), but ha ve too low ma ss and column density to be in- 
compatible with the ICrapsi et al.l criteria. In this light, Sec. 15.3.41 
presents some suggestions for further refinements of criteria. 



5.3.3. Star Formation Efficiency 

Given a typic al duration of the main accretion phase of 
order 10 s vr (iBarsony & Kenvonl Il992t iGreene et all Il994t 
iKenvon & HartmanrJ 1 1995b . and expected velocities for the 
growth of the infalling inner part of YSO en velopes of order 
0.2 km s _1 (e.g., for singular isothermal spheres: fShul 1977b . only 
the mass within a few 10 3 AU from a forming YSO are available 
for its formation. The few 10 3 AU radius of this volume are much 
smaller than the subcore sizes of several 10 4 AU observed by us 
(e.g., Fig. [3] [a]). Thus, only a small fraction of a subcore's mass 
is available for the formation of a given YSO. This fraction is an 
estimate of the subcore's efficiency to form a given YSO. The 
latter is a lower limit to the total star formation efficiency, i.e., 
the mass ratio between the dense cores mass and the total mass 
of all YSOs produced by it. 

We use a radius of 4 200 AU for a quantitative evaluation of 
this efficiency, the maximum infall radius derived for the above 
considered YSO and envelope properties. Some 2% to 20% of 
the total mass of a dense core are contained in a 4 200 AU radius 
aperture centered on its brightest peak (Fig.|4][c]). For other av- 
erage velocities at which the infalling region grows, (vi„), and ac- 
cretion times, T accr , the appropriate aperture radius will be larger 
by a factor « (v; n )/(0.2 kms- 1 )-T accr /(10 5 yr). In case of n oc r 1 
density profiles the aperture mass scales linearly with the aper- 
ture radius. To give an example, recent c2d determinations in- 
dicate a class I lifetime of order 5.3 • 10 5 yr (Evans et al., in 
prep). Since YSOs are supposed to accrete during this phase, 
this would imply (using the above assumptions) that a YSO can 
accrete a mass 5.3 • M4200AU and that cores have a star-formation 
efficiency of order of several 10%. 

Long accretion timescales might, however, imply low ac- 
cretion efficiencies on the size scale of the infalling matter, 
if collapse proceeds in free fall. In this case, YSO accretion 
timescales of several 10 5 yr, combined with average aperture 
masses M4200AU = 0.4 to 0.5 M Q (Sec. 15.1.3b . imply that there 
is enough time to accrete several 0.5 M onto a forming YSO. 
Since, however, the average st ellar mass, (M*) is of order 0.4 M 
(e.g., Table 2 in Kroupa 2002), a significant fraction of this aper- 
ture mass does not seem to end up in the forming star. This frac- 
tion increases with increasing accretion timescale; for the above 
collapsing singular isothermal sphere it is {M m - (M*))/M; n , 
where the infalling mass M m = (T accr /10 5 yr) ■ M4200AU and thus 
is of order 3/4 for r acc , w 4 ■ 10 5 yr. The efficiency could, how- 
ever, still be large if collapse does not occur in free fall. 

In summary, our data indicates a star formation efficiency of 
order 10%, with an uncertainty of a factor of several. Such an 
efficiency is of an order of magnitude exceeding estimates for 
the star formation e fficiency of molecular clouds (< 2%; e.g., 
Leisawitz et al. 1989). This might, e.g., indicate that the star for- 
mation efficiency is not only controlled by processes operating 
on the scales of dense cores, but that also by processes of im- 
portance on larger scales. Also, the low star formation efficiency 
on the scales of cores suggests that dense core mass functions do 
not directly map into stellar initial mass functions: core and YSO 
masses will differ, and it is not clear that their ratio is constant. In 
other words, it is not clear which core mass definition is relevant 
to measure the mass of a core that is available for accretion onto 
forming stars. To give an example, here we discuss three mass 
scales (i.e., M tot , A/50%, an d M^ooau), none of which is proven to 
measure the accretable mass. This questions basic assumptions 
made in some studies relating mass functions of stars and cores, 
in particular since the shape and absolute scale of the mass dis- 
tribution depends on the adopted core mass definition (e.g., Fig. 
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[3j panels [b] and [c]). There is evidence though that the shapes of 
some of t hese distributions might well be coupled to the one of 
the IMF dAlves et aLll2007l) . This would then point at a constant 
efficiency. 

5.3.4. Star Formation Ability 

Figure H(d) indicates that there is a fundamental difference be- 
tween cores actively forming stars (i.e., they contain a YSO) 
and those who don't: starless cores in our sample have aperture 
masses of about 0.1 to 2.0 M Q (also see[3][c]), while all young 
(?boi 5 300 K) YSOs have envelope masses > 0.3 M (also see 
Fig.[8j we exclude the VeLLOs, since their physics may be gov- 
erned by different processes, and CB188, since outflow-core in- 
teraction might have affected the core properties [Sec. 15. 1.411 ). 
The naive interpretation of this observation is that dense cores 
must have M4200AU S 0.3 M in order to be able to actively 
form stars. The full story might be a bit more complex, since 
the present YSO aperture masses will deviate from the ones of 
their natal cores at the onse t of star-formation. For example, 
a collapsing SIS dShul 1 1977h has an aperture mass decreasing 
with time, suggesting that the observed YSO aperture masses 
are lower than the initial ones. 

Thus M4200AU ^ 0.3 M Q appears to be a necessary condi- 
tion for active star formation to occur (defined as embedded 
YSOs being present), but it is not a sufficient one, since there 
are many starless cores above the limiting aperture mass (with 
B18-1 P3, B18-4 P2, and L183 P3, all with M 42 ooau> 0.7 M , 
as extreme examples; also see the discussion in Sec. l5.3.2"l >. Thus, 
cores with M4200AU S 0.3 M Q can actively form stars, but they 
do not need to do so. Note, however, that starless cores with 
M4200AU 0.3 M might well already be in a state of collapse 
and be destined to form a star, even when no embedded YSO has 
formed yet. 

In a similar way, a total mass > 2 M appears to be a nec- 
essary condition for active star formation to be possible. Again, 
there are some starless subcores with total masses much higher 
than those for some star-forming ones. 

It is interesting to test this observation against simple models of 
dense core structure. For this we use models of near-isothermal 
pressure-supported spheres. Using boundary conditions appro- 
priate for typical dense cores in the solar neighborhood (a visual 
shielding extinction of 5 mag from a dense core envelope, a cos- 
mic ray flux of 3 • 10" 17 s _1 , and the solar neighborhood inter- 
stellar radiation field), we evaluate the stability of these spheres 
in case of no support form turbulence or magnetic fields (i.e., 
non-isotherma l Bonnor-Ebert spheres; Kauffmann & Bertoldi, 
in prep., and Galli et al. 2002). 

For fixed total mass, respectively mass at half intensity 
(M5o%), these models predict critical values for parameters like 
the peak column density, so that all objects exceeding these can- 
not be stable against collapse unless supported by additional 
non-thermal pressure. These critical limits (or calculations are 
limited to total masses < 5 M & given a lack of molecular cool- 
ing rates for low densities, leading to a break in the boundaries 
drawn) are indicated in Fig. [4] Within the model, they give the 
limiting properties beyond which a dense core can become un- 
stable and collapse to form a star. If these models have any rel- 
evance, the critical limits should therefore manifest in the distri- 
bution of observed actual dense core properties. 

Notably, as shown in Fig.|4](c) and (d), all young YSOs (ex- 
cluding CB188 CI and VeLLO candidates, as discussed above) 
are indeed consistent with having masses in excess of the mass- 



dependent critical aperture mass for pure thermal pressure, and 
the theoretically derived critical masses seem to describe the ob- 
served necessary condition for active star formation. The models 
also describe the necessary condition on the total mass (Fig.|4][a] 
and [c]). 

If dense cores can be described by the above models, a star- 
less dense core of fixed total mass could evolve towards col- 
lapse by approaching the critical aperture mass (starting with 
sub-critical value) as a quasi-hydrostatic sphere, turn into a col- 
lapsing core once exceeding the critical value, and finally be- 
come a core with active star formation. In this picture starless 
cores with properties exceeding the critical ones for pure ther- 
mal pressure could then be collapsing cores. They could also 
be quasi-hydrostatic equilibria with additional non-thermal pres- 
sure, for which the critical values are beyond those for pure ther- 
mal pressure. If the latter was true, the co-existence of starless 
as well as star-forming cores in the same parameter range sug- 
gests that different cores are supported by different levels of non- 
thermal pressure. 

Note that shifts of the masses by a factor 2, either globally for all 
sources, or for starless cores relative to those with YSOs, only 
marginally influence the above discussion. (Such shifts are pos- 
sible because of the uncertainties in opacities and temperatures, 
as discussed in Appendix[A]) Still, our data would hint at the ex- 
istence for necessary conditions for star formation. Strong YSO- 
to-starless core relative shifts of a factor 2 could, however, lead 
to a clear separation of starless and YSO cores in some diagnos- 
tic diagrams. Then, it might turn out that, e.g., all cores above 
some limiting aperture mass do form stars. This would turn the 
aforementioned necessary condition for star formation into a suf- 
ficient one. Strong global and relative mass shifts by a factor 2 
might also make the applicability of the theoretical justification 
for the limiting curve questionable. 

5.4. YSOs and their Evolution 

Given the number of YSO cores, the c2d MAMBO survey is well 
suited to study evolutionary effects in YSO cores. Our analysis 
also includes the properties of some recently discovered VeLLO 
candidates. 

The c2d MAMBO sample contains 5, 14, and 1 YSOs of 
the classes 0, 1, and II, respectively. This is, e.g., in between the 
0-to-I class ratio of 1:10 for Taurus and the Ophiuchus cloud 
dAndre & Montmerlelll994t iMotte et al.|[l998l) and ratios of or- 
der 1 : 1 for Perseus dHatchell et alj2 007) and some isolated cores 
dVisser et alj |2002). Note, however, that our sample is anything 
but unbiased, and thus not ideal to constrain class lifetimes from 
their relative number. 



5.4.1. YSO Luminosities 

To better characterise the central star we use the internal lumi- 
nosity, Li n t, that is due to the YSO instead of the bolometric lumi- 
nosity, i.e., the total radiative power minus that due to interstellar 
heating. This luminosity is, however, taken to be identical to the 
bolometric one, except for some candidate VeLLOs for which 
detailed estimates exist fr om other studies (Bourke et alJ l2006; 
Kauffmann et al., in prep. : lYounget al.ll2004al) . 

The approximation of the internal luminosity by the bolo- 
metric one is justified since the derived bolometric luminosities 
are usually dominated by the power of the embedded source, and 
not by the heating processes in the interstellar medium (like in- 
teraction with cosmic rays or absorption of the interstellar radi- 
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ation field). The latter contribution can be estimated for the ma- 
terial within 4 200 AU from the density peak. A heating power 
0.16 L Q ■ (M420oau/M q ) is required to heat the matter to a typical 
interstellar temperature of 10 K. For M4200AU = 0.4 M , which is 
typical for the YSO cores in our sample (Sec. 15. 1 ."31 . this power 
is 0.06 L Q . 

In some sources, however, the bolometric luminosity derived 
by us is dominated by the power due to interstellar heating pro- 
cesses. For dust temperatures of at least 10 K this power exceeds 
the inferred bolometric luminosities of the YSOs in LI 52 IF 
PI, L673-7 PI, and L1251A P3 and P5 (also depending on 
whether relying on IRAS or Spitzer data), also demonstrating 
that our method to derive Lboi can miss part of the actual radia- 
tive power. This is also reflected in these sources' unusually large 
submillimetre-to-bolometric luminosity ratios. In such cases the 
bolometric luminosity is not a good measure of the power input 
by the young star. This leads to uncertain heating corrections in 
mass estimates of YSO envelopes (Sec. 13.2.3b . The related un- 
certainties in the mass estimates are, however, not significant, 
given that the low derived luminosities imply small corrections. 

5.4.2. Envelope Masses and Luminosities 

YSOs are expected to evolve in luminosity and envelope proper- 
ties during the star formation process. Three basic properties to 
characterise a YSO are its age, as believed to be roughly mea- 
sured by the bolometric temperature or the infrared class, the lu- 
minosity, and the amount of matter surrounding the star. The lat- 
ter two should depend on the age since they, e.g., depend on the 
accretion rate and the stellar mass and radius, which all evolve 
with time. Figure [6] shows the relation between luminosities and 
envelope masses. The latter are approximated by the mass within 
4 200 AU radi us from the star since t his property has been stud- 
ied in the past dMotte & Andrell200ll) . 

Surprisingly, no clear evolutionary trend is obvious in this 
diagram, except for the significant decrease in M4200AU towards 
the class II phase; the crowding of VeLLO candidates at low lu- 
minosities is due to their definition. Class and class I sources 
do essentially occupy the same region in the parameter space, 
with the class sources possibly having for given internal lumi- 
nosity slightly higher aperture masses than class I sources. 

In particular, the classes and I are no t separ ated by the 
class boundary suggested by Mott e & Andrei (1200 ll) . M4200AU ~ 
O.lM Q (L/L )°- 6to 10 . This law reflects the conceptual defini- 
tion of class sources a s YSOs which hav e accreted less than 
half of their final mass (lAndre et al.l[l9 93). Thus the envelope 
mass in class sources, as approximated by M4200AU, must ex- 
ceed the stellar mass. De pending on the e xact form of the YSO 
mass-luminosity relation. TMotte & Andrei suggest this boundary 
to b e given by the abov e relation (see I Andre & Montmerldl 19941 
and I Andre et al.ll2000l for further details). This boundary is in- 
deed observed to well separate class and class I sources in 
the Ophiuchus mol ecular cloud complex and some isolated cores 
dVisser et~ al. 2002), bu t it fai ls to do so for sources in the Tauru s 
cloud dMotte & Andrell2001h and Perseus (lHatchell et al.ll2007l) . 

To be precise, within the observational errors all class 
sources in the c2d MAMBO sources do indeed have aper- 
ture masses exceeding 0.1 M (L/L ) 06 to , just as expected. 
However, the class I YSOs in our survey do not generally 
fall short of the conceptual luminosity-dependent mass limit 
for class sources. In this respect they are different from the 
aforementioned class I sources in Ophiuchus that fulfill this 
limit. For comparison , like ours, a lso class I YSOs in Taurus 
dMotte & Andre1l2001l) and Perseus dHatchell et al.ll2007l) do not 
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Fig. 6. The relation between the internal luminosity and the mass 
within a peak-centered 4 200 AU radius aperture for YSO cores. 
The internal luminosity is the radiation due to the embedded 
star, i.e., the bolometric luminosity minus the power from in- 
terstellar heating. Markers of different shade and shape refer to 
YSOs belonging to different infrared SED classes assigned on 
basis of their bolometric temperature (see legend). The dashed 
lines indicate the conceptual boundary between class sources 
and objects in later evolutionary stages for a spectrum of stel- 
lar mass-luminosity relations, M 42 ooau ~ 0.1 M (L/L ) 06 t0 I . 
Interestingly, class and class I sources covered by the c2d 
MAMBO survey are not separated by the conceptual class 
boundaries. 



conform to this conceptual boundary. This suggests that the stars 
in our sample reside in regions that are more similar to regions 
of distributed star formation, like Taurus, than to regions of clus- 
tered star formation, like Ophiuchus. 

The failure of the conceptual class boundary to separate 
sources of class and I and the absence of clear evolutionary 
trends between them in Fig. [6] suggests that either the concept 
of infrared classes, our general understanding of YSO evolu- 
tion as a linear sequence of states in the Lmt-A^4200AU param- 
eter space, or both, need some revision. Also, M4200AU might 
simply not be well suited to gauge the mass available for ac- 
cretion. This would, however, require accretion to proceed in a 
manner very different from the collapse of self-similar isother- 
mal spheres: here steadily M4200AU decreases with time, and in- 
dependent from the final size of the infalling envelope (even if 
eventually > 4 200 AU) M4200AU would at least generally probe 
the mass evolution of the envelope. Whatever, the fact that sim- 
ilar studies of different regions find different relations between 
YSO and envelope properties suggests that the environment in- 
fluences these relations. 



5.4.3. Bolometric Temperatures and Luminosities 

Figure Q presents the bolometric temperatures and luminosities 
for sources in the c2d MAMBO survey. It shows that for a given 
star the luminosities are higher and the bolometric temperatures 
are lower when these properties are inferred from IRAS data 
instead of Spitzer observations. This is just as expected, given 
IRAS' better spectral coverage of the dust emission peak at sub- 
millimeter wavelengths. 
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Fig. 7. Bolometric luminosities and temperatures as derived 
from IRAS and Spitzer (markers of different shade) data. VeLLO 
candidates are highlighted by circles around the markers, and the 
YSO class regimes are delimited by vertical dashed lines and la- 
bels. Solid lines indicate characterizations of the same star based 
on different data. For a given source the IRAS data usually in- 
dicates lower bolometric temperatures and higher luminosities 
than the Spitzer data, just as expected given IRAS' better spec- 
tral coverage at long wavelength. Note the range in luminosities 
for a given class. This indicates that the accretion rates, masses, 
radii, or all these properties largely vary within the sample, and 
also a given class, if the luminosities are largely due to accretion. 



Fig. 8. Bolometric temperatures and mass within 4 200 AU ra- 
dius for the YSOs in the c2d MAMBO sample. Candidate 
VeLLOs are highlighted by circles around the markers, and the 
YSO class regimes are indicated by vertical dashed lines and 
labels. The general trend of increasing bolometric temperature 
with decreasing aperture mass might be an evolutionary effect. It 
could, however, also indicate that low bolometric temperatures 
are partially due to reddening by extinction and re-emission at 
longer wavelength, since both should scale with aperture mass. 
VeLLOs notably have the lowest bolometric temperature for 
given aperture mass, but this is likely simply reflects their defi- 
nition as objects of low luminosity. 



This diagram shows that on average Lboi increases with in- 
creasing rboi. This trend is in particular pronounced when look- 
ing at the lower envelope of the data points. However, it is likely 
that this is to a significant extent an observational bias, since 
faint sources of high bolometric temperature have lower mil- 
limetre intensities and are thus less likely to be included in our 
MAMBO-based source selection. For the same reason in general 
sources with high bolometric temperatures are less likely to be 
included. This is at least partially responsible for the sharp drop 
in the number of sources at > 300 K. 

Class sources are found to have luminosities spanning 
across more than two orders of magnitude. This demonstrates 
the diversity of star formation activity in the c2d MAMBO 
survey. Note that, if sources of similar T\, Q \ have similar ages 
and have luminosities dominated by accretion luminosity (i.e., 
Lint ~ GM+M^IR*, where G, M+, 7?*, and M accr are the con- 
stant of gravity, the stellar mass and radius, and the accretion 
rate), then the class sources in the c2d MAMBO survey must 
differ significantly in their mass-to-radius ratios, their accretion 
rates, or both. Enoch et al. (in prep.) find similar luminosity 
ranges for YSOs in Perseus, Serpens, and Ophiuchus. 



However, the large column densities suggested by large aper- 
ture masses might lead to low bolometric temperatures because 
of reddening: photospheric YSO emission might be absorbed at 
short wavelengths and be re-emitted by the envelope at longer 
wavelengths. In that case bolometric temperatures would not un- 
ambiguously characterize YSO ages, but would partially be due 
to reddening not related to the time evolution of the source. Our 
sample is, however, too small to evaluate this issue comprehen- 
sively. 

Within our survey, VeLLO candidates appear to have for 
given aperture mass the lowest bolometric temperature observed. 
However, this likely reflects the definition of VeLLOs as objects 
of low internal luminosity, since for fixed millimetre flux (or 
A^4200Au) the bolometric temperature decreases with decreasing 
luminosity. 

Note that all the "usual" young YSOs with r bo i < 300 K 
have M4200AU > 0.3 M , just as already discussed in Sec. 15. 3. 41 
Also note the aperture masses of M4200AU Z 0.5M Q for the YSOs 
in L1228 PI and Bern48 PI are unusually large, given their bolo- 
metric temperatures of > 300 K. It is not clear how such sources 
of a relatively late evolutionary stage can maintain such massive 
envelopes, if our mass estimates are correct. 



5.4.4. Envelope Masses and Bolometric Temperatures 

Figure [8] presents bolometric temperatures and aperture masses 
for the c2d MAMBO sources. This reveals that on average the 
bolometric temperature decreases with increasing aperture mass. 
This could be an evolutionary trend. The coolest sources are 
supposedly the youngest sources, which are expected to be sur- 
rounded by significant amounts of circumstellar matter. Our ob- 
servations are thus consistent with expectations. 



5.4.5. Assignment of Infrared Classes 

The above discussion often assumed that the bolometric temper- 
ature uniquely parameterizes the YSO age. This is not obvious, 
and possibly not true. We briefly summarize hints on this issue 
from the c2d MAMBO data. 

Figure [9] presents the relation between bolometric tempera- 
ture and SED slope near 10 fim wavelength. The slopes derived 
from IRAS and Spitzer data (ffj^™ and , respectively) do 
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Fig. 9. SED slope versus the bolometric temperature. See Fig. [7] 
for the symbols used. Note that IRAS data usually implies SED 
slopes larger than, and bolometric temperatures smaller than, 
those derived from Spitzer data, just as expected from IRAS' 
better spectral coverage at long wavelengths. Note that class as- 
signments on basis of bolometric temperature and spectral slope 
do agree. 



Fig. 10. Submillimetre-to-bolometric luminosity ratio versus the 
bolometric temperature. See Fig.|7]forthe symbols used. Vertical 
dashed lines indicate ranges for YSO classes assigned on basis 
of 7boi, while the horizontal dotted line (described by inclined 
labels) gives the boundary for the class assig nment based on the 
luminosity ratio following An dre et ai](ll993l) . It's disagreement 
with Tboi-based class assignments is likely caused by positive 
biases in our evaluation of the submillimetre luminosity. 



slightly differ in their definition (see Sec.|4]l. In 4 out of 5 cases 
> for sources detected by both missions. Most 

sources in our survey have spectral slopes of 0.0 < a < 1.5. We 
find a general trend of increasing a with decreasing T\, Q \ across 
the whole sample. No such trend is, however, visible when look- 
ing at class I sources only (class sources are usually not ex- 
pected to have particular values of a). Note that class assign- 
ments on basis of bolometric temperature or spectral slope are 
consistent with each other within our sample. Given the expected 
uncertainties, all sources in our sample are consistent with hav- 
ing Tboi < 650 K (i.e., being a class or class I source) when 
a > (note, however, that classification schemes never required 
class to have a > 0), and 650 < r bo i/K < 2880 (i.e., be- 
ing a class II source) when > a > —2, just as postulated in 
class definitions based on 7b i and a dChen et al.l I 1995b lLadal 
|1987|) . Similar trends are found in the c2d SED studies of the 
Ophiuchus, Perseus, and Serpens molecular cloud complexes by 
Enoch et al. (in prep.). 

Also note that all IRAS sources not detected by MAMBO 
(Table|7]l are consistent with having negative spectral indices and 
therefore belonging to late YSO stages (if they are YSOs and 
not background galaxies). This fits well into a picture in which 
YSOs with negative spectral slopes are evolved and not anymore 
surrounded by dust. 

Our submillimetre-to-bolometric luminosity ratios (Fig. [TO) 
also see Tables [8] and |9]l and Tboi-based class assignments are 
not well in accord with the original in itial observa t ional crite- 
rion to identify class sources given bv lAndreet all (119931) . i.e., 
a luminosity ratio » 0.005; using the latter criterion, all but 
one of our YSOs would be in the class category. Note, how- 
ever, that the submillimetre part of our SEDs is the worse con- 
strained one and that the luminosities are uncertain by a factor 
2(Sec. l4~2l >. Anyway, submillimetre luminositie s of 0.01 to 0.1 
for ? boi ~ 70 K have been observ ed before dYoung & Evansl 
120051) and also IVisser et"aT] d2002l) found that essentially all 
YSOs in their sample would be of class if t he class is as- 
signed based on a luminosity ratio > 0.005 (see lHatchell et alj 



120071 for a further example; note that both studies do, like us, 
have SEDs not well sampled at submillimetre wavelengths). The 
role of the submillimetre-to-bolometric luminosity ratio for clas- 
sification purposes thus needs to be further studied with better 
wavelength coverage. 

Also, the bolometric temperature might not uniquely param- 
eterize age, as the aforementioned trend of decreasing bolomet- 
ric temperature with increasing aperture mass might indicate. 

It has also been speculated that actual class I sources have 
SEDs like class sources if the circumstellar disk-e nvelope 
structure is seen edge-on (e.g., IWhitnev et al.l 120031 for an 
overview). Interestingly, the class source in IRAS05413 P3 
appears to have such a particular geometry; here the outflow 
has an estimated inclination of > 85° with respect to the 
line of sight, suggesting that a YSO disk is seen nearly edge- 
on (IClaussen et al.l 1 1998b . This geometry might explain why 
IRAS05413 P3, being a class sou rce marginally v i olatin g the 
conceptual mass-luminosity limit bv lMotte & Andrei (1200 ll) . has 
an extreme position in Fig. [6] 

In summary, class assignments based on 7boi are consistent with 
those based on a, and either is fine to classify sources of a rela- 
tively late evolutionary stage. One needs to employ the bolomet- 
ric temperature though if intending to study deeply embedded 
sources, for which the spectral index criteria do not hold. The 
7boi thus appears to be more versatile, but requires a well sam- 
pled SED at long wavelengths. 

In contrast, assignments based on the OT42ooAu-to-L; nt ratio 
are in conflict with those based on 7b i- A Tboi-based classi- 
fication does however lead to a larger spectrum in evolution- 
ary states: 7boi indicates class and class I sources, while the 
'M4200AU-to-Li nt ratio suggests that our sample essentially only 
contains class objects. Since the latter seems somewhat un- 
realistic, the larger range in assigned classes suggests that the 
bolometric temperature is more useful for classification. 

Similarly, the submillimetre-to-bolometric luminosity ratio 
assigns all but one YSO to the class 0. Again, 7b rbased classi- 
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fication appears superior to the one based on luminosity ratios. 
This might, however, in part be due to uncertainties in our SEDs 
at wavelength 100 to 1000 fim and should not be generalized. 
Still, for data like ours luminosity ratios are barely useful. 

Anyway, it is not clear which existing YSO class assign- 
ments best characterize the evolut ionary ages of YSOs. In this 
respect it might be useful (see, e.g.. lRobitaille et a l. 2006) to dis- 
tinguish between YSO classes, that describe and organize obser- 
vational properties such as the bolometric temperature, and evo- 
lutionary stages, that characterize physical properties such as the 
presence of envelopes and disks. 

5.4.6. YSO Offsets 

The spatial resolution of MAMBO and Spitzer allows to test 
whether YSOs actually reside in the very center of their natal 
dense core. The discovery of offsets would be very interesting, 
since present YSO evolutionary scenarios do not predict them. 
As shown in Table [T3l we do, however, not find offsets at a sig- 
nificant level. 

For this search the YSO positions are taken to be those of 
Spitzer counterparts listed in Table|9] Otherwise the positions of 
the 2MASS counterparts given in Table [8] are used. For YSOs 
with non-c2d Spitzer data only, Gaussian curves were fitted to 
the Spitzer images (see notes to Table [T31 , The dust intensity 
peak positions are derived as the intensity-weighted mean posi- 
tion of the pixels above 90% peak intensity. The offsets are given 
in Table [T3| The uncertainty in the dust emission peak position 
due to noise is derived from Monte-Carlo experiments with ar- 
tificial noise, as explained in Sec. 13.21 The typical pointing un- 
certainty of 3" has to be further added to derive the typical total 
dust emission position uncertainty. 

We do not find significant YSO offsets in our sample. In a 
first cut we find two YSOs with offsets exceeding the expected 
uncertainty by a factor 2. However, for one of them (L1251A 
P3) the pointing corrections were up to 5'.' 9. For the other one 
(LI 172 A PI) the position offset exceeds the uncertainty by a fac- 
tor 2.2; this is marginally significant, and should be followed up 
at higher resolution, but not conclusive evidence for a YSO leav- 
ing its natal core. Our data thus indicates that generally YSOs 
reside in the centers of their envelopes, with core- YSO relative 
motions not e.g. significantly affecting the accretion histories. 
This is in accord with studies of core-YSO-offsets in Perseus 
( J0rg ensen et al.l 120071 : Enoch et al., in prep.) and Ophiuchus 
(J0rgensen et al., in prep.). 

Note, however, that our data does not reveal the offset be- 
tween the L1014 PI colu mn density peak and the embedded can- 
didate VeLLO found by iHuard et all f2006) on basis of extinc- 
tion maps, although th eir offset of 10" to 15" should be easily 
detected by MAMBO. IHuard et"aT] speculate that this might be 
because of internal heating of the L1014 core, which is supposed 
to shift the dust emission peak away from the column density 
pea k and towards t he embedded source (see Young et aL 2004a 
and lWu et alj2007l for SCUBA and SHARC n data on L1014, re- 
spectively). Their extinction map is, however, not well sampled 
towards the VeLLO position, which introduces uncertainties to 
their map. 

5.5. Formation and Evolution of VeLLOs 

The c2d MAMBO survey provides an opportunity to investigate 
how candidate VeLLO natal cores are different from starless and 
YSO cores. For the first time it allows a direct comparison of 



Table 13. YSO offsets from dust emission peaks. For each dust 
emission peak, and related subcore, associated with a YSO listed 
in Tables|8]and|9]we give the chosen reference source (preferen- 
tially seen by Spitzer) and the observed offset and position angle 
(i.e., position of star w.r.t. the dust emission peak). 



Field 


Peak 


Position Reference 


Offset 






(see notes) 


arc sec 


L1521F 


PI 




5.0 


±0.6 


B18-1 


rP5i 


SSTr2d T043215 4+242859 


1.1 


±0.6 


B18-4 


Pi 


IRAC 


2.2 


±0.6 


TMC 1 - 1 C 


P3 


SSTr2d T044138 8+255627 


3.3 


±2.5 


TMC-1 


P2 


SSTr2d T0441 12 7+254635 


3.1 


±0.6 


JRAS05413 


P2 


MIPS 


3.8 


±0.5 




P3 


MIPS 


3.6 


±0.3 


CB188 


PI 


SSTc2d J192014.9+1 13540 


2.3 


±0.7 


L673-7 


PI 


SSTc2dJ192134.8+l 12123 


2.1 


±0.9 


LI 148 


PI 


SSTc2d J204056.7+672305 


7.6 


±2.7 


L1082A 


PI 


2MASS J20531346+6014425 


2.7 


±0.9 


L1228 


PI 


SSTc2d J205712.9+773544 


0.8 


±0.1 


Bern48 


PI 


SSTc2d J2059 14.0+782304 


0.8 


±0.3 


L1172A 


P2 


MIPS 


7.3 


± 1.4 


LI 177 


P2 


IRAC 


2.9 


±0.1 


L1021 


PI 


2MASS J21212751+5059475 


2.5 


±0.8 


L1014 


PI 


SSTc2d J212407.5+495909 


2.2 


+ 1.3 


L1251A 


P3 


SSTc2d J222959.5+751404 


8.8 


± 1.3 




P4 


SSTc2d J22303 1.8+75 1409 


1.2 


±0.2 




P5 


SSTc2d J223 105.6+75 1337 


4.3 


±0.8 



Notes: Position for B18-4 PI (average of Gaussian fits to all IRAC 
bands) is 04:35:35.4, 24:08:20.2 (J2000.0); for IRAS05413 P2 (from 
fits to MIPS band 1) it is 05:43:46.264, -01:04:43.68; for IRAS05413 
P3 (from MIPS band 1) it is 05:43:51.367, -01:02:52.76; forL1172AP2 
(from MIPS band 1) it is 21:02:21.1 18, +67:54:20.17; and for LI 177 P2 
(from IRAC band 4 only, where the significant scattering is supposedly 
lowest) it is 21:17:38.650, +68:17:33.27. 



the dust emission properties of all three kinds of objects. We 
here discuss whether, and if yes how, VeLLO cores differ from 
starless and YSO cores. 

5.5.1 . Masses and Total Radii of VeLLO Cores 

Figures[3]and|4]show that VeLLO candidate cores cannot be dis- 
tinguished from starless cores in terms of mass, column density, 
and total effective radius. In particular, three of the VeLLO cores 
have column densities, masses within the 4 200 AU aperture, and 
masses within the 50% intensity contour — all properties naively 
associated with a good chance for star formation — much lower 
than the starless peaks B18-1 P3 and B18-4 P2. Note that this 
observation is in conflict with some naive expectations for prop- 
erties of starless and YSO cores: a picture of dense core evolu- 
tion in which any dense core begins to actively form stars once 
it exceeds some fixed limit in mass, density, or both, is inconsis- 
tent with our data. Section |5 .3 .41 discussed this and other data in 
the context of necessary, but not sufficient, conditions for active 
star formation. 

Note that 3 out of 4 VeLLO candidates are consistent with 
fulfilling these conditions (i.e., in Fig. |H only the candidate 
VeLLO in LI 148 resides in a region entirely devoid of normal 
YSOs, and thus also misses the theoretical star formation crite- 
rion by a large margin). In this respect the notion that VeLLOs 
reside in cores not expected to be able to form stars might be 
wrong. However, it might be that some VeLLOs can conceptu- 
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ally be characterized as YSOs forming in cores that only barely 
fulfill the necessary conditions for active star formation, as ob- 
served for 3 out of the 4 candidate VeLLOs in our sample. The 
latter might suggest a special route of star formation for some 
VeLLOs. 

The VeLLO candidate in LI 148 PI poses a real challenge to 
present star formation scenarios, since it does not fulfill the nec- 
essary conditions for star formation to occur. This source might 
thus either hint at a particular mode of star formation, where the 
conditions derived in Sec. l5.3.4l do not apply, or this YSO (recent 
IRAM Plateau de Bure imaging indicates either a jet or a disk; 
Kauffmann et al., in prep.) is in an advanced evolutionary stage 
in which the natal core's structure does not resemble the initial 
one. 



5.5.2. Peak Radii of VeLLO Cores 

Most VeLLO candidate cores have peak density profiles flatter 
than g oc r~ 2 , since their effective radius at 70% peak intensity is 
(with the exception of the candidate VeLLO in L673-7 PI) sig- 
nificantly larger than the one of a g oc r~ 2 density profile smeared 
by a 20" beam (i.e., 10'.'7). It is beyond the scope of the present 
paper to explore whether this flattening is because the density 
profile becomes similar to p oc r~^ 2 , a s found for infalling en- 
velopes around YSOs (see Evans 2003 for a review), or is due 
to a region of spat ially constant density, as f ound for starless 
dens e cores (e.g. , Ward-Thompson et alJll999h . Note, however, 



which Crapsi et al. (2005) suggested as an indicator of an ad- 
vanced evolutionary stage of a core. This might bear important 
hints on the nature of VeLLOs. 

5.5.3. VeLLO Offsets from Column Density Peaks 



that lCrapsietaU J2004I) show that MAMBO maps for L1521F 
are consistent with a functional form for the central flattening as 
found for starless cores. The high reso lution extinction data for 
L1014 obtained by Hua rd et all (H006) has not been analysed in 
such a fashion. 

Compared to starless cores, candidate VeLLO natal cores do, 
however, have unusually steep dust emission intensity profiles, 
as e.g. probed by the effective radius of a dust emission peak 
at its 70% intensity level. For given aperture mass, VeLLO cores 
(with the possible exception of the uncertain radius of LI 148 CI) 
do have smaller effective radii than most starless cores (Fig. [5]), 
even when correcting for steepening due to internal heating by 
the VeLLO (for this we first use the luminosity-dependent beam- 
averaged dust temperature [Eq. [8] to predict the peak intensity 
for dust of 10 K, and then evaluate the effective radius at 70% of 
this intensity). 

The unusually steep dust intensity profiles of candidate 
VeLLO natal cores thus indicate unusually steep column density 
profiles. This distinguishes them from most starless cores in our 
sample. Steep radial profiles in the density, and therefore also in 
the column den sity, are indeed expected fo r dense cores tempo- 
rally near (e.g., iMcKee & Hollimanlll999l) or after (see iLarsonl 
2003 for a review) the onset of gravitational collapse. Density 
gradients in VeLLO cores are consistent with this theoretical pic- 
ture. VeLLO candidates also seem to have unusually large cen- 
tral densities, when compared to starless cores (Fig. 0. This is 
again in line with the above theoretical expectations. 

The steeper density profiles and higher central densities of 
candidate VeLLO natal cores suggest that they are temporally 
closer to the onset of gravitational collapse than most starless 
cores. It thus appears that VeLLO cores are structurally different 
from most starless cores in that the ir density structure is more 
evolved. This was first suggested bv lHuard et alj (120061) . on ba- 
sis of surface-to-center density contrasts for L1014 PI that ex- 
ceed those of previously studied starless cores. Also, all candi- 
date VeLLOs (possibly excluding LI 148 PI) have effective radii 
at the 70% dust emission peak intensity level below 4 800 AU, 



iHuard et al.1 (120061) found that the candidate VeLLO in LI 014 
is offset by 10" to 15" from the apparent natal column density 
peak, L1014 PI. If L1014-IRS does indeed drift away from its 
natal core, and therefore away from the mass reservoir avail- 
able for accretion, this migh t explain the low r ate of accretion 
onto L1014-IRS inferred bv lBourke et al.1 d2005l) . Low accretion 
rates could in turn explain the low luminosities of VeLLOs. They 
would also suggest low, and possi bly substellar, final masses for 
VeLLOs with low accretion rates (Bo urke et al.l20 05). However, 
as already discussed in Sec. 15.4.61 we do not find any significant 
offsets. Our data might not rule out such offsets though, since 
the MAMBO data might not be suited to search for offsets (Sec. 
l5A6l >. 

6. Summary 

We present a dust continuum emission imaging study of a 
comprehensive sample of isolated starless and YSO-containing 
dense molecular cores (including 4 candidates of the enigmatic 
VeLLOs), the c2d MAMBO survey. This survey is technically 
different from previous ones in that it is 

- more sensitive to faint (few mJy per 11" beam) extended 
(> 5') emission than in previous surveys thanks to modern 
instrumentation and a newly devised data reduction strategy 
(Sec. l2.3.TT i. resulting in 

- very rich structure compared to existing surveys, since we 
detect faint extended emission from sources otherwise only 
detected in their most prominent intensity peaks. 

We combine our MAMBO data with Spitzer and IRAS data to 
fully exploit our dust emission maps (Sec. This also reveals 
4 previously unreported YSOs (TablefTTTi. including one VeLLO 
candidate and 2 class sources. Our analysis (Sec.|5]l first deals 
with the general structure and evolution of dense cores (Secs. lBTTI 
tol53T>: 

- characteristic values for the mass within peak-centered aper- 
tures of 4200 AU diameter are similar to those of singu- 
lar isothermal spheres, suggesting that cores prefer a near- 
critical state (with respect to gravitational collapse; Sec. 
l5T3l andFig.l3lrcl); 

- subcores (i.e., extended substructure within cores) are usu- 
ally elongated, implying that they are not hydrostatic objects 
supported by isotropic pressure (Sec. 15.1.41 and Fig.|3][d]); 

- rough density estimates indicate central H2 densities be- 
tween 10 4 and 10 6 cm 3 while the density frequency dis- 
tribution suggests that cores in our sample are not in free fall 
(Sec.[533]andFig.[5]i; 

- comparison with criteria by ICrapsi et al.l d2005l) reveals 5 
candidate evolved cores (plus 3 candidates of lesser quality), 
out of which 3 are reported here for the first time, but the 
failure of these criteria to highlight candidate VeLLO cores 
with active star formation casts some doubt on them (Sec. 
E321 

- rough estimates for the dense core mass available for accre- 
tion onto a forming YSO yield values of order 10% of the 
dense core total mass (with an uncertainty of a factor of sev- 
eral), suggesting that the star formation efficiency is of the 
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same order (which questions studies relating the mass func- 
tions of cores and stars; Sec. 15.3.31 and Fig.|4][c]); 

- there are necessary (but not sufficient) conditions which a 
dense core must fulfill in order to be able to actively form a 
star (i.e., a core fulfilling these can contain a YSO, but does 
not need to), which can be understood in the framework of 
quasi-static dense core models (Sec. 15.3.41 and Fig. |4j»; and 

- diverse biases can and do introduce spurious features into the 
data, requiring care when interpreting the observations. 

We then specifically turn to young stellar objects in our sample, 
where we in particular compare various methods to assign SED 
classes (Sec. l5.4b : 

- class assignments based on the bolometric temperature are 
in conflict with those based on the comparison of YSO lumi- 
nosities and envelope masses (Sec. 15.4.21 and Fig.|6]i; 

- classes assigned on the basis of 7boi and SED slopes agree 
well, i.e., the data is consistent with a % for Tboi S 650 K 
(Sec. l5A5l andFig.l9T): 

- among the class assignment methods used in the present 
study, only the one based on bolometric temperatures im- 
plies a reasonable range of evolutionary stages in our sample 
and is also well applicable to deeply embedded sources (Sec. 
l5A5li : 

- disagreement of the 7boi-based class assignments with those 
from YSO luminosities and envelope masses suggests that 
the environment has an impact on YSO evolution since 
these assignments agree in other regions (like the Ophiuchus 
molecular cloud; Sec. 15 .4.21 and Fig.[6]i; 

- there are no significant offsets between YSOs and their natal 
dense cores, implying that core-star relative motions do not 
have a significant influence on YSO evolution via the accre- 
tion history (Sec. 15.4.61 and Table [T3l: and 

- assuming that the bolometric temperature roughly character- 
izes YSO ages, the large range of bolometric luminosities 
(two magnitudes) in early YSO stages implies a large vari- 
ation of mass-to-radius relations, accretion rates, or both, in 
our sample. 

The combined information on cores in general and YSOs in par- 
ticular can then be used to assess the relevance of some charac- 
teristics of candidate VeLLO natal cores (Sec. 15.5b : 

- 3 out of 4 VeLLO candidate cores are consistent with fulfill- 
ing the above necessary conditions for active star formation, 
suggesting that the notion that VeLLOs form in cores not ex- 
pected to be able to form stars is wrong, but some do only 
barely fulfill the requirements for active star formation, sug- 
gesting that some VeLLOs might form in cores just barely 
able to produce stars, possibly via a particular mode (Sec. 
13311 ): 

- candidate VeLLO natal cores have density profiles steeper 
than and central densities larger than most starless cores, 
suggesting that VeLLO candidate cores are structurally dif- 
ferent from starless ones (Sec. 15.5.2k and 

- given our failure to detect core-VeLLO offsets, low accre- 
tion rates because of core-star motion as an explanation for 
the low VeLLO luminosities are not consistent with our ob- 
servations (Sec. 15.5.31 ). 

In this respect our study — the first one to compare candidate 
VeLLO cores based on a homogeneous dataset — helps to sig- 
nificantly further our understanding of the nature of VeLLOs. 
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Fig. 1. Maps of the sources in the c2d MAMBO survey after smoothing to 20" resolution. Coordinates are given for J2000.0. 
The sources are roughly ordered in sequence of increasing right ascension. Solid and dotted black contours are drawn for positive 
intensities and intensities < mJy beam -1 , respectively. Contours start at zero intensity and are spaced by ±3 mJy per 11" beam; 
this corresponds to the noise level in the map center times a factor m 2 to 3 (Table|2]i. Above 15 mJy beam -1 , indicated by a thick 
solid contour, the contour spacing is 5 mJy beam -1 . Gray shading also reflects these intensity variations. White contours indicate 
subcore boundaries, while crosses indicate dust emission peak positions. Horizontal bars indicate the angular scale corresponding to 
a scale of 0.2 pc w 40 000 AU at the adopted core distances. Circles and squares mark the positions of all infrared stars detected by 
IRAS, and Spitzer stars also detected as MAMBO peaks, respectively. If detected by several missions, IRAS Point Source Catalogue 
(PSC) positions are used instead of those from the Faint Source Catalogue (FSC), and Spitzer positions are used instead of those 
listed from IRAS. The IRAS PSC source 20410+6710 in LI 148 is not marked since missing Spitzer counterparts indicate that this 
source is an artifact. [See the alternative URL in the astro-ph comments, and the journal article, for higher resolution maps.] 
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Fig. 1. continued. 
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Fig. 1. continued. 




Fig. 1. continued. 




Fig. 1. continued. 



Table 3. General data on subcores. Bracketed and starred subcore designations indicate subcores with uncertain properties due to artifacts and 
subcores with associated YSOs or VeLLOs, respectively. For every subcore the table lists the position of its "center of gravity", the subcore area 
and effective radius, the integrated flux density and mass for the whole subcore as well as its area above 50% of its peak intensity, and the maximum 
signal-to-noise ratio within the subcore boundaries. 
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_L / . I iy _!_ u.uu 


1 6 1 1 + 07 

I W .11 _!_ vy . vy / 


8.44 ± 





.20 


7.77 ± 0.19 


43 


B18-4 


Cl* 


04 35 36 

\Jl ~J «y V_7 


.7 


+24 09 05 


.0 


24.1 


166.1 


23260 


1 2 54 + 1 1 


7 05 + 34 

/ • vy _!_ vy . ^7 r 


4.03 ± 





.17 


2 27 + 1 1 


42 


B18-5 


Cl 


04 35 50 

U i »y »y »y vy 


.4 


+24 09 20 

i z— v./ y z— vy 


.6 


11.1 


112.6 


15768 


2 92 + 04 

Z— . y z— 1 v./ . v./ 


2 69 + 04 

i-.uy _!_ \y.\y~ 


1.63 ± 





06 


1 51 + 05 


22 


TMC-1C 


Cl 


0441 21. 


,5 


+26 02 55 


.2 


14.3 


128.0 


17921 


4.19 ± 0.05 


3 86 + 04 

. vy vy _!_ u.ut 


2.38 ± 





08 


2 1 9 + 08 

^— . 1 y _!_ vy . vy vj 


22 




C2 


04 41 31 


7 


+26 04 06 

T^z— Ay ut \y Vy 


.5 


16.9 


139.4 


19511 


a 33 + 06 

UpJ y -!- >y • >y *y 


5 84 + 05 

«y . vy ~ _!_ vy . vy ^/ 


3.06 ± 





,25 


2 82 + 23 

Z— . V.J Z— 1 V./ . z— ^ / 


28 




C3 


04 41 35, 


9 


+26 00 54 

y^z- vy vy vy >y~ 


1 


26.7 


175.0 


24505 


11 A3 + o 07 


10 71 + 07 

j. \y . / j. _!_ vy . vy / 


7.00 ± 





.30 


6 45 + 27 

vy . 1 »y _!_ \y . / 


28 


TMC-1 


Cl 


04 40 59, 


.6 


+25 53 46 


j 


5.9 


82.0 


11479 


1 64 + 05 


i 51 + 05 


0.95 ± 





.10 


88 + 09 

vy . vy vy _!_ vy . vy y 


10 




C2* 


04 41 13 


,3 


+25 46 31 

1 Z— „ / 1^ V ) ~ ' 1 


.7 


14.5 


128.9 


18046 


5 05 + 09 


2 80 + 05 

^- . vy vy _!_ \y.\y,y 


0.51 ± 





.08 


28 + 04 

vy vy _!_ u.ut 


25 




C3 


04 41 16 


,0 


+25 49 1 1 


.6 


32.7 


193.5 


27092 


13 82 + 1 1 


12 73 + 10 

± _ * / *j _!_ vy . \_ vy 


7.77 ± 





62 


7.16 ± 0.57 


21 




C4 


04 41 34 


9 


+25 44 20 

1 Z- l l z— V./ 


.7 


17.2 


140.2 


19634 


6 40 + 09 

*y . i \y _!_ \j ,\j j 


5 90 + 09 

-J • _!_ vy . vy y 


3.85 ± 





.19 


3 55 + 17 

y -!- \y . j. / 


13 




C5 


04 41 42. 


.3 


+25 41 17 


,4 


13.3 


123.6 


17302 


5 72 + 1 


5 27 + 09 

~->.— 1 _!_ vy . vy y 


3.83 ± 





.23 


3 53 + 21 


12 




C6 


04 41 52 


,6 


+25 38 00 

\ i-i^J *J U \J\J 


.2 


19.6 


149.8 


20976 


5 81 + 09 

i u i _!_ \j * vy ^y 


5 35 + 09 

*y*bj*y 1 v./ . vy y 


3.26 ± 





,22 


3 00 + 21 

* vy vy _!_ vy.^y. 


12 


LI 507 A 


Cl 


04 42 39. 


.3 


+29 43 56 

\ i-t j i .y .y vy 


.6 


16.0 


135.2 


18933 


3 86 + 07 

«-> . u vy _!_ vy .vy / 


3 55 + 06 

_!_ vy . vy vy 


2.24 ± 





.16 


2 07 + 1 5 

. vy / _!_ \y . j. »y 


12 




C2 


04 42 39, 


.6 


+29 41 06 


.2 


14.3 


128.0 


17916 


2.83 ± 0.06 


2.60 ± 0.06 


1.28 ± 





,13 


1.18 ± 0.12 


11 




C3 


04 43 03 


1 


+29 43 35 


.4 


8.7 


99.6 


13944 


1.75 ±0.06 


1.61 ±0.05 


0.88 ± 





.14 


0.81 ±0.13 


8 


L1582A 


Cl 


05 31 55 


.3 


+ 12 32 24 


.2 


5.1 


76.7 


30683 


1.16 + 0.04 


8.71 ±0.30 


0.74 ± 





.06 


5.59 ± 0.46 


11 




C2* 


05 32 03, 


1 


+ 12 30 43 


.2 


15.1 


131.4 


52554 


3.97 + 0.05 


29.86 ± 0.40 


1.48 ± 





,12 


11.16 ±0.91 


22 


IRAS05413 


Cl 


05 43 32, 


,7 


-01 07 27 


.6 


4.0 


67.3 


30278 


0.76 ± 0.05 


7.23 ± 0.49 


0.39 ± 





.05 


3.72 ±0.49 


6 




C2* 


05 43 46, 


.4 


-01 04 49 


.4 


17.2 


140.2 


63097 


4.08 ± 0.07 


18.24 ±6.32 


0.27 ± 





.04 


1.20 ±0.42 


32 




C3* 


05 43 51, 


.6 


-01 02 51 


.2 


10.6 


110.4 


49690 


3.11 ±0.04 


7.92 ±0.11 


0.39 ± 





,01 


0.99 ± 0.02 


94 




C4 


05 44 00, 


1 


-01 05 00 


.3 


3.5 


63.7 


28656 


0.59 ± 0.03 


5.64 ± 0.28 


0.43 ± 





.04 


4.09 ± 0.42 


8 


L1622A 


[Cl] 


05 54 21, 


,3 


+01 50 34 


.0 


1.8 


45.4 


20418 


0.52 ± 0.04 


4.90 ± 0.36 


0.37 ± 





.06 


3.51 ±0.55 


6 



Table 3. continued. 



Core 




a 


6 


A 








M tot 


^50% 


M 5Q % 


cbeam 
max 






J2000.0 


J2000.0 


arcmin 2 


arcsec 


AU 


Jy 


M G 


Jy 


M 






rP9*i 
L^z j 


OS SA 9Q 8 
W9 9t Z7.0 


-1-01 zL7 9Q S 
+U1 m / Z7. J 


9 

z.u 


ZL8 9 


91 6QA 
Z lOVt 


78 4- 07 
W. / it W.W / 


1 A6 4- 6^ 
/ .to it W.09 


07 4- 00 
W.W / it w.ww 


69 4- OS 
W.OZ it W.W9 


O 




p^ 


OS SA A 9 6 

Wj Jt- tZ.O 


-1-01 SzL ^1 zL 


91 ^ 

Zl .J> 


1 S6 


70999 
/ wzzz 


S AQ t 06 
J.t-7 it w.wo 


S9 97 4- SQ 
jz.z / it w. jy 


1 .TJ It W. 1 1 


1 ^ 69 t 1 9Q 

1 J .OZ It 1 .Z" 


98 
zo 




PzL 


OS SA S9 6 

WJ Jt JZ.O 


-1-01 S7 zLQ ^ 
tUI 9 / j 


1 s 1 

1 J.l 


1 J 1 . J 


SQ0Q7 


% OS 4- OS 
9.W9 it W.W9 


9Q 00 4- A6 
zy.ww it w.to 


1 10 4-011 

1 . 1 y ± w. 1 1 


11 ^9 4_ 1 00 
1 1 OZ it 1 .WW 


1 7 
1 / 




PS 


OS SS 09 Q 
WJ JJ \)L.y 


1 09 00 9zL (\ 
-v\jL WW Z'H-.O 


A 1 


7^ 9 
/ J.Z 


^9Q96 
9ZVZO 


77 4- 0^ 
v). 1 1 ± W.W9 


7 4_ ^0 
/ .9W it WOW 


A^ 4- OS 
W.tj it W.W9 


A OQ 4- AQ 
t.w" it \j.my 






T 1 8^ 


P1 


1 S SA 00 S 
1 J Dm WW. J 


09 S9 98 ^ 

— WZ jZ Zo.O 


8 1 
0. 1 


7O. J 


1 0SQ7 

1 W9t» / 


9 OQ 4- OS 
L.\jy ± W.W9 


1 1 Q 4- 0^ 
1 . 1 7 ± W.W9 


8S 4- 11 
W.o9 it W. 1 1 


A8 4- 06 
W.to it W.WO 


1 6 
lO 




P9 

lz 


1 S SA 08 6 
1 J Jt Wo.O 


09 S9 0^ % 
WZ JZ W9.9 


1 1 n 


1 7 


1 S^69 
1990Z 


6 7S 4- 06 
O. / 9 it W.WO 


% 8A 4- 0^ 
9.ot it W.W9 


1 Q9 4. 17 
1 .7Z ± U.l / 


1 OQ 4- 1 

1 .\jy ± w. 1 w 


^8 

JO 




P^ 


1 S SA OQ 1 
1 J Jt Wt». 1 


09 zLQ 1 7 8 
— wz <+y 1 / .0 


1 Q 


1 ZL7 7 


1 69S9 
IOZjZ 


A 8Q 4- 08 

0.07 it W.WO 


% Q9 4_ OS 
O.yL ± W.W9 


a A a _i_ 99 
j.tj it w.zz 


1 QS 4- 1 9 

1 .7J it W. IZ 


1 8 
1 




tpzli 


1 S SzL 1 9 9 
1J Jt YL.L 


09 S8 99 1 
— WZ JO zz. 1 


1 ^ 1 


1 99 ^ 
1ZZ. J 


1 ^AS8 
1 9tjo 


a OA 4- 07 
j.Wt it W.W / 


1 1% 4. OA 
1 . / j it w.wt 


1 77 -u 96 
1 . / / it W.ZO 


1 00 4- 1 s 

1 .WW it W. 19 








ps 


1 S SzL 1 % S 

1J Jt 1JO 


09 SS 9S S 

WZ JJ LD.D 


0. 1 


8^ 6 


Q1 QS 

y 1 yj 


1 97 4- OA 
1 .Z / it w.wt 


79 4- 09 
W. / Z it w.wz 


7S 4- 10 
W. / 9 it W. 1 W 


A^ 4- 06 
W.tj it W.WO 


1 
1 w 






1 S SzL 9S 9 

1J Jt ZJ.Z 


—09 SzL 98 zL 

— WZ Jt ZO.H- 


A S 


79 


7Q99 

/ 7ZZ 


70 t OA 
W. / W it w.wt 


AO t 09 
W.tW It w.wz 


A9 t 06 

W.tZ it W.WO 


9A t 0^ 

W.Zt it W.Wj 




t a^8 

l_-t JO 


pi 


18 1 A OzL 7 

1 O It- Wt. / 


07 08 99 8 
— w / wo zz.o 


1 J.Q 


1 ^ Q 

1 JJ.7 


^61 AS 
90 it j 


9 91 4- OS 
Z.Z 1 it W.W9 


7 S6 4- 1 7 
/ .90 it w. 1 / 


1 ^94-019 
1 OZ it W. IZ 


A S^ 4- A1 
t.jj it W.tl 


Q 

y 


T AQ9 


pi 

V^l 


18 1 S zL8 1 

1 O 1 J to. 1 


0^ zLS 1 1 f\ 
W9 HO 1 1 -O 


99 ft 
ZZ.O 


1 fil 1 


A^AQA 
t jtyt 


S ^6 4- OS 
9. 90 it W.W9 


18 ^S t 17 
1 0.99 t w. 1 / 


9 1^4-0 10 
Z. 1 j it w. 1 w 


7^1 4_ ^S 
/ . j 1 ± W.99 


^A 
9t 


PR1 88 




1 Q 90 1 S 1 


-1-1 1 ^S SO ^ 

Til J J JU.J 


8 7 
0. / 


QQ Q 

yy.y 


9QQS6 


1.JJ it W.Wj 


9 fS\ 4_ 07 

Z.O J it W.W / 


90 t 09 
W.ZW It w.wz 


^Q 4- 01 

WO" it W.Wl 


AO 
tw 


T 67^ 7 

L*U / j- I 


pi * 


1 Q 91 ^S S 

IV Zl JJO 


+ 11 zi vy . 1 


1 a n 


1 Of, 8 
IZO.o 


^809Q 


9 -21 1 OS 
Z.j 1 it W.W9 


8 ^8 4-0 91 
o.9o it W.Z1 


90 4- OS 
W.ZW it W.W9 


79 4- 09 
W. / Z it w.wz 


1 7 
1 / 




rP9i 


1 Q 9 1 zLzL 

I7 Zl tt.W 


1 7 S8 1 

+ 11 1 / jO. 1 


0. / 


1 no n 


9QQQ1 

Z777 1 


1 OQ 4- OS 

1 .\jy ± W.W9 


A 69 4- 90 
t.OZ it w.zw 


A6 4- OQ 
W.tO it W.W;? 


1 Q6 4- ^8 

1 .yyj ± woo 


7 
1 






1 Q 91 S^ 7 

1" Zl J J. 1 


-1-11 90 OS Q 
tii zw w j . y 


9 7 


zL 


1 669Q 


^0 t 09 
W.JW it w.wz 


1 9S t 1 1 
1 .ZJ It w. 1 1 


1 S t 0^ 
W. 1 J it W.Wj 


6^ t 1 A 

W.O J it w. 1 1 


7 


T 67S 

L*D / j 


pi 


1 Q 9^ S^ 
17 Zj JJ.U 


iii 07 zL^ 

til u / ^+9.W 


1 c 

1 . J 


ZL9 1 
^fZ. 1 


1 9699 
1ZOZZ 


19-1-0 01 
W. 1Z it W.W1 


S1 4- OS 
W.9 1 it W.W9 


1 4- 09 
W.1W it W.WZ 


A9 4- 08 
w.tz it w.wo 


s 


t 1 1 no 

L 1 1 WW 


P1 * 


90 ^6 1 Q 6 
ZW jO I7.O 


iAQ S^ 19 ^ 
tOj 99 1Z.U 


1 8 % 


1 ZLzL f 


S78AQ 
9 / Ot7 


^ S^ 4- 06 
9.99 it W.WO 


19. W9 it W.V9 


0^1 4- OA 
WO 1 it w.wt 


1 n + or 
1 . 1 9 it w.wo 


^9 
9Z 


T 1 04 1 9 


rpi l 


90 ^7 09 ^ 

ZW J / WZO 


+ S7 zL1 91 7 

T J / *t 1 Z 1 . / 


1 7 


zlzt 1 

L r L r. 1 


1 76AA 
1 / ott 


96 t 09 
W.ZO It w.wz 


1 cn 4_ 1 s 

1 J I U.U 


1 6 t 09 
w. 10 it w.wz 


1 1 7 4_ 1 S 
1.1/ It W. 1 J 


u 




rP9i 


90 ^7 08 ^ 
ZW 9 / Wo. 9 


-i-S7 97 (\ 

TJ / J7 Z / .O 




88 7 
OO. / 


^SAQ1 
99tV 1 


1 1 9 4_ OA 
1 . 1Z it w.wt 


8 AA 4- ^0 
o.tt it WOW 


SA 4- 07 
WOt it W.W / 


A 07 4- S9 
t.W / it WOZ 








p^ 

V J 


90 ^7 1 ^ S 
ZW 9 / 19.9 


_uS7 zL7 SO S 
tj / *+ / 9W.9 


8 S 

O.J 


Q8 Q 
70.7 


^QS6A 
97»90t 


9 A9 4- OS 
Z.tZ it W.W9 


1 8 9^ 4- ^S 
lo. Zj it W.99 


1 1 S 4- OQ 

1 . 19 it w.wv 


8 67 4- 71 
o.O / it W. / 1 


1 7 
1 / 




PzL* 


90 ^7 91 
ZW J 1 z i .w 


+ S7 zLzL 08 "K 

T J / *+*+ WO. J 


99 1 


1 SQ 
1 jy.v 


6^61 9 

O JU 1Z 


8 97 4- 07 

.Z / It W.W / 


99 SQ t Q6 

ZZ.J7 It W."0 


QS t 09 
W."J it w.wz 


9 60 t 11 
Z.OW It w. 1 1 


79 
/ z 




tpsi 
IVjJ 


90 ^7 99 
ZW 9 / ZZ.W 


-i-S7 SO ^9 1 
tj / jW 9Z. 1 


zl 


79 8 
/ z.o 


9Q1 1 8 
Ly 1 1 


1 1 Q 4- OS 
1 . 1 y ± W.W9 


8 Q1 4- ^6 
o.T'i it woo 


7Q 4- 10 
w. / y ± w. 1 w 


S Q1 4- 7S 
j.7l it W. / 9 


Q 

y 


T 1 1 A8 

L 1 Ito 


pi * 


90 zLO zL8 
ZW tW to.W 


1 fJi 99 07 f\ 
to / zz w / .0 


7 9 
/ .z 


7U.0 


9QA^O 
Lym jw 


77 4- 0^ 
W. / / it W.W9 


% 19 4-0 1 % 
9. 1Z it W. 1 9 


^A 4- 06 
WOt it W.WO 


1 ^6 4- 9S 
1 OO it W.ZJ 


7 




P9 


90 zL1 08 7 
ZW 1 1 wo. / 


+67 90 01 zL 
tu / ZW W 1 . 1 


1 9 Q 


1Z1.U 


^QS1 A 
jyj it 


1 7A t OA 
1 . / 1 m w.wt 


8 6S t 91 

O.OJ It W.Z1 


87 t 08 

W.O / it W.WO 


A ^9 4- A9 
toz It w.tz 


1 
1 w 


T 1 1 SSF 

L 1 1 J jC 


rpi l 
IV ij 


90 zL^ SA 7 
ZW tj Jt. 1 


to / j / jl.j 


zl 8 


7zL ^ 


9A1 ^9 

Zt 1 jZ 


77 4- OA 
W. / / it w.wt 


% 89 4- 18 
j.oz ± u.10 


AO 4- 07 
W.tW it W.W / 


1 Q7 4_ ^A 
L.y / ± wot 


7 




rP9i 
IVZj 


90 zLzL 01 Q 

ZW tt Wl 


4-67 zLQ (\ 
to / my sj 


J.7 


OO. J 


91 609 
Z IOWZ 


S1 4-0 0^ 
W.9 1 it W.W9 


9 SS t 1 s 
ZO j it W. 19 


AO 4- OS 
W.tW it W.W9 


1 Q6 4- 96 
1 .yo ± w.zo 


s 

J 




rp^i 
IVjj 


90 zLzL OQ zL 

ZW *+*+ W".t 


■A-f\l zL1 ^8 

tu / *+l JO.W 


8 


JU.J 


QQ01 

77UI 


10+001 
w. 1 W it W.Wl 


S 1 +0 07 

WO 1 it W.W / 


OS t 01 
W.Wj it W.Wl 


9^ 4- 07 

W.Z J it W.W / 


u 


T 1 089P 9 


rpi l 


90 zLQ SS 
zw my 99. w 


4-60 1 S zLQ 7 
tOW 19 *+y . / 


if, n 
zo.u 


1 79 7 

1 / Z. / 


6Q079 
Ot»W / Z 


f. A1. _i_ 1 1 

O.t j it w. 1 1 


A8 ^S 4- 80 
to.99 it W.oW 


1 Q8 4- ^A 
1 .70 it wot 


1 A QO 4- 9 SS 
it.yw it zoo 


1 ^ 
1 9 


T 1 089P 


pi * 


90 S1 98 % 

L\J J 1 ZoO 


4-60 1 8 ^zL 7 
tow 1 jm. 1 




1 m 9 

1W / .z 


A98Q^ 
tZo"9 


9 ^A 4- 06 
Z.jt it W.WO 


1 9Q 4- 1 89 
1 \j.Ly ± 1 .oz 


A1 4-0 0^ 
W.tl it W.W9 


1 81 4- ^9 
1 .0 1 it WOZ 


91 
Zl 


T 1 089 A 


PI* 


90 S^ 1^7 

ZW J J 1 J. / 


4-60 1 zL zLS 
tow it tj.w 


zl zl 


71 
/ 1 .w 


98A09 
zotwz 


QQ 4- 0^ 
yj.yy m w.wj 


% f.Q 4_ AA 
J.O" It w.tt 


9A t 09 
W.Zt It w.wz 


8Q t 1 1 
W.O" it w. 1 1 


96 
zo 




P9 


90 S^ 97 Q 
ZW 99 L 1 .y 


4-60 1 zL ^ 8 
tow it jj.o 


1 fl 9 
1W.Z 


1 08 ^ 
lUo. J 


A^98 

t j jZo 


9 AH 4_ OA 
Z.t / it w.wt 


18 S6 4- 9Q 
1 000 it \j.Ly 


98 4- 00 
W.Zo it w.ww 


9 OQ 4- OA 
L.\jy ± w.wt 


6^ 
09 




V J 


90 S^ SO zL 
ZW 99 JW.t 


4-60 OQ AS 1 

tow \jy t j. 1 


1 f. 7 

10. / 


1 ^8 9 
1 Jo.Z 


SS97Q 

99Z / y 


% 1 A 4- OS 
9 . 1 1 it W.W9 


19 81 4- 79 
iz.o 1 it w. / z 


1 8 4- 01 

w.io it w.wi 


7S 4- OA 
W. / 9 it W.Wt 


9Q 

£,y 




CA 


90 S^ S1 Q 

ZW J J J1.7 


4-60 1 7 1 7 Q 
tow 1 / 1 / 


8 9 
o.z 


Q7 
y 1 .yj 


^8808 

JOOWO 


1 0% t OA 
1 .Zj It w.wt 


Q 9A 4- 98 

".Zt It W.ZO 


AO t OA 
W.tW It w.wt 


"K OS t ^0 
J.WJ it wow 


1 "K 
1 j 


T 1998 


PI* 


90 S7 1 ^ 9 
ZW 9 / 1 9.Z 


4-77 ^S zLzL 
t / / j j tt.w 


zn ^ 

T-J. J 


111 7 
zzz. / 


AASAS 
tt jt j 


n 9i 4- 1 

1 j .Zt it w. 1 w 


11 1 9 4_ 08 
1 1 . iz it w.wo 


S6 4- 01 
WOO it W.Wl 


A7 4- 1 
W.t / it W.Wl 


1 07 
1 W / 


R^rr»A8 
JjcriltO 


PI * 


90 SQ 1 ^ Q 
zw jy ij.y 


4-78 9^ OA A 
t / Zj wt.t 


I J.J 


1 ^ ZL 


9668Q 


A 1 7 4- OS 
t. 1 / ± W.W9 


9 06 4- 09 
Z.WO it w.wz 


A^ 4- 01 
W.tj it W.Wl 


9S 4- 00 
W.ZJ it W.WW 


1 SQ 

lJ7 


L1172A 


CI 


21 02 20.5 


+67 51 08.8 


3.8 


65.6 


18905 


0.82 ± 0.05 


3.21 ±0.18 


0.51 ± 0.11 


2.01 ± 0.41 


7 




C2* 


21 02 22.8 


+67 54 03.1 


12.4 


119.0 


34269 


3.04 + 0.06 


7.63 + 0.32 


0.77 + 0.12 


1.93 + 0.08 


21 


LI 177 


CI* 


21 17 39.0 


+68 17 30.3 


11.5 


114.8 


33063 


2.93 ± 0.05 


4.76 + 0.20 


0.47 + 0.01 


0.76 ± 0.03 


104 


L1021 


CI 


21 21 44.5 


+50 57 51.5 


8.8 


100.5 


25114 


1.02 + 0.03 


2.98 ± 0.09 


0.64 ± 0.08 


1.87 + 0.24 


8 


L1014 


CI* 


21 24 07.2 


+49 59 05.0 


4.7 


73.3 


18326 


0.85 ± 0.03 


1.93+0.07 


0.29 ± 0.03 


0.66 ± 0.06 


16 


LI 103-2 


CI 


21 42 14.2 


+56 43 34.1 


8.9 


101.0 


25254 


1.53 + 0.04 


4.48+0.13 


0.92 + 0.10 


2.70 ± 0.30 


8 



Table 3. continued. 



Core 




a 

J2000.0 


6 

J2000.0 


A 
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Table 4. Geometry data for subcores. See Table[3]for an explanation of the subcore designations. For every subcore the table lists the major and 
minor axis at the 50% peak intensity level and their ratio, the position angle of the major axis (east of north), the filling factor for an ellipse defined 
by the aforementioned axis lengths and orientations, and the effective radius at the 70% peak intensity level. 
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Table 4. continued. 
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1 u ± + 


n 

u. 


88 
,00 


_|_ 


n 


09 

.OJ 


11 + 04 
1 1 .u ± U.+ 


9908 + 118 

JZ70 ±110 


r A79 n 

\-XS 1 J- 1 


P1 * 


S6 -i- 64 
JO ± O+ 


1 689Q + 1 QOQO 


9A + 1 
JO ± 1W 


1 0784 + 9141 
lu / 0+ ±91 + 1 


1 


.6 


+ 


1 

1 . 


,9 


QO + 90 
y\J ± zw 


n 


88 
.00 




u 


96 
.ZO 


19 0+1 7 
1Z.U ± 1 . / 


9A09 + S1 9 
JOUJ ± J 1 J 




rP9i 


1 88 4- 19 

1 OO ± 1 J 


S64Q1 + 97QQ 
JO+71 ± 9/ yy 


1 AA + 1 S 
1UU ±19 


4Q780 + 4447 
+7 / Ou ± +++ / 


1 


.1 


+ 


n 

V). 


1 
. 1 


SO + 97 

Ju ± 9 / 


n 


j j 


_|_ 


n 


06 
-uO 


99 7 + S 1 
zz. / ± J . 1 


6806 + 1 S99 
OOUO ± 1JZJ 




rP9i 

L^9J 


110 4- 18 

1 lu X 1 O 


9988S + S9Q9 

JZOOJ ± JJ7J 


S9 + 1 9 

JZ ± 1 J 


1 S4S4 + 9899 
19+9+ ± 9099 


2 


.1 


+ 


n 

V). 


Zl 
.+ 


S6 + 1 

JO ± lu 


n 

u. 


<< 

, j j 


_|_ 


n 


1 1 
. 1 1 


17 7 + 90 
1 / . / ± j .7 


S9Q8 + 1 1 76 

JZ70 ±11/0 


I A7S 


pi 


191 4- 19 
191 ± 19 


9Q99S + 9QQS 
97ZJJ ± J77J 


AO + A 


1 8001 + 1 8A9 
lOuul ± 10UZ 


2 


.2 


± 


n 

V). 


.9 


1 9S + S 

1ZJ ± J 


n 


69 
.OZ 


_i_ 
± 




1 6 
. IO 


94 8 + S 9 
Z+.O ± J .z 


7446 + 1 S74 
/ ++O ± 1 J / + 


i 1 1 on 

1-.1 1UU 


PI * 


64 4- 1 1 
D^T ±11 


9S760 + 491 Q 
Z9 /OU ± +Z 1 y 


94 + 9 
9+ ± Z 


1 941 4 + AQ1 
19+1+ ± U7l 


1 


.9 


+ 


n 

V). 


,9 


91+6 
9 1 ± O 


n 


7 1 
, / 1 


_|_ 


n 


06 
.OO 


11 S + 4 
1 1 . J ± U.+ 


ZLAOS + 1 A1 
+OUJ ± 101 


I 1 041 9 


rpi l 


88 4- 8 
OO ± 


9S1 SQ + 91 OS 
99197 ± 91UJ 


S1 + 8 
9 1 ± 


90967 + 9994 

ZuZu / ± JZZ+ 


1 


.7 


+ 


n 

V). 


9 


8+14 
± 1+ 


n 


79 
, / z 


_|_ 


n 

u 


00 

.U7 


17 8 + 90 
1 / .0 ± J .u 


71 08 + 1 1 88 
/ 1 uo ± 1100 




rP9i 


1 66 + 1 9 
lOO ±19 


669SQ + S941 
OOZJ7 ± 9Z+1 


1 94 + 1 9 
19+ ± 1Z 


S9SQ7 + 4791 
9997 / ± + / 9 1 


1 


.2 


+ 


n 

V). 


9 

,z 


1 66 + 14 
lOO ±1 + 


n 


40 




n 


OS 
.09 


9S 8 + 9 Q 
ZJ.O ± Z.7 


1 0998 + 1 1 49 
1U9Z0 ± 1 1+Z 




P9 

^9 


1 64 4- S 
10+ ± 9 


6SSQ7 + 1 QOQ 

OJJ7 1 ± I7U7 


AQ + 1 
\}y ± lu 


97489 + 9Q1 4 
Z/+0Z ± J7I + 


2 


.4 


± 


n 

V). 


,9 


1 77 + S 

1 / / ± 9 


n 


78 
, / 


_i_ 
± 


n 


07 


97 Q + 9 Q 
Z / .7 ± 9.7 


1 1 1 70 + 1 S71 
1 1 1 /U ± 1 J / 1 




P4* 


ZlQ 4- 1 
+7 X 1 


1 Qzt99 + 9S7 
1 7+ZZ ± J J 1 


41 + 1 
+ 1 ± 1 


1 6S8S + 404 
luJOJ ± +W+ 


1 


.2 


+ 


n 

V). 


n 


S7 + 9 
J / ± J 


n 


08 
,70 


_|_ 


n 

U 


01 
.ul 


14 9 + 04 

l + .Z ± U.+ 


S6QS + 149 

JO7J ± i+J 




rpsi 

L^9J 


1 88 4- 1 6 

1 OO ± lu 


7S090 + 690S 
/ JUZU ± OZUJ 


AO + Q 

vy ± " 


97SSS + 9490 

Z / JJJ ± J4ZU 


2 


.7 


+ 


n 

V/. 


Zl 
.+ 


1 7S + 9 

1 / J ± J 


n 

u. 


79 
, / z 


_|_ 


n 


06 
.UO 


99 Q + 6 1 

9Z.7 ±0.1 


1 91 60 + 9496 
1J107 ± z+zo 


T 1 1 48 

1_*1 1+0 


pi * 


99zL 4- 99 
99+ ± ZZ 


1 0870Q + 7099 
IWO /U7 ± /UJZ 


71+19 
/ 1 ± 1Z 


99099 + 9Q4A 
ZJUZJ ± 97+0 


4 


.7 


+ 


n 

V). 


7 

. / 


4S + 9 
+9 ± 9 


n 


10 

,9 7 




n 
u 


06 
.UO 


1 S 9 + 6 9 
1J.Z ± o.z 


ZLQ9Q + 901 4 
+7Z7 ± ZU1 + 




P9 


901 +98 
Z7l ± zo 


Q44S8 + 0118 
7++ J ± "110 


89 + 1 1 
oz ±11 


96700 + 9696 
ZO / uu ±9090 


3 


.5 


+ 


n 

V). 


zl 

.+ 


SO + 9 

Ju ± Z 


n 


70 


_|_ 


n 


06 
-UO 


98 + 4 8 
JO .u ± +.0 


1 9944 + 1 S67 

1ZJ++ ± 1 Ju / 


I 1 1 SSF 


rpi l 


1 S8 + 14 
1J0 ± 1 + 


S1 97S + 4491 
919/9 ± ++9 1 


Q7 + 1 9 

y 1 ± 1Z 


91 9A9 + 97Q9 
9 1 909 ± 9 / 79 


1 


.6 


± 


n 

V). 


,9 


47 + Q 

+ / ± 7 


n 


4S 
,+9 


_i_ 
± 


n 


08 
.UO 


91 S + 4 7 
Zl .9 ± +. / 


607 1 + 1 S49 
07 / 1 ± 1 J+J 




rP9i 


1 79 + 90 
1 / z ± zu 


SS81 9 + 64Q1 
99o 1 9 ± U+7 1 


8A + 7 
OD ± / 


97QQA + 9941 
Z / 77O ± Z9+1 


2 


.0 


+ 


V). 


9 


80 + S 
Ou ± J 


n 


67 
,0 / 


_|_ 


n 


08 
.UO 


99 S + S 8 
JJ.J ± J.O 


1 08Q1 + 1 88S 
1U071 ± 100J 




rP9i 


80 + 1 


9AOQS + 9999 


98 + 8 
90 ± 


1 9908 + 9S91 

IZZuO ± ZJZ1 


2 


.1 


+ 


n 

V/. 


zl 

.+ 


99 + 99 

JJ ± ZJ 


n 


99 
, j j 


_|_ 


n 


1 8 
. 1 


8 4 + 94 
.+ ± J .+ 


9799 +1191 
Z / ZZ It 1 1 z 1 


I 1 089P 9 


rpi l 


998 + 9S 

Zjo ± Z9 


0^1 A9 4_ 1 01 49 
"9109 ± 1W1+Z 


1 90 + 19 
IZu ±19 


4808S + S914 
+0u09 ± 9Z1 + 


2 


.0 


+ 


n 

V). 


1 

. 1 


84 + 7 
0+ ± / 


n 


SQ 

.9 7 




u 


OS 
.U J 


96 + 6 9 
JO.U ± O.J 


1 441 9 + 9S9Q 

1++1J ± ZJJ7 


I 1 089P 


P1 * 


S8 + S 

JO ± 9 


99989 + 9070 
ZJZo9 ± ZU/U 


49 + 4 
+z ± + 


1 A778 + 1 44S 
1 / / ± 1 ++9 


1 


.4 


+ 


n 


1 

, 1 


1 67 + 1 1 
IO / ± 11 


n 


OS 

.79 




n 


OS 
.U9 


1 A + Q 

IO. U ± U.7 


641 A + 9AS 
O+IO ± JOJ 


T 1 089 A 


P1* 


48 + S 
+0 ± 9 


1 0977 + 91 78 

1 7Z / / ± Z 1 / 


90 + 1 
Ju ± 1 


19141 + S80 

1Z1 + 1 ± JOU 


1 


.6 


± 


n 

V). 


9 

,z 


1 99 + 4 
1 zz ± + 


n 

u. 


Q8 
.70 


_|_ 


n 


09 
.UJ 


11 7 + 04 

11./ It U.+ 


4679 + 140 
+u /z ± 1+7 




P9 
*^Z 


90 + 1 
JU ± 1 


1 90AA + 900 


9S + 1 
Z9 ± 1 


1 0070 + 904 
1 \J\J l\J± ZW+ 


1 


.2 


+ 


n 

V). 




99 + 6 
9Z ± O 


n 


00 

77 


_l_ 


n 


01 
.Ul 


Q 7 + 9 
7. / ± U.J 


986Q +111 
J0O7 ±111 




pa* 

V J 


99 + 9 
99 ± 9 


1 99^^ +1010 
19999 ± lulu 


98 + 1 
Zo ± 1 


1 1 9A0 + 9Q4 
1 1 90u ± 97+ 


1 


.2 


+ 


n 

V). 


1 

, 1 


97 + 1 
9 / ± 1U 


1 

1 . 


00 
An) 




n 


09 
.UZ 


Q 7 + 4 
7. / ± U.+ 


986Q + 1 66 
J0O7 ± 1 OO 




P4 


1 89 + 1 S 
ioz ± 19 


79704 + S891 
/ Z / U+ ± 9 9 1 


99 + 1 9 
9Z ± 1Z 


1 9891 + 4800 
1Z091 ± +0uu 


5 


.7 


+ 


1 

1 . 


1 

. 1 


46 + 4 
+0 ± + 


n 


8S 
,09 


_l_ 


u 


. 1 J 


91 A + 9 7 
Zl .0 ± J . / 


86S9 + 1 4Q4 
OOJZ ± 1+7+ 


I 1 998 


P1 * 
^1 


99 + 1 
9Z ± 1 


A9S0 + 1 79 
099U ± 1 / Z 


97 + 1 
Z / ± 1 


S4Q7 +119 
9+7 / ± 1 1 9 


1 


.2 


± 


n 


n 


19+7 
1 9 ± / 


1 


00 
An) 


_l_ 


n 


1 
.Ul 


10 1 +00 
1U.1 ± u.u 


901 9 + 
ZU 1 J ± u 


Rf>rn48 

DClll+O 


P1* 


9S + 
ZJ ± u 


ZlQA.1 + 99 
+"U 1 ± JZ 


99 + 1 
zz> ± 1 


ZLSOO + 1 94 

tJw ± 1Z+ 


1 

1 


1 


+ 


n 

u 


n 

A) 


1 80 + 9S 
1 Ou ± ZJ 


1 
1 


00 


_l_ 


n 


01 
.Ul 


8 4 + 01 
.+ ± u. 1 


1 67S + 1 Q 

10 / J ± 17 


I 1 179 A 
1^1 1 /Zr\ 


P1 
V^l 


1 07 + 10 
1U 1 ± 1U 


9087S + 9QQS 
jUo / 9 ± Zjyyj 


1 00 + 11 
1 uu ±11 


987A9 + 9907 
ZO / OZ ± 99u / 


1 


.1 


+ 


n 

V). 


9 
,z 


1 1 Q + 48 
1 17 ± +0 


n 

u. 


70 




n 


1 9 
. 1 9 


91 9 + A 9 
Z 1 . J ± o.z 


61 9S + 1 78S 
OIZJ ± 1 / OJ 




P9* 
*^z 


QO + Q 


9SQAQ + 9A44 
Z9"U" ± ZO++ 


74 + Q 
/+ ± y 


91 9A1 + 9479 
Z1901 ± Z+ / 9 


1 


.2 


+ 


n 

V). 


1 

. 1 


70 + 4Q 

/ u ± +7 


n 


89 
,oz 




n 


07 
.u / 


91 8+9 1 
z 1 .0 ± z. 1 


6981 + AOS 
ozo 1 ± ouo 


LI 177 


CI* 


28 ±0 


8142 ± 122 


26 ± 


7457 ± 103 


1 


.1 


+ 





.0 


105 ±5 


0. 


,96 


± 





.02 


9.7 ± 0.2 


2786 ± 70 


L1021 


CI 


216+12 


54043 ± 3097 


115 ± 12 


28869 ± 2922 


1 


.9 


+ 


0. 


.2 


42 ±6 


0. 


.69 


± 





.07 


35.6 ±5.8 


8899 ± 1443 


L1014 


CI* 


62 ± 10 


15519 ±2484 


49 ±3 


12255 ±851 


1 


.3 


+ 


0. 


,2 


160 ± 18 


0. 


,92 


± 





,07 


18.2 ± 1.5 


4551 ±378 


LI 103-2 


CI 


273 ±7 


68328 ± 1785 


81 ±8 


20275 ± 2006 


3 


.4 


+ 


0. 


,4 


56 ±2 


0. 


.69 


± 





.08 


36.0 ±5.0 


8995 ± 1255 


L1251A 


[CI] 


173 ± 15 


51784 ±4605 


98 ±24 


29293 ±7108 


1 


.8 


+ 


0. 


,5 


103 ±11 





,54 


± 





,12 


26.7 ± 7.2 


8013 ±2149 




C2* 


147 ± 10 


44246 ±3114 


56 ±2 


16868 ± 703 


2 


.6 


± 


0. 


,2 


81 ±3 


0. 


,59 


± 





.08 


19.7 ± 1.1 


5923 ± 332 



2 

> 

CD 

o 
•§ 



C/3 

►a 



Cl. 

in 
3 
B- 

Q 
g 

a. 

n 

o 



Table 4. continued. 



Core 




arcsec 


a 

AU 


arcsec 


b 

AU 


a /b 


P.A. 
deg 


/ 


arcsec 


AU 




C3* 


37 ± 1 


10998 ± 292 


35+1 


10459 ±312 


1.1 ±0.0 


71 ±21 


1.00 ±0.01 


11.8 ±0.2 


3554 ± 69 




C4* 


158 ± 12 


47291 + 3503 


41 ±3 


12350 ± 759 


3.8 ±0.3 


135 ±3 


0.80 ± 0.05 


18.1 ± 1.3 


5429 ± 376 


LI 197 


CI 


200 ± 72 


59945 ±21604 


59 ± 11 


17669 ± 3436 


3.4 ±0.6 


127 ±4 


0.48 ± 0.06 


25.0 ±2.9 


7493 ± 861 



Table 5. Data on dust emission peaks. Bracketed and starred peak designations indicate peaks with uncertain properties due to artifacts and peaks 
with associated YSOs or VeLLOs, respectively. For every peak the table lists the peak position, the associated subcore, if applicable, the peak flux 
density per 11" beam (after smoothing to 20" resolution) and the related column density and extinction, the flux density and radius for an aperture 
of 4200 AU radius, and the flux densities for the c2d standard apertures between 20" and 120" diameter. 



Core 


a 


6 


Asso. FjSf 


AT(H 2 ) 


A v 


F4200AU 


M4200AU 


F 20" 


F 40" 


F 80" 


F 120" 




J2000.0 


J2000.0 


mjy beam 1 


10 22 cm" 2 


mag 


mjy 


M 


mJy 


mJy 


mjy 


mjy 



L1355 


rpn 


02 52 41.2 


+68 57 01 


rcn 


11 ± 2 




TP21 


02 52 58 1 


+68 53 17 


rc2i 


13 + 2 




TP31 


02 53 18.9 


+68 57 25.5 


rc3i 


14 + 2 




TP41 


02 53 40.3 


+68 53 55.5 




11 ± 2 


L1521B-2 


PI 


04 23 30.8 


+26 38 07.0 


CI 


15 ± 1 




P2 


04 23 37.1 


+26 40 23.5 


rc3i 


9 ± 1 




rP3i 


04 23 42.3 


+26 41 26.5 


rc3i 


6 ± 1 


L1521F 


PI* 


04 28 39.4 


+26 51 34.0 


CI* 


66 ± 2 


L1521-2 


PI 


04 29 26.9 


+27 00 53.5 


CI* 


8 ± 1 




P2 


04 29 37.4 


+26 58 05 5 


C2 


19+1 


L1524-4 


PI 


04 30 07 5 


+24 25 48.0 


CI 


19+1 




rp21 


04 30 15.9 


+24 24 20.5 


rc2i 


11 + 2 




TP31 


04 30 22.3 


+24 27 26.0 


rc3i 


10 ± 2 


B18-1 


PI 


04 31 45.4 


+24 32 44.5 


CI 


13 ± 1 




P2 


04 31 47.4 


+24 31 55 5 


CI 


15+1 




P3 


04 31 55.6 


+24 32 55.0 


CI 


49 ± 1 




P4 


04 32 08 9 

V 7 1 ^' jl— UU>7 


+24 30 21.0 




9 + 2 




rP5*i 


04 32 15.4 


+24 29 00.5 


rc2*i 


> 65 ± 


B18-2 


pi 


04 32 30.0 


+24 50 17.0 


Cl 


9 + 2 




P2 


04 32 41.3 


+24 51 30 5 


C2 


14+1 


TMC-2 


PI 


04 32 43.8 


+24 23 14.5 


Cl 


45 ± 1 




P2 


04 32 45 3 


+24 24 31.5 


Cl 


41 ± 1 


B18-4 


PI* 


04 35 35.2 


+24 08 17.5 


Cl* 


71 ± 2 




P2 


04 35 37 7 


+24 09 17 


Cl* 


71 ± 2 


B18-5 


PI 


04 35 51.0 


+24 09 27.5 


Cl 


25 ± 1 


TMC-1C 


PI 


04 41 21.8 


+26 02 32.5 


Cl 


27 ± 1 




P2 


04 41 32.2 


+26 03 42.5 


C2 


36+1 




P3* 


04 41 38.7 


+25 56 25.0 




12 ± 2 




P4 


04 41 39.2 


+26 00 12.5 


C3 


39+1 


TMC-1 


PI 


04 40 59 5 


+25 53 54 5 


Cl 


21 ± 2 




P2* 


04 41 12.8 


+25 46 33.5 


C2* 


52 ± 2 




P3 


04 4 1 13 8 


+25 49 49 5 


C3 


36 + 2 




P4 


04 41 37.2 


+25 43 38.5 


C4 


30 ±2 




P5 


04 41 43.9 


+25 41 22.0 


C5 


31 ± 3 




P6 


04 41 51.9 


+25 38 09.5 


C6 


24 ± 2 


L1507A 


PI 


04 42 30.5 


+29 39 22.5 


C2 


10 ± 2 




P2 


04 42 38.3 


+29 43 48.5 


Cl 


20 + 2 




P3 


04 42 38.9 


+29 45 54.5 


Cl 


15 ±2 




P4 


04 42 41.0 


+29 41 07.5 


C2 


18 ±2 




P5 


04 43 02.0 


+29 43 41.5 


C3 


17 ±2 


L1582A 


PI 


05 31 49.1 


+ 12 32 00.0 




14 ±2 




P2 


05 31 56.1 


+ 12 32 24.5 


Cl 


17 ±2 




P3 


05 32 01.1 


+ 12 30 32.5 


C2 


30 ± 1 




P4 


05 32 05.4 


+ 12 31 14.5 


C2 


22 ± 1 



0.72 


+ 





16 


0.87 


+ 





.16 


0.95 


+ 


0. 


.12 


0.76 


+ 


0. 


.16 


1.00 


+ 


0. 


.09 


0.62 


+ 


0. 


07 


0.43 


+ 


0. 


.09 


0.51 


+ 


0. 


.10 


1.26 


+ 


0. 


.09 


1.30 


+ 


0. 


.09 


0.76 


+ 


0. 


.10 


0.70 


+ 


0. 


.16 


0.89 


+ 


0. 


07 


1.01 


+ 


0. 


.07 


3.30 


+ 


0. 


08 


0.59 


+ 


0. 


.12 


0.58 


+ 


0. 


11 


0.91 


+ 


0. 


.09 


2.98 


+ 


0. 


.07 


2.73 


+ 


0. 


.07 



7.6 ± 1.7 

9.3 ± 1.7 
10.1 ± 1.3 
8.1 ± 1.7 

10.6 ± 0.9 
6.6 ±0.8 
4.6 ± 1.0 

5.4 ± 1.0 
13.4 + 0.9 
13.9+1.0 
8.1 ± 1.1 
7.4+1.7 

9.5 + 0.8 

10.7 ± 0.8 
35.1+0.8 
6.3+ 1.3 

6.2+1.2 
9.7+1.0 
31.7 + 0.7 
29.1+0.7 



4 


.75 


± 


0. 


11 


50.5 + 1.2 


1 


.68 


± 


0. 


08 


17.9 + 0.8 


1 


.78 


± 


0. 


08 


18.9 + 0.9 


2 


.42 


± 


0. 


09 


25.8 + 0.9 


2 


.63 


± 


0. 


09 


28.0+ 1.0 


1 


.40 


± 


0. 


15 


14.9+1.6 


2 


.39 


± 


0. 


11 


25.4+1.2 


2 


.00 


± 


0. 


16 


21.3 + 1.7 


2 


.08 


± 


0. 


18 


22.2+1.9 


1 


.63 


± 


0. 


14 


17.3 ± 1.5 





.69 


± 


0. 


12 


7.4+1.3 


1 


.35 


± 


0. 


11 


14.4+1.1 


1 


.03 


± 


0. 


13 


11.0+1.4 


1 


.20 


± 


0. 


11 


12.8+1.1 


1 


.14 


± 


0. 


13 


12.2+1.4 





.92 


± 


0. 


14 


9.7+1.5 


1 


.16 


± 


0. 


11 


12.3 + 1.1 


2 


.01 


± 


0. 


09 


21.3 ± 0.9 


1 


.44 


± 


0. 


09 


15.3 + 0.9 



83 + 14 
102 ± 14 
115+ 11 
84+15 
259+ 11 
153 + 10 
78+12 
972 ± 17 
135 + 12 
331 ± 11 
343 ± 12 
204+ 13 
115 + 21 
219 + 9 
274 + 9 
808 ± 10 
61 ± 16 
> 360 ± 50 
116+15 
223 ± 12 
835 + 9 
849 + 9 
974 ± 19 
1282 ± 14 
462 ± 10 
525 ± 10 
705 + 11 
71 ± 17 
794 ± 12 
406 ± 19 
684+ 18 
668 ± 15 
569 ± 20 
605 ± 23 
460+ 18 
116+15 
369 ± 14 
173 ± 17 
308 ± 14 
289 ± 17 
43 + 6 
48 + 5 
82 + 4 
58 + 4 



0.01 
0.01 
0.01 
0.01 
0.01 
0.01 



0.16 + 0.03 
0.19 + 0.03 
0.22 ± 0.02 
0.16 + 0.03 
0.24 ± 0.01 

0.14: 
0.07: 
0.77: 
0.12: 
0.30: 
0.32: 

0.19 + 0.01 
0.11 + 0.02 
0.20 ± 0.01 
0.25 ± 0.01 
0.74 ± 0.01 
0.06 ± 0.01 
> 0.18 ±0.02 
0.11 ±0.01 
0.21 ± 0.01 
0.77 ± 0.01 
0.78 ± 0.01 
0.55 ± 0.03 
1.18 ±0.01 
0.43 ± 0.01 
0.48 ± 0.01 
0.65 ± 0.01 
0.04 ± 0.00 
0.73 ± 0.01 
0.37 ± 0.02 
0.38 ± 0.01 
0.62 ± 0.01 
0.52 ± 0.02 
0.56 ± 0.02 
0.42 ± 0.02 
0.11 ±0.01 
0.34 ± 0.01 
0.16 ±0.02 
0.28 ± 0.01 
0.27 ± 0.02 
0.32 ± 0.05 
0.36 ± 0.03 
0.61 ± 0.03 
0.44 ± 0.03 



31 + 7 

36 + 7 

40 + 5 
28 + 7 
38 + 4 
24 + 3 
17 + 4 
173 + 6 
21 + 4 
50 + 4 
50 + 4 

30 + 4 

31 + 7 

37 + 3 

41 ±3 
124 ±3 
23 ±5 
> 204 ± 
26 ±5 
35 ±4 
111 ± 3 
101 ±3 
195 ±6 
172 ±5 

64 ±3 
63 ± 3 
89 ±4 

38 ±5 
97 ±4 
55 ±6 
143 ±6 
91 ±5 
80 ±7 
83 ±7 

65 ±6 
31 ± 5 
50 ± 5 

43 ±5 
48 ±4 
47 ±6 
40 ±6 

44 ±4 
76 ±4 
54 ±4 



14 



79 ± 14 
91 ± 13 

110 ± 10 
84 ± 14 
120 ± 8 
73 ±6 
46 ±8 
537 ± 11 
59 ±8 
161 ±7 
164 ± 8 
98 ±9 
61 ± 14 
104 ±6 
128 ±6 
421 ±7 
53 ± 10 
> 293 
56 ± 10 
112+8 
402 + 6 
379+6 
524 ± 13 
643 + 9 
217 ± 6 
245 ±7 
328 ±7 
58 ± 11 
356 ±8 
174 ± 12 
366 ± 12 
307 ± 10 
256 ± 13 
277 ± 15 
199 ± 12 
66 ± 10 
168 ±9 

111 ± 11 

140 ±9 
146 ± 11 
61 ± 12 

141 ±9 
260 ±7 
188 ±7 



30 



198 ± 28 
241 ± 27 
334 ±21 
246 ± 29 
438 ± 16 
256 ± 13 
113 ± 16 
1476 ± 23 
216 ± 17 
556 ± 15 
605 ± 16 
324 ± 18 
206 ± 29 
349 ± 12 
446 ± 13 
1200 ± 13 
90 ±21 
> 480 ± 85 
172 ± 20 
369 ± 16 
1388 ± 11 
1488 + 11 
1557 ± 27 
2068 ± 19 
809 ± 13 
885 ± 14 
1197 ± 15 
58 ± 23 
1350 ± 16 
613 ±25 
1053 ± 25 
1147 ±20 
999 ± 27 
986 ± 30 
766 ± 24 
143 ±21 
653 ± 19 
296 ± 23 
507 ± 18 
470 ± 23 
125 ± 25 
493 ± 18 
837 ± 15 
632 ± 15 



337 ± 44 
528 ± 42 
695 ±31 
449 ± 45 
832 ± 23 
473 ± 19 
249 ± 25 
2471 ± 34 
416 ±25 
1124 ±22 
1190 ±25 
634 ± 27 
468 ± 45 
704 ± 19 
945 ± 19 
1992 ± 20 
210 ±32 
> 584 ± 473 
334 ± 30 
641 ± 24 
2694 ± 17 
3130 + 17 
3126 ±41 
3841 ± 29 
1566 ± 19 
1740 ± 21 
2306 ± 22 
76 ± 35 
2659 ± 24 
1164 ±38 
1900 ± 37 
2246 ± 29 
1960 ±41 
1997 ± 45 
1553 ± 36 
256 ±31 
1335 ± 28 
576 ± 34 
1030 ± 27 
863 ± 34 
302 ± 40 
904 ± 28 
1477 ± 23 
1283 ± 23 



Table 5. continued. 



O 



Core 




a 


S 


Asso. 


^beam 


N(K 2 ) 


Ay 


F 4200 AU 


-^4200 AU 


^20" 


^40" 


^80" 


^120" 








J2000.0 


J2000.0 




mJy beam" 1 


10 22 cm" 2 


mag 


mJy 


M Q 


mJy 


mJy 


mJy 


mJy 




ip a cnC/1 i o 
1KAoUj41 j 


n 1 
rl 


nc A-y A 

Uj 4j jZ.4 


m no 00 ^ 
— Ul U/ Z5.J 


L.1 


lj + J 


n oc __. n 1 
U.V5 ± U. 1 / 


1 A C , 1 Q 

1U.J ± 1.5 


in _i_ a 
J / ± 


n qa _i_ n nA 
U.JO ± U.UO 


QQ _i_ O 
JO ± / 


1 no _i_ 1 a 
1U5 ± 14 


oon _i_ on 

Z5y ± zy 


CAO _i_ A A 

jy 1 ± 44 






do* 
rz 


AC An AC. A 

Uj 43 40.4 


m n/1 /in ^ 
— Ul U4 4U.J 


LZ 


1CI _i_ 1 
3V ± 1 






nn _i_ 1 
yy ± j 


n a c _i_ n 1 /i 
U.4j ± U.14 


1 nn _i_ q 
1UV ± J 


OAQ _i_ O 

Zoj ± / 


on/i _i_ 1 a 
/U4 ± 14 


1 000 _i_ 00 
1Z55 ± ZZ 






Fj 


ac /i c 1 

Uj 43 j l.o 


m no ^0 n 
— Ul UZ JZ.U 


Lj 


1 O/I _i_ 1 

1Z4 ± 1 






111 _i_ 1 
j j / ± j 


n nc _i_ n n 1 
U.Vj ± U.U1 


QAA _i_ A 
jOV ± 4 


An/i _i_ 
OV4 ± / 


1 O/I O _i_ 1 c 

1Z4Z ± lj 


MIA _i_ OQ 

1 / j4 ± Zj 






F4 


nc /ii An n 
Uj 43 OU.U 


m n^ m ^ 
— Ul Uj ULj 


L4 


11 + 1 


n oq _i_ n 1 n 
U. / J + U. 1U 


_i_ 1 1 
/. / ± 1.1 


OC _i_ A 

Zj ± 4 


n O/i _i_ n n/i 
U.Z4 ± U.U4 


OA 1 A 

jU ± 4 


nn _i_ q 

yu ± 


ono 1 1 

zy / ± 1 


COC _i_ OO 

J / J ± Z / 




T 1 AOO A 

L 1 oZZA 


rr> 1 ~\ 


ac c/i on o 
Uj j4 ZU.Z 


1 m cn /io n 
+U1 jU 4Z.U 


1 1 1 


1 n _i_ 1 
Ly ± j 


1 00 _i_ n 1 n 

l.Zo + u.iy 


1 1 A _i_ O 1 

1 J.O ± Z.l 


A C _i_ O 

4j ± / 


n /io _i_ n no 
U.4Z ± U.U/ 


cn _i_ 
JU ± 5 


1 QA _i_ 1 A 

1 jy ± 10 


QOA _i_ QC 

J lo ± JJ 


AC/1 _i_ CO 

0j4 ± j / 






tdoi 
|FZJ 


n^ m 11 n 
Uj j4 jl.U 


■ n 1 /io n/i n 

+U1 4V U4.U 


LLZ J 


1 C 1 Q 
lj + J 


1 no -1- n 1 
1 .UZ ± U. 1 / 


1 n q -1- 1 c 

iu.y + 1.5 


"31-1-0 

Jl ± / 


n Qn -1- n aa 
U.jU ± U.UO 


Q< -1- 1 
JJ ± / 


1 QQ -1- 1/1 
1 Jj ± 14 


QAQ -1- Q 1 
jUo ± Jl 


'nQA 4. /I 
jjO ± 45 






dq 

rj 


nc ca /in q 

Uj j4 4U. j 


, ni ^/i ns < 
+U1 J4 Uo.j 


Lj 


11 -J- 1 
J / + 1 


a^ .j- n no 
Z.4j + U.Uo 


OA 1 -1- n Q 
ZO. 1 + \j.y 


1 1 a 

ol ± J 


n 00 -1- n no. 
U. / / ± U.Uj 


QA -1- A 

yO ± 4 


OQQ -1- O 

zyy ± / 


QA^ j- 1 ^ 

o0j ± 1j 


1 A A/1 -1- OO 
1404 ± ZZ 






D/1 


n^ ^o n 

Uj j4 jZ.U 


+U1 J / Jl.j 


r^/i 
L4 


1 Q 1 1 
Ly + 1 


1 og -1- n no 
l.Zo + U.Uo 


1 1 1 _i_ n q 
1 j. / + U.5 


/IO -1- 1 

4Z ± J 


a /in -1- n nQ 
U.4U ± U.Uj 


/1Q -1- 1 
4o ± J 


1 00 j- 1 


CA/1 -1- 1 Q 
J04 ± I J 


I 1 1 q , OA 

I I iy ± ZU 







D'n 
rJ 


n^ ^/i <"3 /i 
Uj j4 jj.4 


, ni co ac a 
+U1 J5 4j.U 


L4 


OO _i_ 1 
ZZ ± 1 


1 A c , n no 
1.4j + U.U5 


1 c -1- n a 
l j.j + U.y 


Al -1- 3 
4/ ± J 


n /j c 1 n a q 
U.4J ± U.Uj 


J4 ± 4 


1 QO -1- O 
lo / ± / 


AQQ _i_ 1 A 

4yo ±14 


QOQ _i_ OO 

o5y ± zz 






DA 

ro 


AC CC A/1 

Uj jj U4.j 


1 no nn /in c 
+UZ UU 4U. J 


r^c 
Lj 


1 _l_ 1 

1Z ± 1 


a 00 1 n 1 n 
U.5j ± U.1U 


q q _i_ 1 n 
5.5 ± l.U 


OA _i_ A 

ZO ± 4 


A OC _i_ A A/1 

U.Zj ± U.U4 


jl ± 4 


1 nc _i_ q 
IUj ± 


OA 1 _i_ 1 O 

zyi ± 1 / 


CO/1 _i_ OA 

jZ4 ± Zo 


g 


T 1 QQ 

LIoj 


D1 
t! 


1 < ^/i nn q 
lj j4 UU.o 


no ^0 O/i c, 
— UZ jZ Z4.j 


LI 


0< j- 

Zj ± Z 


1 aq -1- n 1 n 
i.oy + u.iu 


10 n± 1 1 

lo.U ±1.1 


AQO -1- 1 O 
J / ± 1 / 


A QA -1- A A 1 
U.JO ± U.U1 


A/1 -1- A 
04 ± 4 


ono -1- o 

zu / ± y 


/Uj ± lo 


1 QOA -1- OO 
1 J 10 ± Z / 


ffman 




tdot 


1 c c/i nA /i 
lj j4 U0.4 


no cn /in n 
— UZ J5* 4V.U 


LC4J 


1 /i _i_ 
14 ± Z 


n nc _i_ n 1 /i 
U.5*J + U.14 


1 n 1 _i_ 1 c 
1U. 1 + l.j 


1 C/1 _i_ O/I 

1 j4 ± Z4 


a aa _i_ n ni 
U.UV ± U.U1 


A< _i_ A 

4j ± O 


AA _i_ 1 O 

OO ± 1Z 


1 OO 1 OC 

1 / / ± Zj 


/1AA _i_ QO 

4oU ± 35 




DQ 


1 c C/i r\o o 
1 J j4 Uo. 1 


no co 10 c. 
— UZ jZ jo.j 


LZ 


C A 1 1 
JO + 1 


q oa 1 n nn 
J. /O ± U.UV 


/in 1 1 1 n 
4U. 1 + l.U 


1 Q 1 A 1 1 C 
1 jlO ± 1 J 


n oc 1 n n 1 
U. / J ± U.U1 


1 A O 1 A 

L4-Z ± 4 


A OQ 1 O 

4oj ± O 


1 A OA 1 1 A 
14Z0 ± lo 


O A 1 A 1 O/I 

Z4iy ± Z4 


n et a 




n 1 
F4 


1 c c/i nn /i 
1j j4 UV.4 


no A Q CO A 

— UZ 45 Jo.U 


f^i 
Lj 


1A _i_ O 

j4 ± Z 


on _i_ n 1 
z.zy + U.1Z 


O/i 1 _i_ 1 1 
Z4.J ± l.j 


1 no 1 _i_ on 
1UZ1 ± ZU 


n cq _i_ n m 
U.Jo ± U.U1 


OC _i_ c 

5j ± J 


q 1 n _i_ in 
j 1 U ± 1 U 


1 1 Q Q _i_ Ol 
1 1 JJ ± Zl 


OO 1 O _i_ QO 

ZZ1 / ± jZ 






|C C/l 1 A A 

1j j4 1U.4 


no cn ^c n 
— UZ jV jj.U 


rr^/in 
LL4J 


to 1 

15 + Z 


1 1 n _i_ n 1 /i 

i.iy ± u.14 


1 O _i_ 1 c 
1Z. / + l.j 


111 _i_ OQ 

jj / ± Zj 


n 1 n _i_ n n 1 
u.iy ± U.U1 


/in _i_ a 
4V ± 


1 O/I _i_ 10 

1Z4 ± 1Z 


QO 1 _i_ OC 

j / 1 ± Z j 


onn _i_ qo 

ouy ± J / 






rO 


1 c C/i 1 q A 

lj J4 lj.0 


no ^ /i /i n 
— UZ jj 44. U 


Lj 


1 A j. 
10 ± Z 


1 no -1- n 1 n 
l.Uo + U.IU 


11 ^ -L. 1 1 

1 l.j ± 1.1 


QQ/1 _i_ 1 O 
J5*4 ± 1 / 


n 00 -1- A Al 
U.ZZ ± U.U1 


/IO -1- /I 
4Z ± 4 


1 O/l _i_ Q 

1Z4 ± y 


All -l. 1 Q 
4JZ ± lo 


QC/1 _i_ OO 
oj4 ±11 


2 




["DOT 

[F/J 


1 C C/1 OC Q 

lj j4 Zj.o 


no C/i 'sn c 
— UZ j4 jU.j 


LL6J 


I 1 _i_ 

I I ± Z 


n oa 1 n n 
U. / o + U.l J 


Q 1 __. 1 A 

5.1 + 1.4 


1 Q/1 _i_ Ol 

1V4 ± Zl 


nii 1 n ni 
U. 11 ± U.U1 


QQ _i_ C 

JJ ± J 


A/1 _i_ 1 1 

04 ± 1 1 


O 1 Q 1 OQ 

Ziy ± Zj 


Al A __. Q C 
4/1 ± JJ 


> 
S 


T /1QO 

L4j5 


D1 
t! 


1 C 1 A AA Q 
15 14 UO.o 


no nc ai < 
— U/ U5 4j.j 


1 

LI 


1 n j_ 1 
1U ± 1 


n aq -1- n nc 
U.0V + U.Uo 


/.4 + U.5 


^0 -j- < 
J / ± J 


A OA -1- A AO 
U.ZU ± U.UZ 


OQ -1- 1 

Zo ± j 


00 -1- O 
5Z ± / 


QOQ -1- 1 Q 

jZj ± Ij 


A/1/1 -1- OA 
044 ± ZU 


CD 


T /1Q0 
L4yZ 


D1 
t! 


1 Q 1 C AQ A 

15 iJ 4o.4 


na /i^ /ig ^ 
— Uj 4j 4o.j 


LI 


OQ -j— 1 
Ly ± 1 


1 Q/i 1 n nA 
1.5*4 + U.UO 


on -1- n a 
ZU. / + U.O 


1 A/1 _i_ A 
104 ± 4 


A ^A -1- A A1 
U.JO ± U.U1 


OO -1- O 

/z ± z 


OAQ -1- ^ 
Z0j ± J 


QAQ -l. 1 n 

y4y ± iu 


1 OOQ -1- 1 < 

1 /Zj ± 1 j 


O 

S 




tdot 
|FZJ 


IO 1 C CI A 

15 1j j1.4 


AQ C,A AA C 

— Uj jU UU.j 




in 1 t 
1U ± Z 


n AO 1 A K 

U.OV + U.l j 


O A _i_ 1 c 

/.4 + l.j 


Al _i_ 1 n 
4/ ± 1U 


A 1 A _i_ A A Q 

U. Io ± U.Uj 


Q 1 __. A 
Jl ± O 


CQ __. 1 Q 
JJ ± lj 


110 , OQ 

1 lo ± Zo 


1 OA _i_ C 1 

1 ly ± jl 


r^R 1 qc 
Lr>lo5 


D 1 * 
rl 


1 on 1^1 
Ly ZU 1 J. 1 


,11 QC /IO C 

+ 11 J J 4Z. J 


r^i * 
Ll 


JJ + 1 






1/10 -j- ^ 
14Z ± J 


A OQ -1- A A 1 
U.Zo ± U.U1 


1 nn -1- 
1UU ± z 


OIO 1 C 

Z15 ± J 


c,/in 4. 1 n 
J4U ± 1U 


Q/IO -1- 1 < 

o4Z ± 1 j 




T 1 

LO / J- / 


D 1 * 

rl 


1 n 1 "2 /i a 
iy ZI 34.0 


1 1 1 01 in ? 

+ 1 1 zi iy.j 


Ll 


OA _i_ O 

ZO + z 






1 nA _i_ a 
1U0 ± 


a qa _i_ n m 
U.jy ± U.U1 


OQ _i_ /I 

/ j ± 4 


1 0/i _i_ n 
1 /4 ± y 


11 A _i_ 1 Q 

j /4 ± lo 


A OO _i_ OA 

45Z ± Zo 


TO 




rZ 


lOOl QQ Q 

Ly LL jy.y 


,11 on 11 ^ 
+ 11 ZU ZJ. j 


Ll 


1 _i_ 1 
1Z ± 1 


n cn -1- n no 
U.5U + U.UV 


q < _i_ 1 n 

5.J + l.U 


^O -1- ^ 

jZ ± j 


A OO -1- A AO 
U.ZZ ± U.UZ 


1A -j- /I 
J4 ± 4 


cn -1- q 
5U ± 


O^O -j- 1 A 
ZjZ ±10 


^AQ -j- 1A 
jOJ ± Z4 


O 




rPQ 1 


1001 /ion 
Ly Zl 4Z.U 


_i.11 12 /in n 
+ 11 15 4V.U 


|LZJ 


I 1 1 1 

II + 1 


n 00 -1- n 1 n 
U. /Z ± U. 1U 


0-u.11 

/ . / ± 1.1 


IO _i_ A 
J5* ± 


A 1 A -1- A AO 
U.10 ± U.UZ 


QA -1- /I 
JU ± 4 


A'n -1- Q 

Oj ± 


ZUj ± 1 / 


/IOQ -1- OA 
4Zo ± ZO 


Spitz 




|F4J 


1Q01 CQ /I 

Ly ZL JJ.4 


,11 on nA n 

+ 11 ZU UO.U 


LljJ 


t A 1 O 

1U ± z 


n aq -1- n 1 n 
U.Oo + U.IU 


/.j ± 1.1 


Al _i_ A 
4j ± 


A 1 Q _i_ A AO 
U. lo ± U.UZ 


Qn -1- a 
JU ± 4 


OA -1- Q 

10 ± y 


1 CA -1- 1 Q 
150 ± lo 


Q/1Q _i_ OO 
J4J ± Z/ 


T AO^ 


D1 
rl 


1 O OQ ^A n 
iy Zj JO.U 


, 1 1 (\h on ^ 
+ 1 1 U / ZU. j 


1 

Ll 


^ _l_ 1 

J + 1 


n ^n -1- n no 
U.jU + U.U / 


1 _i_ c\ n 
j.Z + u. / 


1 Q -1- /I 
IV ± 4 


A AQ -1- A AO 
U.Uo ± U.UZ 


11 i 1 

1 J ± J 


QC 1 A 

JJ ± 


QO -1- 1 1 

yz ±11 


1 </1 -i- IO 
1 J4 ± 1 / 





t 1 i nn 
LI 1UU 


D 1 * 

rl 


On QA 1 Q 

ZU jo iy.z 


+oj jj 10. j 


r 1 1 * 
Ll 


/1/1 J_ 1 
44 ± 1 






1 Qn __. /i 
1 jU ± 4 


A /I Q _i_ A A/1 

U.4o ± U.U4 


1 K __. /I 

IZj ± 4 


OQ/1 __. O 

zy4 ± / 


OA1 _i_ 1 c 

/ol ± lj 


1 OAC __. OO 

IZOj ± ZZ 



0. 


L1U41-Z 


TD1 1 

|F1J 


on qo ao i 
zU j / UZ. / 


■ Co /i 1 qo n 
+j / 41 jZ.U 


LL1J 


11-j-O 
1 1 ± Z 


n 00 -1- n 1 n 
U. /Z + U.IU 


0-u.11 

/. / + 1.1 


Q/1 1 A 

J4 ± 4 


A OA -1- A AQ 
U.ZO ± U.Uj 


QQ 1 /I 

JJ ± 4 


OA -1- Q 

10 ± y 


OOQ -1- 1 Q 
ZZJ ± lo 


A AO -1- OA 
4U/ ± Z0 






["DOT 


on qo no o 
ZU j / Uo.y 


+ J / JV JO.J 


|L2J 


1 1 __. 

lj ± Z 


n on __. n 1 n 
u.yu ± U.IU 


Q A _i_ 1 1 

y.o ±1.1 


/1 1 __. /I 
41 ± 4 


A Q 1 _i_ A A Q 

U.Jl ± U.Uj 


QQ __. A 

jo ± 4 


1 nn __. q 

iuu ± y 


QOO _i_ 1 O 
jZ / ± Ll 


A1 A __. OA 
01 ± Zo 


mall Clo 




Fj 


on ^o n q 
zU j / lz.j 


1 co /io en n 
+J / 4/ JU.U 


Lj 


on _i_ 
zy ± Z 


1 no __. n 1 1 

i.y / + u.11 


on o __. 1 1 

zu.y ±1.1 


Q/1 1 C 

o4 ± j 


A AQ _i_ A A/1 
U.Oj ± U.U4 


on _i_ a 
ly ± 4 


oqc _i_ n 

zjj ± y 


OC 1 _i_ 1 Q 

/Jl ± lo 


1 QQO _i_ OO 

1 jj / ± Ll 




D/1 * 

F4 


on in 01 1 
ZU J / Z 1 . 1 


1 co /i /i 1 1 n 
+J / 44 1 J.U 


r^/i * 
L4 


1 1 /i _i_ 1 
114+1 






iia _i_ /i 
jj4 ± 4 


A AQ _i_ A AQ 

U.Vj ± U.Uj 


Qnn _i_ /i 
jUV ± 4 


ocn _i_ 
5 jU ± 


OQQA _i_ 1 

ZUoO ± 1 / 


11 A A _i_ OC 

jjo4 ± Zj 




|FjJ 


on ^o 00 n 
ZU J / ZZ.U 


1 co cn 1 n 
+J / jU 1 /.U 


LljJ 


15 ± Z 


1 00 1 n n 
l.ZZ ± U.l J 


1 1 n __. 1 /i 
1 J.U ±1.4 


/IO _i_ A 

4/ ± 


A QC _i_ A A/1 

U.jj ± U.U4 


/I C _i_ c 

4j ± j 


1 -in _i_ 1 1 

i4y ± 11 


ACiA _i_ OQ 

4y4 ± Zj 


QAO _i_ QC 
OO / ± JJ 


c 
a. 


T 1 1 /I Q 
LI 14o 


D 1 * 
rl 


on /in ^A 
ZU 4U jO.Z 


_i_AO OO ^G ^ 
+0 / ZZ Jo.J 


r^i * 
Ll 


Q -1- 1 
V ± 1 






Q^ -1- /I 
JJ ± 4 


A 1 A -L. A AO 
U. 14 ± U.UZ 


OA -1- Q 
ZO ± J 


^Q -1- O 

jy ± / 


1 OO -1- 1 A 

1 1 z ± 14 


OQA -1- OA 

zyo ± ZU 


g 




DO 

rZ 


on /ii 1/1/1 
ZU 41 14.4 


_i_ao on A ^ ^ 
+0 / ZU 4j.J 


r^o 
LZ 


1 O -1- 1 
1Z ± 1 


n oo -1- n no. 
U. ly + U.Uo 


c /i -1- n c 
5.4 ± U.5 


A Q -1- /I 
45 ± 4 


A O/I -1- A AO 
U.Z4 ± U.UZ 


Qn -1- Q 
JU ± J 


QQ -1- A 

yj ± 


QQ^ -1- 1 Q 
JJJ ± iJ 


C/1Q 1 OA 
J4j ± ZU 


0- 


T 1 1 

LI 1 J jL 


TD1 1 
|F1J 


on /i/i nn ^ 
ZU 44 UU. J 


_i_ao in 11 ^ 
+0 / j / jj.j 


|L1J 


lO-i-O 

1Z ± z 


n C3 _i_ n 1 1 
U.5j + U.l 1 


C C -1- 1 1 
5.5 ± 1.1 


/1Q -1- A 
4V ± 


A O/I -1- A AQ 
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Table 5. continued. 
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Appendix A: Dust Emission Properties 

This appendix summarises formulas to convert the observed 
flux density of dust emission into column densities and masses. 
These are evaluated for the standard assumptions for dust prop- 
erties made by the c2d collaboration. The discussion presented 
here is intended to serve as a future reference for the c2d project. 



A.1. Molecular Weight 

As one usually wishes to express column densities in terms 
of particles per area, the molecular weight needs to be intro- 
duced into the equations. The molecular weight per hydrogen 
molecule, fin 2 , is defined via fin 2 m nN(H2) = M; here M is 
the total mass contained in a volume with //(H2) hydrogen 
molecules, and oih is the H-atom mass. It can be calculated 
from cosmic abundance ratios. For hydrogen, helium, and met- 
als the mass ratios are M(U)/M ~ 0.71, M(He)/M * 0.27, and 
M(Z)/M ~ 0.02 , respectively, where M = M(U) + M(He) + 
M(Z) (Cox 2000). As //(H) = 2//(H 2 ), M(H) = m H N(H), 



m 2 



Al 



2M 



2M 



m H N(H 2 ) m H //(H) M(H) 



2.8. 



(A.l) 



Note the difference to the mean molecular weight per free parti- 
cle, yU p , defined via ju p muN = Ai, where N ~ Af(H 2 )+Af(He) for 
gas with all H in molecules. As the helium contribution is domi- 
nated by 4 He, and each H2 molecule has a weight of two hydro- 
gen atoms, //(H 2 ) = M(H)/(2m H ) and //(He) = M(He)/(4m H ). 
Thus 



M 



m H [MH 2 ) + yV(He)] 
M 



m H [M(H)/(2m H ) + M(He)/(4m H )] 
M/M(H) 



[l/2 + M(He)/(4M(H))] 
2.37; 



(A.2) 

(A.3) 

(A.4) 
(A.5) 



the classical value of /i p = 2.33 holds for an abundance ratio 
N(H)/N(He) = 10 and a negligible admixture of metals. 

Both molecular weights are applied in different contexts. The 
mean molecular weight per free particle, // p , is e.g. used to eval- 
uate the thermal gas pressure, P - gk^T/fip, where T, g, and 
£b are the gas temperature, density, and Boltzmann's constant, 
respectively. The molecular mass per hydrogen molecule, yUn 2 , is 
needed below to derive particle column densities. 

A.2. Radiative Transfer 

A.2.1. Equation of Radiative Transfer 

The intensity emitted by a medium of temperature T and of op- 
tical depth t v at the frequency v is given by the equation of ra- 
diative transfer, which reads 



I v = B v (T)(\-z r ") 



(A.6) 



in the case of local thermal equilibrium. Here B v is the Planck 
function. For the optical depth, 



gds , 



(A.7) 



where k v is the (specific) absorption coefficient (i.e., per mass) or 
dust opacity. If most hydrogen is in molecules, the optical depth 
can be related to the column density of molecular hydrogen, 



N Hl - f « H , ds = f — - — ds = 1 f K v g ds , (A. 



8) 



and thus 



HH 2 mnK y 



(A.9) 



where «h 2 is the particle density of hydrogen molecules. This 
reads in a usable form 



Nh 2 = 2.14-10 25 cnT 2 r. 



0.01 cm 2 g- 1 



(A. 10) 



where the chosen numerical value of k v is characteristic for 
wavelength A w 1 mm. As typical H 2 column densities are be- 
low 10 23 cm 2 , t v is expected to be by far smaller that 1. Thermal 
dust emission in the (sub-)millimetre regime is therefore opti- 
cally thin. 

For optically thin conditions the equation of radiative trans- 
fer can be simplified: 



l V ~ By(T)Ty. 



(A. 11) 



As the optical depth is related to the column density (Eq. IA.10l . 
Eq. ( |A. 1 It relates the observed intensity to the column density. 

A.2.2. The Planck Function 
The Planck function reads 
2hv 3 1 



By(T) 



c 2 e hv/(k B T) _ I 



(A. 12) 



in which c is the speed of light and h is Planck's constant. In the 
Rayleigh- Jeans limit, hv <sc k%T, this simplifies to 



2v 2 

B v (T)=—k B T. 



However, the limiting condition, which reads 



A » 1 .44 mm 



(iok) 



(A. 13) 



(A. 14) 



in useful units, is under typical dust temperatures of about 10 K 
not fulfilled for observations at about 1 mm wavelength. The 
exact value of the Planck function is 

3 



B V (T) = 1.475 • 10~ 23 W nT 2 Hr'sr" 1 ( -^—) 

\GHz/ 

1 



e 0.0048(v/GHz)(r/10 K) _I _ ^ 

in useful units. 

A.3. Observed Quantities 

The received flux per beam is related to the intensity by 



(A.15) 



-/ 



IyPdQ, 



(A. 16) 



where P is the normalised power pattern of the telescope (i.e. 
max[P] = 1). Defining the beam solid angle as the integral 



PdQ, 



(A. 17) 
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one can derive a beam-averaged intensity, 



which reads 

(A.18) A^ H2 = 2.02-10 20 cm- 2 (e 1 - 439 W mm >"P : 



The beam solid angle can be conveniently approximated for tele- 
scopes with a beam profile similar to a Gaussian function, 



P(6) = e 



-e 2 /(28l) 



(A. 19) 



(where the angle 6 gives the distance from the beam center). The 
parameter #o is related to the half power beam width of the tele- 
scope, y i a 



00 = 



7 HPBW 



V8hi(2) 

For these idealisations the beam solid angle is 



/"•CO /-*c> 

Q A = In I P(6»)6>d6» = 27r I 
Jo Jo 



= 2nB\ 



e -o 2 l(2el) 0d0 



~ 41n(2) UHPBW - 
Using the conversion 

1 ai-csec =4.85- l(T 6 rad, 
one obtains 

Q A = 2.665 -lO-^srf^f 
\ arcsec / 



(A.20) 



(A.21) 

(A.22) 
(A.23) 



(A.25) 



The parameter #hpbw does not need to be the real telescope 
beam, but depends on the calibration of the data. To give an ex- 
ample, some software packages (e.g., MOPSIC) apply scaling 
factors to the data (i.e., Fjf am ) when spatially smoothing a map, 
so that the beam width to that the calibration refers is equal to 
the spatial resolution of the map after smoothing. 

For the idealisations made here the average intensity derived 
from Eq. dA. 18b is a good approximation to the actual intensity 
only if the source has an extension of the order of the main lobe 
of the telescope. Otherwise radiation received via the side lobes 
may have a significant contribution to F^' im , and the idealisa- 
tion of the telescope beam by Gaussian functions is too simple. 
Inclusion of the side lobes increases Q A , and thus reduces (/ v ). 
The nature of the average intensity derived in Eq. dA. 1 8b is there- 
fore comparable to the one of main beam brightness tempera- 
tures used in spectroscopic radio observations. Note that, how- 
ever, spatial filtering of bolometer systems effectively reduces 
sidelobes. Unfortunately, this reduction depends in details on the 
mapping pattern and the source source geometry. It can therefore 
not be included here. 



A.4. Derivation of Conversion Laws 
AAA. Flux per Beam and Column Density 



Equations (IA.9I IA. Ill IA. 18b relate optical depth, column den- 
sity, intensity, and the observed flux per beam with each other. 
Rearrangement yields 



^beam 



Q. A fi Hl mHK v B y (T) ' 



(A.26) 



7beam 



' \mm/ 

Id 



mJy beam 1 / \ 10 arcsec / 



(A.27) 



(0.01 cm 2 g- 1 
in useful units. 

A.4.2. Flux and Mass 



The mass is given by the integral of the column densities across 
the source, 



M - /j Hl m H J n h, dA. 

Substitution of Eqs. ( IA.9IIA.1U yields 

M = f I v dA. 

k v B v {T)J 



(A.28) 



(A.29) 



The surface element dA is related to the solid angle element d£2 
by dA = d 2 dQ, where d is the distance of the source. Thus 



M : 



7 v dQ 



d 2 F v 



k v B v (T)J k v B v (T) 7 
(A. 24) where F v — J I v dQ is the integrated flux. This reads 



(A.30) 



M 



0.12 M (e 



1.439(/l/mmr 1 (7'/10 K)" 







(0.01 cm 2 g' 
in useful units. 



, (— ) 

Jy / \ 100 pc / \mm/ 



(A.31) 



A.5. Dust Temperature and Opacity 

Dust temperatures of order 10 K have be en predicted fo r dense 
cores since some time. In the past years IZucconi et alj d200ll) 
estimated dust temperatures in the range 11 to 6 K for so- 
lar neighborhood c ores. The gas an d dust te mperature in thes e 
should be similar (Goldsmith" 2001; also see iGalli et al.ll2002h . 
Gas temper atures of order 10 K are indeed common for dense 
cores (e.g., iTafalla et al.l 120041) . Such temperatures even pre- 
vail in dense cores residing in molecular clouds for ming stel- 
lar clusters ( i .e., in Perseus; Rosolowskv et al. 2007). Finally, 
ISchnee et al.l d2007bl) recently mapped the dust temperature dis- 
tribution in TMC-1C using data partially also included in the 
present work. They derived dust temperatures of 1 3 to 5 K, sim- 
ilar to gas temperatures recently derived for LI 544 (Craps fet al.l 
|2007|) . These temperatures are lower than thos e typically in- 
ferred fro m IRAS (> 15 K; ISchnee et al.l l2005l) and ISO data 
(> 10 K; iLehtinen et al.l 120051) . possibly because these mis- 
sions primarily probed very extended emission because of large 
beams. We therefore here adopt a temperature of 10 K in absence 
of protostellar heating. 

Figure I A. 1 1 illustrates the uncertainty in mass estimates 
due to uncertainties in the temperature (based on Eq. IA.3 lb . 
For MAMBO this shows that for actual dust temperatures of 
12 to 8 K our assumption of a dust temperature of 10 K will 
lead to an error in the mass estimate of factors 0.7 to 1 .5, respec- 
tively. These increase to factors of 0.6 to 3 for temperatures of 
14 to 6 K. The uncertainty significantly decreases with increas- 
ing wavelength. 
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100 1000 
Wavelength [/Ltrn] 



Fig.A.l. The ratio of the actual mass to the one derived when 
assuming a dust temperature of 10 K. Solid lines indicate the 
ratio for various actual dust temperatures. Vertical labels indi- 
cate the observing wavelength of several existing and upcom- 
ing bolometer arrays (for completeness also including Spitzer 
MIPS). The related uncertainty in mass estimates decreases sig- 
nificantly with increasing wavelength. 



Dust opacities are p resently not well constr a ined. W e here adopt 
dust opacities from Ossenko pf & Henn ing (1994) that hold for 
dust with thin ice mantles coagulating for 10 5 yr at an H-density 
of 10 6 cirT 3 , respectively (keeping the product of time and den- 
sity constant) for 10 6 yr at an H2-density of 5 • 10 4 cirT 3 , which 
appear to be reasonable conditions. At the MAMBO wavelength 
of 1.2 mm their opacity of 0.01 cm 2 per gram of interstel- 
lar matter is in between va lues used an d suggested by other 
groups (from 0.005 cm 2 g~' llMotte et alJ ll998h to 0.02 cm 2 g" 1 
IlKriigel & Siebenmorgenll994H ). 

The opacities qu oted above are consisten t with observational 
opacity constraints, van der Tak et al. ( 1999) conclu de that only 
coagulated dust from Ossenkopf & Henningl d 19941) provides a 
match to their extensive dust and molecular line emission data 
set. They probed a massive young star with a heated envelope 
though. Anyway, Shirley et al. (submitted) derive the same for a 
less luminous YSO in B335. 
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Fig. B.l. Summary of mapping parameter spaces accessible by 
various techniques. The non-highlighted area at af SO urce/4ff < 1 
is accessible using conventional techniques. This is the realm 
in which any past and present bolometer arrays, including to- 
tal power devices like SCUBA II and LABOCA, operate. If one 
whishes to remove skynoise, and does not use an iterative ap- 
proach for this, this further limits the accessible parameter space 
to lisource/^an-ay < 1- Failure to comply with these limits re- 
sults in artifacts of various kind. The area with light highlight- 
ing is accessible when using our iterative restore algorithm, 
which by design also includes iterative skynoise filtering. This 
significantly increases the accessible mapping parameter space. 
However, as indicated by the arrow, the convergence speed de- 
creases when increasing d source /{ e ff. Maps occupying the area 
with strong highlighting cannot be reduced without artifacts. 



A.6. Standard c2d Conversion Factors 
A.6.1. Conversion from Dust Emission 

Table IA.1I summarises the c2d standard conversion fac- 
tors for masses and column densities from dust emis- 
sion. The wavelengt h -depe ndent dust opacities are from 
lOssenkopf & Henningl (|1994|) and hold for dust with thin ice 
mantles coagulating for 10 5 yr at an H-density of 10 6 cm -3 . We 
adopt a dust temperature of 10 K. This choice for the dust tem- 
perature and opacity are the standard assumptions made by the 
c2d collaboration. The values listed in Table lA. fl are thus thought 
to serve as a standard reference within the collaboration. 

A.6.2. Conversion to Extinction 

The c2d standard conversion factor between column densities 
and visual extinction, 

N Hl = 9.4 • 10 20 cnT 2 (Ay/mag) , (A.32) 



is taken from Bo hlin et~ai] dl978l) . They combined measure- 
ments of H2 and Hi from the Copernicus satellite for lightly red- 
dened stars to get 

([N m + 2N Hl ]/E(B - V)> = 5.8 ■ 10 21 cnT 2 mag" 1 . (A.33) 

For a standard total-to-selective extinction ratio Ry = Ay/E(B - 
V) = 3.1 this yields the above conversion factor. 

Appendix B: Data Reduction Details 

B.1. Mapping Parameter Space 

The equations on limitations of data reduction approaches, Eqs. 
([TJO, effectively span parameter spaces accessible by different 
reduction methods. We explore these limits in Fig. IB. II 

For brevity, we express the effective map size, i SC!m + { c hop, 
by ( e g in the following. The vertical limit at d somce /d dndy = 
1 directly follows from Eq. ©. In principle, it continues to 
<4ource/^eff — * 00 ■ However, Eq. (O anyway limits the realm of 
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Table A.l. Standard c2d conversion factors for masses and column densities from dust emission. For each bolometer camera 
employed by the c2d collaboration we list the effective wavelength, the half power beam width, the dust opacity at the effective 
wavelength (per gram of ISM), the conversion factors between intensity and column density, and between integrated flux density 
and mass (this for a distance of 100 pc), respectively. 



Camera 


A 


$HPBW 


Ky 




M/F v (d 




fim 


arc sec 


2 -1 

cm - g 


cm 2 (mjy beairT 1 ) -1 


M Jy- 1 


SHARC ii 


350 


8.5 


0.101 


7.13 • 10" 


0.031 


SCUBA 


450 


7 


0.0619 


1.42- 10 19 


0.041 




850 


15 


0.0182 


1.34 • 10 20 


0.18 


BOLOCAM 


1120 


31 


0.0114 
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unproblematic data reduction along this axis. To obtain an ex- 
pression for this limit, we divide Eq. (O by d. drr . dy and solve for 

^source /^array • At 

4source/4ff > 1, the region accessible to iterative 
techniques is limited by 

^source ^ ^scan ^chop ^array • (B . 1) 

Beyond this limit, no receiver of the bolometer array moves off 
source, and sky brightness information is impossible to recon- 
struct. Equation (Q~|) sets the indicated limits on the convergence 
speed. 

B.2. Illustration 

Figure |B72"1 summarizes the data reduction improvements due to 
our iterative method. To do so, we created artificial raw data 
(for the LI 103-2 region of our survey) towards a circular source 
of 500" diameter and spatially constant intensity. We included 
beam smoothing, and added correlated and uncorrelated noise at 
various amplitudes. In order to roughly match the observed de- 
crease of the noise level during skynoise filtering, the correlated- 
to-uncorrelated noise ratio is chosen to be 3. 

Figure IB .21 (a) presents a map obtained from noiseless 
raw data, using classical methods (i.e., plain reconstruction of 
the intensity distribution from the chopped bolometer data). 
Inspection reveals intensity depressions of > 30% within the 
source. Since such biases are additive, these depressions can off- 
set intensities of smaller objects located in the map center, if 
such are present. Within the inner area of 400" diameter, the 
mean intensity has dropped to 80% of the true value. 

Figure IB. 21 (b) is constructed like the map in panel (a), but 
was derived from raw data with superimposed artificial noise. 
We use noise-to-source intensity fractions of 90% and 30% 
for the correlated and uncorrelated noise. The simulated source 
thus corresponds to very faint extended features in our maps. 
Correspondingly, classical data reduction (i.e., reconstruction 
without skynoise filtering) gives a map full of artifacts (note the 
intensity gradient within the source) in which the signal-to-noise 
intensity ratio is of order 1 . In practice, such a map would not be 
usable. 

Figure |B~2l (c), however, shows that iterative map reconstruc- 
tion including skynoise filtering allows to reliably extract source 
structure from the data used for panel (b). Since most of the cor- 
related noise has been removed, the signal-to-noise intensity ra- 
tio is of order 3 here. Furthermore, the artifacts present in panels 
(a) and (b) are gone: the mean intensity within the inner 400" 
diameter is 101% for the map shown. 
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Fig. B.2. Demonstration of our bolometer data reduction ap- 
proach, using artificial maps centered on a source of 500" diam- 
eter. Panels (a) and (b) present maps reconstructed using classi- 
cal methods. Panel (c) displays a map derived using our method, 
which was derived by iterative processing of the data shown in 
panel (b). It is obvious that, even when using noisy raw data, our 
new approach yields maps with an artifact level better than those 
derived from noiseless data using classical methods (note de- 
pressions of > 30% in panel [a], and their absence in panel [c]). 
The arrows in panel (a) indicate the length (600") and mean 
orientation of the two co-added maps. Circles and their separa- 
tion show the array diameter and the chop throw. Labels indicate 
the presence of noise (a superposition of correlated and uncorre- 
lated noise at levels of 90% and 30%, respectively, of the model 
source peak intensity), or its absence. [See the alternative URL 
in the astro-ph comments, and the journal article, for higher 
resolution maps.] 
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